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Werden Kondensatoren über einen Metalldraht 
' entladen, so verdampft dieser explosionsartig unter 
starker Lichtemission, wenn die in den Kondensatoren 
gespeicherte Energie ein Mehrfaches der zur Ver- 
| dampfung des Metalls nötigen Energie beträgt. Diese 
sog. Drahtexplosionen wurden bereits unter verschie- 
‚denen Gesichtspunkten (Erzeugung sonnenähnlicher 
Absorptionsspektren, Lichtquelle hoher Intensität, 
Auftreten sehr hoher Temperaturen, Stoßwellen) 
untersucht und dabei der zeitliche Strom-, Spannungs- 
oder Lichtverlauf gemessen, Spektren aufgenommen 
| und zeitlich aufgelöste Bilder des Vorganges mit Kerr- 
| Zellenverschlüssen, Drehspiegelanordnungen oder 
| Trommelkameras gemacht!. Die Ergebnisse lassen 
sich leider zum großen Teil schlecht miteinander ver- 
gleichen, da die Versuchsbedingungen in wesentlichen 
| Punkten -(gespeicherte Energie, Stromstärke und 
Stromanstieg, Drahtmaterial) stark voneinander ab- 
wichen. 
In der vorliegenden Arbeit wird über Versuche 
| berichtet, bei denen Kupferdrähte mit 0,035 bis 
0,24 mm Durchmesser und 70 bis 125 mm Länge 
| durch Stromstöße mit Stromstärken bis zu 3,5 - 10% A 
| und Stromanstiegen bis etwa 10!! A/sec in Luft von 
| Atmosphärendruck verdampft wurden. Dabei wurde 
' der zeitliche Verlauf des Entladestromes, der Strom- 
änderung und der am Versuchsdraht liegenden Span- 
nung gemessen. Aus diesen Größen kann der zeitliche 
Verlauf der von den Drähten aufgenommenen Energie, 
des Drahtwiderstandes und des magnetischen Druckes 
_ berechnet werden. Die Ergebnisse erlauben Aussagen 
| über die Aufheizung des Drahtmaterials und die 
Unterbrechung der metallischen Leitung sowie die 
Größenordnung der im Inneren des Drahtes auftreten- 
_ den Drucke und Temperaturen. 


I. Versuchsanordnung 
1. Entladekreis 


Abb. 1 zeigt die Schaltung des Entladekreises und 
der Meßanordnung. Die beiden parallel geschalteten 
Kondensatoren hatten eine Gesamtkapazität von 

1,3 @F und wurden auf Spannungen zwischen 25 und 
55kV aufgeladen. Die Induktivität des Entlade- 
kreises betrug 0,55 uH. 


2. Meßanordnung 


! Die am Versuchsdraht auftretenden Spannungs- 
impulse, deren Spitzen zum Teil über 100 kV be- 


1 Das Literaturverzeichnis am Ende dieser Arbeit be- 
"schränkt sich im wesentlichen auf neuere Veröffentlichungen. 
- Eine ausführliche Zusammenstellung der bis 1955 erschienenen 
_ Arbeiten findet sich bei Con [5]. 
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trugen, wurden durch einen. ohmschen Spannungs- 
teiler von 2kQ auf einige kV untersetzt, über ein 
Koaxialkabel in die abgeschirmte Meßkabine über- 
tragen und dort durch einen zweiten ohmschen Span- 
nungsteiler von 3kQ auf die Eingangsspannung des 
Oszillographen untersetzt. Die Verschleifung eines 
Rechteckimpulses durch die beiden Spannungsteiler 
betrug 2 - 10°® sec. 


Abb.1. Schaltskizze des Entladekreises und der Meßanordnung. Ci,. 

Kondensatoren, F Funkenstrecke, D Versuchsdraht, K Koaxialkabel, 

O0 Oszillograph, M abgeschirmte Meßkabine, R,_, Spannungsteilerwider- 

stände, R,_, Abschlußwiderstände, R,,s Meßwiderstand für Strom bzw. 

Stromänderung, R, Aufladewiderstand, Z,, AN Spule zur Messung des Stromes 
bzw. der Stromänderung 


Der zeitliche Verlauf des Entladestromes und 
seiner Änderung wurde mit zwei Spulen gemessen, die 
jeweils einen der beiden Rückleiter umgaben und mit 
induktionsarmen Meßwiderständen abgeschlossen wa- 
ren. Die Induktivität Z, der Spulen und ihre Win- 
dungszahl N sowie die Meßwiderstände R,, waren so 
aufeinander abgestimmt, daß in dem einen Meßkreis 
die Bedingung 


(* = Spulenstrom) (1) 


erfüllt war und damit die am Meßwiderstand auf- 
tretende Spannung U, dem Entladestrom / propor- 
tional war — nämlich 

nr (2) 


m N 


und im anderen Meßkreis die umgekehrte Bedingung 


\ di 
Bm ti >L, 5: (3) 


so daß die Spannung U,, der Stromänderung propor- 
tional war — nämlich 


U, =. (4) 


KAyusec 0 


ce 4 
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Für beide Spulen mußte außerdem die Bedingung 
erfüllt sein, daß ihre Eigenfrequenz groß gegen die 
Frequenz der zu messenden Impulse war. 

Bei einigen Versuchen wurde zusätzlich der zeit- 
liche Verlauf der Lichtemission spektral unzerlegt mit 
einem Multiplier gemessen und, ebenso wie der zeit- 
liche Verlauf von Spannung, Strom und Stromände- 
rung, auf einen Oszillographen gegeben, wo er photo- 
graphiert wurde. 


II. Ergebnisse 


1. Zeitlicher Verlauf von Spannung, Strom 
und Stromänderung 


Bei den vorliegenden Untersuchungen wurden 
Kapazität und Induktivität des Entladekreises kon- 
stant gehalten, während die Kondensatorauflade- 
spannung und damit die zur Verfügung stehende 

Energie sowie Länge und 
w” Durchmesser der Kupfer- 


“ drähte in der aus Tabelle 1 
2 zu entnehmenden Weise 
Sera TEE ss variiert wurden. Da einer- 
Y2- seits die Schreibgeschwin- 
KM digkeiten der zur Verfügung 
stehenden Oszillographen 


nicht genau übereinstimm- 
ten, andererseits ein ge- 
nauer Vergleich der Oszillo- 
gramme von Drahtexplo- 
sionen, die unter gleichen 
äußeren Versuchsbedingun- 
gen nacheinander erfolgt 
waren, eine sehr gute Repro- 
duzierbarkeit des Kurven- 
verlaufes zeigte, wurde der 
zeitliche Verlauf von Span- 
nung, Strom und Strom- 
änderung jeweils nachein- 
ander mit dem gleichen 
Oszillographen aufgenommen. Um die Genauigkeit zu 
erhöhen, wurde jede Größe, z.B. der Spannungsver- 
lauf, mehrmals unter den gleichen Versuchsbedin- 
gungen gemessen. 

Als typisches Beispiel zeigt Abb. 2 den zeitlichen 
Verlauf der am Versuchsdraht liegenden Spannung, 
des Entladestromes und der Stromänderung während 
der Verdampfung eines Kupferdrahtes in zeitlicher 
Zuordnung. 

Wie der Stromverlauf zeigt, ist etwa nach einer 
Mikrosekunde der Draht soweit verdampft, daß die 
metallische Leitung unterbrochen ist. Der zu diesem 
Zeitpunkt einsetzende zweite Stromimpuls entspricht 
einer Gasentladung, in der die Kondensatoren ihre 
Restenergie abgeben. Im Spannungsverlauf tritt 
während der Stromunterbrechung eine deutliche 
Induktionsspitze auf, die ein Mehrfaches der Konden- 
satoraufladespannung betragen kann. 


7 


9 
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Abb.2. Zeitlicher Verlauf von 
Spannung am Draht (obere 
Kurve), Entladestrom (mittlere 
Kurve) und Stromänderung (un- 
tere Kurve) während der Ex- 
plosion eines Kupferdrahtes mit 
0,24mm Durchmesser und 124mm 
Länge. Kondensatoraufladespan- 
nung 50 kV 


2. Zeitlicher Verlauf 
der von den Drähten aufgenommenen Energie 


Es soll zunächst gezeigt werden, wie der zeitliche 
Verlauf der von den Drähten aufgenommene Energie 
aus dem gemessenen Verlauf des Entladestromes /(r), 
der Stromänderung dI/dr und der am Versuchsdraht 
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liegenden Spannung U(T)gem sowie meßbaren kon- 
SDR Größen berechnet werden kann. Dazu soll die, 
Zeitpunkt { betrachtet werden. Die von den Konden-) 
ide bis zu diesem Zeitpunkt Ü a: A 


ht; io ebsfeilndt wird, 


B.0) = [0,0 9,021. 


4 [tma 


ist, kann man Gl. (5) auch in der Form 


DO.) = Veto) 


t 
2 ao EEE Ki. zoljrena 


t 
genommen, während der Betrag R,- [ I(t)? dr i 


t z 
Entladekreis außerhalb des Versuchsdrahtes als Wirk- 
energie umgesetzt wurde und die Energie en (1)? 
im Magnetfeld steckt (C = Kapazität, L,, R, = In- 
duktivität und Widerstand des Entladekreises, L = 
Induktivität des Versuchsdrahtes, U, = Kondensator- 
spannung). Man erhält also für die vom Draht bis zur 
Zeit t aufgenommene Energie die Gleichung 


:— 


t 
— RS Ia)edr. 
to 


t t 
Bi) = U.) Amar — al ma 


6) 
Eine zweite Möglichkeit, #(t) zu berechnen, bietet 1 
Gleichung 5 
t 
= (Ile) ülddr 

h k 
Hier macht es einige Schwierigkeiten, den ea 
Anteil der am Ver suchsdraht liegenden Spannung U(t) 
zu bestimmen. Er kann erstens aus dem u 
Spannungsverlauf Ugem (T), der sich aus einem ohm- 
schen und einem induktiven Anteil zusammensetzt, 
nach der Gleichung 

i 


Dt) — Ugem (tl) — (LM): 


berechnet werden. Das Glied M -dI/dr berücksichtigti 
dabei die Spannung, die durch den Entladestrom /(t) 
in der vom Spannungsteiler und seinen Zuleitungen 
gebildeten Schleife induziert wird. Zweitens kann man - 
u) aus dem Verlauf von Strom und Stromänderung 
sowie der Kondensatoraufladespannung U,(t,) be- ; 
stimmen. Die Kondensatorspannung zu einem Zeit- 
punkt £ ist nämlich einerseits 
D, (6) = 


t 
U.lt)— 5 f Ie)az, () 
to 


andererseits ist sie 


T,() = (R+R,)Il)+(L+L,) 4 


a0) 
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1 

er'ohmsche Anteil der zum Zeitpunkt t am Versuchs- 
-aht liegenden Spannung ist also 
)-R- I) = U.) - 

| lar 11) 
| dt T=t k 


‘ D 

| Zune R,I()— (L,+L) 
ur praktischen Auswertung wurden die Spannungs-, 

rom-und Stromänderungskurven in 20 bis 30 Punkte 

ufgeteilt. Dann wurde der zeitliche Verlauf der am 

'ersuchsdraht liegenden ohmschen Spannung punkt- 

eise nach Gl. (8) und (11) berechnet, ebenso der zeit- 

che Verlauf der von den Drähten aufgenommenen 

nie nach GI. (6) und (7). Für die Zeit des starken 

omanstieges zu Beginn der Entladung konnte 

llerdings die ohmsche Spannung nicht berechnet 

rerden, da sie dort nur einen kleinen Bruchteil der 

ıduktiven Spannung beträgt, und somit die Gln. (8) , 
nd (11) eine Differenz von fast gleichgroßen Gliedern 

arstellen. Auch für die aufgenommene Energie 

önnen deswegen in diesem Bereich keine Werte an- 

egeben werden. Für die späteren Zeitpunkte stim- 

hen die nach den beiden GIn. (8) und (11) berechneten 

‚pannungswerte gut miteinander überein, ebenso die 

Inergiewerte, die nach Gl.(6) und (7) berechnet 

yurden. Da sowohl die ohmsche Spannung als auch 

lie aufgenommene Energie nach zwei Methoden be- 

echnet wurden, die verschiedene und voneinander 
habhängige Meßgrößen benutzen, ist eine gewisse 
Xontrolle der Ergebnisse gegeben. 

Abb. 3 zeigt für verschiedene Kombinationen der 
Tersuchsparameter die auf diese Weise berechneten 
Verte der vom Draht pro Masseneinheit aufgenom- 
menen Energie als Funktion der Zeit. Alle Kurven, 
lie zu Drähten mit gleichem Durchmesser gehören, 
ind in der gleichen Strichart gezeichnet. Zur Orien- 
jierung sind — ebenso wie in Abb. 5 — einige charak- 
eristische Energiewerte eingetragen: 


E, = (15 — T,) = 0,465 Wsee/mg 


(12) 


st die Energie, die zur Aufheizung des Kupfers von 
ler Zimmertemperatur 7, —=293° K auf die Schmelz- 
emperatur 7, = 1356° K nötig ist (cd = 0,1041 - 
4,185 Wsee/g - grad ist die mittlere spezifische Wärme 
wischen diesen beiden Temperaturen). Addiert man 
lazu die latente Schmelzwärme HE; —0,21 Ws/mg, so 
srhält man die Energie 


eg E, = E, + By = 0,675 Ws/mg (13) 
die zum Schmelzen des Kupfers nötig ist. 
# E,=E,+0c,(7,—T;) =1,31Wsimg (14) 


st die Energie, bei der die Siedetemperatur beim 
Druck von 1 Atm 7, —2570° K erreicht wird (6, = 
0,1248 - 4,185 Wsec/g - grad ist die mittlere spezifische 
Wärme zwischen T,, und 75) und 


E,=E,+ Ey = 6,1 Wsec/mg (15) 


die zur Verdampfung des Kupfers nötige Energie, 

wobei Ey =4,8 Wsee/mg die latente Verdampfungs- 

wärme ist. 

_ Aus Abb. 3 ist zu ersehen, daß den Drähten bei 

diesen stromstarken Entladungen bis zur Unterbre- 

er der metallischen Leitung ein Mehrfaches der 
nergie E, zugeführt werden kann, die beim Druck 


von 1 Atm zur völligen Verdampfung nötig ist. In 
Übereinstimmung mit den Ergebnissen einer früheren 
Arbeit [20] zeigt sich außerdem, daß die von den 
Drähten bis zur Unterbrechung der metallischen Lei- 
tung pro Masseneinheit aufgenommene Energie durch 
Erhöhung der Kondensatoraufladespannung gesteigert 
werden kann und daß sie, wenn die übrigen Parameter 
konstant bleiben, mit dem Drahtdurchmesser wächst. 
Letzteres gilt aber nur, solange die in den Konden- 
satoren gespeicherte Energie groß gegenüber der 
während des ersten Stromimpulses an den Draht ab- 
gegebenen Energie ist. 
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Abb. 3. Zeitlicher Verlauf der von den Drähten pro Masseneinheit 
aufgenommenen Energie (Bedeutung der Zeichen wie Abb. 5) 


3. Zeit- und Energieabhängigkeit des Drahtwiderstandes 


Der zeitliche Verlauf des Drahtwiderstandes wurde 
nach der Beziehung 
U) 
a — OR 
berechnet, wobei der ohmsche Anteil der am Draht 
liegenden Spannung U (f) wieder nach @l. (8) und (11) 
bestimmt wurde. Am einfachsten und genauesten 
kann der Drahtwiderstand für den Zeitpunkt des 
Strommaximums (die entsprechenden Werte finden 
sich in Tabelle 1, Spalte 6) berechnet werden, da er 
dort wegen dI/dr=0 durch den Quotienten aus ge- 
messener Spannung Ugem (t) und Entladestrom I(f) 
gegeben ist. In Abb. 4 ist für die verschiedenen Kom- 
binationen der Versuchsparameter der zeitliche Ver- 
lauf des Drahtwiderstandes, der zum besseren Ver- 
gleich jeweils auf den Kaltwiderstand R, bezogen ist, 
dargestellt. Auch hier und in Abb. 5 sind alle Kurven, 
die zu Drähten mit gleichem Durchmesser gehören, 
in derselben Strichart gezeichnet. Die beiden einge- 
zeichneten horizontalen Linien stellen das aus Tabellen 
entnommene Widerstandsverhältnis für den Schmelz- 
punkt dar. 

Aus den in Abb. 3 und 4 wiedergegebenen Ergeb- 
nissen erhält man die Energieabhängigkeit des Draht- 

4* 


(16) 
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widerstandes, die in Abb. 5 dargestellt ist. Sie stimmt 
bis zur Energie E, gut mit den aus bekannten Daten 
berechneten Werten überein, die durch die dick ein- 
gezeichnete Linie verbunden sind, und zeigt keine 
Abhängigkeit von den Versuchsparametern. Bei 
höheren Energien dagegen ist der Drahtwiderstand 
nicht mehr allein eine Funktion der aufgenommenen 
Energie. 

Das in Abb. 4 und 5 bei den Drähten mit 0,24 mm 
Durchmesser und dem Draht mit 0,10 mm Durch- 
messer, 70 mm Länge und 50 kV Aufladespannung zu 


De & 
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Verdampfung von Kupferdrähten 

| 
Magnetfeld eines starken Stromes /, der gleichmä 
über den Querschnitt eines Zylinders mit Radiu 
verteilt ist! und parallel zur Achse fließt, erzeugt im 
stationären Zustand den zur Achse hin gerichteten 


Druck 
PS] (m 


Setzt man für r den ursprünglichen Radius ein, so 
ergeben sich aus den gemessenen Stromkurven für d 
Zeitpunkt des Strommaximums magnetische Drucke 
zwischen 3-10%° und 27.10° Bar? (s. Tabelle 1% 
Spalte 5). Sie nehmen mit dem Drahtquerschnitt zur 


) 
| 


1 I[Weber]? 
r [em]? 


dyn 


cm? 


TT 
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Abb.4. Zeitlicher Verlauf des Drahtwiderstandes R/R, (Bedeutung der 
Zeichen wie Abb. 5) 


beobachtende Absinken des Widerstandes gegen Ende 
des ersten Stromimpulses rührt daher, daß dort eine 
Gasentladung einsetzt (vgl. Abb. 2). 


4. Magnetischer und gaskinetischer Druck 


Von verschiedenen Autoren [14], [17], [20] wurde 
bereits darauf hingewiesen, daß bei der Verdampfung 
von Drähten durch intensive Stromstöße das Draht- 
material durch den magnetischen Druck des vom 
Entladestrom erzeugten Feldes weit über die Siede- 
temperatur beim Druck von 1 Atm hinaus zusammen- 
gehalten wird. Auf Grund der bereits vorher beschrie- 
benen experimentellen Ergebnisse, nämlich des zeit- 
lichen Verlaufes von Entladestrom und aufgenommener 
Energie, kann näherungsweise berechnet werden, wie 
weit das flüssige Drahtmaterial über die Siedetempe- 
ratur beim Druck von 1 Atm hinaus aufgeheizt wird, 
ehe es stark zu expandieren beginnt und wann diese 
Expansion einsetzt. 

Zunächst sollen die magnetischen Drucke berech- 
net werden, die bei den Versuchen auftraten. Das 


Abb. 5. Drahtwiderstand R/R, in Abhängigkeit von der pro Masseneinheit | 
aufgenommenen Energie E N 


anwächst. In Wirklichkeit lagen die magnetischen 
Drucke etwas niedriger, da die Volumenausdeh- 
nung, die wahrscheinlich hauptsächlich in radialer 
Richtung erfolgte, berücksichtigt werden muß. Der 
höchsten für das Strommaximum berechneten Draht 
temperatur von rund 8000° K entspricht eine Ver- 
größerung des Volumens von 7,1 auf 12 cm®/Mol 


1 Die Voraussetzung der homogenen Stromverteilung war 
bei den hier beschriebenen Versuchen gut erfüllt. Denn bei 
den dieksten verwendeten Drähten mit einem Radius r 
0,12 mm war die höchste auftretende Frequenz f= 0,8 MHz. 
Nach [2] ist dann die Änderung des Widerstandes R, infolge 
der Stromverdrängung 
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wobei ö= (7 ef we n) "= 0,075 das Eindringmaß ist. Mit 
| 

steigendem Widerstand wird der Einfluß der Stromverdrän- 
gung dann sehr rasch noch geringer. 
® 1 Bar = 10% dyn/cm? kann im Rahmen der Genauigkeit 


dieser Abschätzungen gleich 1 Atm gesetzt werden. 
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(s. Abb. 6), so daß — unter Annahme einer rein radia- 
len Ausdehnung — in diesem Fall der berechnete 
magnetische Druck um einen Faktor 1,7 zu hoch ist. 
In den übrigen Fällen ist der Fehler geringer. 

Dem nach innen gerichteten magnetischen Druck 
wirkt der gaskinetische Druck entgegen. Seine Be- 
rechnung scheiterte bisher an der ungenügenden Kennt- 
nis des zeitlichen Energie- und Temperaturverlaufes 
und daran, daß es mit den bekannten thermischen 
Zustandsgleichungen nicht möglich war, das Verhalten 
von realen Flüssigkeiten und Gasen bei hohen Drucken 
richtig wiederzugeben. Dagegen scheint eine von 
Hımpan [3] angegebene thermische Zustandsgleichung 
auch bei hohen Drucken eine gute quantitative Über- 
einstimmung mit der Erfahrung zu liefern!. Sie hat 
die Form 


a 


J : 
u zT zT) 


} (Peer de 


wobei P, V und T die thermischen Zustandsgrößen 
sind, R die universelle Gaskonstante und a, b, ce, d vier 
Stoffkonstanten, mit denen die Gleichung den jeweili- 
gen Verhältnissen angepaßt werden kann. Zur Be- 
stimmung der Stoffkonstanten benötigt man vier von- 
einander unabhängige Gleichungen. Drei davon er- 
hält man wie bei der van der Waalschen Zustands- 
gleichung aus den thermischen Zustandsgrößen für 
den kritischen Punkt. Für die vierte Gleichung müs- 
sen die thermischen Zustandsgrößen eines zweiten 
Punktes bekannt sein. Da der kritische Punkt in 
einem Zustandsgebiet liegt, wo — bei konstant ge- 
haltener Temperatur kleine Druckänderungen 
bereits große Volumenänderungen hervorrufen, soll 
der zweite Punkt in einem Zustandsgebiet liegen, wo 
gerade das Umgekehrte der Fall ist. Damit soll er- 
reicht werden, daß die Zustandsgleichung in einem 
möglichst großen Bereich quantitativ gültig ist. Die 
kritischen Daten des Kupfers (P,—8750 Bar, AZ 
26 cm?/Mol, T,—=9300° K) wurden nach einem Ähn- 
lichkeitsgesetz [17] aus den kritischen Daten des 
Quecksilbers bestimmt. Als zweiter Punkt wurde die 
flüssige Phase des Siedepunktes bei Atmosphären- 
druck gewählt (P,=1,013 Bar, V,=10,2 em?/Mol, 
T,=2833°K)?. Damit ergab sich aus Gl. (18) für 
Kupfer die thermische Zustandsgleichung 


83,144: T 
— 72 9,891804 
34,766381 - 1010 
= (7 8,675056)(T-V+ 2.693918 - 105) 


12 


(19) 


Abb. 6 zeigt die nach dieser Zustandsgleichung be- 
rechneten Isothermen im P-V-Diagramm. 

Da der zeitliche Verlauf des magnetischen Druckes 
in guter Näherung berechnet und der Verlauf der 
Drahttemperatur aus der aufgenommenen Energie 
abgeschätzt werden kann, solange keine starke Ex- 
pansion des Flüssigkeitszylinders zu erwarten ist, 
können nun aus dem P-V-Diagramm die Volumina 


mit sehr gutem 
2 Um zu sehen, wie die 
zwar in einem bestimmten Gebiet liegen soll, aber nicht genau 


-795 festgelegt ist — die thermische Zustandsgleichung beeinflußt, 


Abb. 6.. P-V-Diagramm des Kupfers nach einer Zustandsgleichung von 


HIMPAN 


wurde auch einmal der Schmelzpunkt als zweiter Punkt ein- 
gesetzt. In Abb. 6 sind einige Isothermen, die nach der so 
gewonnenen Zustandsgleichung berechnet wurden, gestrichelt 


eingezeichnet. 


54 M. KEILHACKER: Explosionsartige Verdampfung von Kupferdrähten 


Tabelle 1. Versuchsergebnisse für die verschiedenen Kombinationen der Versuchsparameter 


an Strommaximum Beginn der starken Expansion 

SER s 

173 „u 

a © © 

ee ee E 2 E er) = E 8 
=) = | = 3 DM an 2 = Eh ® E 
En] = Z 8 = =! des BEE 50 P a {iM 

{3} re a = 17} ke} Beit= a8 u © =" I} © = 
Sl 8 & =| =| ep 35 ä =) = =| = 2 
5 Bee = = S 5A 2 =] = a =) Ei 
ee Be 8 = |g8 35 2 3, N S, = 
3 = {77} nn | = = a = A 
& | + F 

d DT; vs 
mm | mm | kV m?/Mol 


4,0 8000 13 


0,10 | 70 | 28 | 8 } | 3300 3500| 0,36 
1100 |40 | 0,22 | 8, { 42 ‚= 3200 3400| 0,26 | 4,25 |8000-8500| 14 
| 125 | 50 26 | 3800 4200| 0,30 | 4,25 8000-8500 14 
70 | 50 0,22 | 44 8500-9000, 17 
0,195 100 | 25 0,90 5 |8000-8500| 14 
5 


125 31,5 


8500— 9000| 14 


0,24 | 70128 
| 100 | 40 
| 125 50 
| 70.) 50 


Spalte 1 2 3 4 5 6 7 


bestimmt werden, auf die sich das flüssige Kupfer bei 
einem Gleichgewicht zwischen magnetischem und gas- 
kinetischem Druck nach und nach ausgedehnt hat. Es 
zeigt sich, daß bis zu Temperaturen, die je nach den 
Versuchsbedingungen zwischen 8000 und 9300° K 
(Tabelle 1, Spalte 11) liegen, die Volumenvergrößerung 
wegen des herrschenden Druckes von etwa 10* Atm 
nur gering war. Erst bei diesen Temperaturen er- 
reichte das Kupfer den Zustand, der im P-V-Dia- 
gramm durch den horizontalen Verlauf der entspre- 
chenden Isothermen gekennzeichnet ist und begann 
stark zu expandieren. Je stärker der magnetische 
Druck war, desto weiter mußte das Drahtmaterial 
aufgeheizt werden, bis die starke Expansion einsetzte 
(s. Tabelle 1, Spalte 5 und 11). Bei den Drähten mit 
0,24 mm Durchmesser, bei denen die magnetischen 
Drucke am größten waren, dürfte der Übergang zur 
starken Expansion ziemlich kontinuierlich erfolgt 
sein und zwar bei der kritischen Temperatur von 
9300°K. 


. III. Die einzelnen Phasen der Verdampfung 
von Drähten durch intensive Stromstöße 


1. Aufheizung bis zum Schmelzpunkt und Schmelzen 


Es soll nun zusammengefaßt werden, was sich auf 
Grund der experimentellen Ergebnisse über die Vor- 
gänge bei der Verdampfung von Drähten durch inten- 
sive Stromstöße sagen läßt. 

Der Beginn des Stromimpulses ist durch das Zün- 
den der Funkenstrecke gegeben, die für diese Betrach- 
tungen als idealer Schalter angesehen werden kann. 
Der Stromanstieg ist im ersten Moment 


dI D;(to) 

drin +Z 
Er betrug bei einer Gesamtinduktivität L,-+L von 
etwa 0,65u.H je nach den Aufladespannungen T, (t,) 
4 1010 bis 8,5 - 101% A/sec. Die Kondensatorspannung 
verteilt sich zu Beginn, da die ohmsche Spannung 


(20) 


r=I, 


noch vernachlässigbar ist, entsprechend den Induktivi- 


8500—9000| 17 


1,10 4,4 

0,85 | 4,5 9000— 9300| 20 

0,80 | 4,5 9000— 9300| 20 

0,68 | 4,65 9300 26 
8 ) 10 11 12 


täten auf den Entladekreis und am Versuchsdrah 
liegt die Spannung 


L 
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Da der Strom zunächst mit nahezu konstanter Ge 
schwindigkeit ansteigt, zeigt die an den Versuchs 
drähten gemessene Spannung zu Beginn einen horn 
zontalen Verlauf (s. Abb. 2). 
Der Draht wird durch die aufgenommene Joule=' 
sche Energie sehr rasch aufgeheizt. Die Stromvertei 
lung ist praktisch homogen (vgl. Anmerkung S. 52) | 
Dementsprechend wird auch die Temperaturvertei 
lung homogen sein, da die Wärmeverluste zunächst 
vernachlässigbar sind. Die Drahttemperatur T{ 
kann bis zum Schmelzpunkt gut nach der Beziehung 


1 
1 
| 
1 
| 


| 
| 
| 
j 
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berechnet werden, wobei j die Stromdichte in A/em 
ist. Für Kupfer haben die Stoffkonstanten folgen 
Werte: 

Spezifischer Widerstand o=1,7.10%Q.cm, 
sein Temperaturkoeffizient &, = 4,3 : 10°? grad’. 
Mittlere spezifische Wärme 
bis zum Schmelzpunkt 
Spezifisches Gewicht 


c = 0,1041 cal/g - grad. 
g = 8,93 g/em?. 


Da der Strom am Anfang der Entladung etwa mit 
konstanter Geschwindigkeit ansteigt, und der Schmelz 
punkt schon einige Zeit vor dem Strommaximum er 
reicht wird (vgl. Tabelle l, Spalte 8), kann man ir 
erster Näherung 


t er 
[ira = dr e=l 2: 
A — 


1 Diese Beziehung erhält man aus der kalorischen Zu 
standsgleichung Z(tl)=m + c{T(t)— T(t,)}, wobei E(t) 
vom Draht bis zur Zeitt aufgenommene Energie ist. 


2 
Dotto) -; (28 
(L&+Z) 


5] 

setzen, wobei d der Drahtdurchmesser in em ist. Die 
Kondensatoraufladespannung U,(t,) ist in Volt ein- 
zusetzen, die Induktivitäten L; und L in Henry. Aus 
61. (22) erhält man dann mit der Stoffkonstanten für 
Kupfer die Gleichung 

| 10 | _e) N%. 

OT) = 233 10 (rn) = 1.293. (4) 


‘Für den Zeitpunkt tg, zu dem im Draht die Schmelz- 
temperatur T,=1356° K erreicht wird, erhält man 
in dieser Näherung die Beziehung 

3 

| Delle) \3 

| en] 

‘die einen etwas zu kleinen Wert ergibt, da dj/dr kon- 
'stant angenommen wurde, während es in Wirklichkeit 
"etwas abnimmt. 

Der Schmelzpunkt wurde jeweils in ?/, bis ®/ı der 
‚Zeit erreicht, in der der Strom auf seinen maximalen 
"Wert anstieg. Der Drahtwiderstand war dann etwa 
'6- R, (R, = Drahtwiderstand bei 20°C) und stieg 
während des Schmelzens, das nur einige 103 sec in 
"Anspruch nahm, sehr rasch auf 12 R, an. Ein Knick 
im Verlauf der am Draht gemessenen Spannung und 
der Stromänderung (s. Abb. 2) ist diesem raschen An- 
"steigen des Drahtwiderstandes beim Schmelzen zuzu- 
ordnen. Der Schmelzvorgang war einige Zeit vor Er- 
reichen des Strommaximums beendet (vgl. Tabelle 1, 
Spalte 8). 


(25) 


2. Aufheizung des flüssigen Drahtmaterials bei nahezu 
konstantem Volumen, überhitzter Flüssigkeitszustand 


Nach dem Schmelzen wird das flüssige Metall 
immer schneller aufgeheizt. Im Bruchteil einer 
Zehntel-Mikrosekunde — ebenfalls noch vor dem 
 Strommaximum — wird in den Drähten die Tempe- 
ratur von 2570° K erreicht, bei der das Kupfer unter 
"Normaldruck zu verdampfen beginnt (s. Tabelle 1, 
‘Spalte 8). Infolge der Trägheitskräfte und des hohen 
magnetischen Druckes (in der Größenordnung von 
- 10: Atm) wird der flüssige Kupferzylinder ohne wesent- 
liche Volumenvergrößerung weiter aufgeheizt und das 
‘Kupfer befindet sich im Zustand einer überhitzten 
"Flüssigkeit. Wie in Abschnitt II.4 abgeschätzt 
_ wurde, setzt erst bei Temperaturen zwischen 8000 und 
9300°K von der Drahtoberfläche her eine starke 
- Expansion des Kupfers ein. Vorher dürfte der Kupfer- 
zylinder einigermaßen homogen aufgeheizt werden, 
da die Stromverteilung homogen ist und Strahlungs- 
_ verluste in diesem Temperaturbereich wegen der kur- 
zen Aufheizzeiten noch keine wesentliche Rolle spielen. 
“Man könnte deshalb bis zum Beginn der starken 
Expansion den zeitlichen Temperaturverlauf aus der 
“von den Drähten aufgenommenen Energie berechnen, 
wenn die spezifische Wärme des flüssigen Kupfers in 
_ diesem Temperaturbereich bekannt wäre. Leider ist 
sie nur bis 1500° C bekannt. Oberhalb 1500° C kann 
“ man die zu einem bestimmten Zeitpunkt im Draht 
erreichte Temperatur abschätzen, wenn man zwei 
_ extrem voneinander abweichende Annahmen über die 
spezifische Wärme des überhitzten flüssigen Kupfers 
macht. Einmal soll angenommen werden, daß die 
} spezifische Wärme oberhalb 1500° © nicht von der 
Temperatur abhängt und konstant den für 1500° 6 
geltenden Wert (Cs =0,12 cal/g- grad) hat. Man 
erhält dann Maximaltemperaturen (in Tabelle 1 mit 


| 
| 
I: 
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Tr, bezeichnet). Als zweites soll die Annahme gemacht 
werden, daß die spezifische Wärme über die bis 
1500°C bekannten Werte hinaus weiter linear an- 
steigt. Dies ergibt Temperaturen (in Tabelle 1 mit 7,,, 
bezeichnet), die etwa die untere Grenze der tatsächlich 
erreichten Temperaturen andeuten. Der Unterschied 
zwischen T,,, und 7, beträgt bei 4000°K etwa 10%, 
bei 10000°K etwa 30%. 

Im Rahmen der Genauigkeit dieser Abschätzung 
kann man sagen, daß die Temperaturen von 8000 bis 
9300° K, bei denen die explosionsartige Expansion 
beginnt, 2 bis 15 - 10°® sec nach dem Strommaximum 
erreicht werden (Tabelle 1, Spalte 2 und 9). 


3. Von außen nach innen fortschreitende 
starke Expansion, weiteres Aufheizen des Drahtinneren, 
Stromunterbrechung 


3.1. Modellvorstellung 

Die Expansion setzt von der Drahtoberfläche her ein 
und pflanzt sich mit Schallgeschwindigkeit zur Achse 
hin fort. In dem von der Verdünnungswelle erreichten 
Gebiet nehmen Dichte und Leitfähigkeit rasch ab. Der 
Strom verteilt sich deshalb nicht mehr homogen über 
den Drahtquerschnitt, sondern zieht sich mit dem Fort- 
schreiten der Verdünnungswelle immer mehr auf die 
Drahtachse hin zusammen. Das Drahtinnere, das von 
der Verdünnungswelle noch nicht erreicht wurde, wird 
bei konstantem Volumen weit über die kritische Tempe- 
ratur hinaus aufgeheizt, während außerhalb dieses 
Gebietes infolge der adiabatischen Expansion eine 
Temperaturerniedrigung eintritt. Es bildet sich ein 
starkes Temperaturgefälle. Hat die Verdünnungs- 
welle die Drahtachse erreicht, so beginnen auch dort 
Dichte und metallische Leitfähigkeit abzunehmen. 
Sobald die Dichte im Drahtinneren soweit gesunken 
ist, daß die Leitungselektronen nur noch auf Grund 
des Tunneleffektes von einem Atom zum anderen 
übergehen können, wird die metallische Leitung rasch 
unterbrochen (nach Davıp [17] ist dies bei einer 
Dichte des Kupfers von 0,5 bis 1,0 g/cm? der Fall, 
während die Dichte bei Zimmertemperatur 8,92 g/cm? 
beträgt). 

3.2. Expansion und Stromunterbrechung 

3.2.1. Bxperimentell bestimmte Abschaltzeiten 

Die Zeit zwischen dem Beginn der starken Expan- 
sion (Zeitpunkt t,) und der Unterbrechung der metal- 
lischen Leitung (Zeitpunkt t,) soll als Abschaltzeit 
bezeichnet werden. Da t, aus den experimentellen 
Ergebnissen abgeschätzt werden kann (Werte s. Ta- 
belle 1, Spalte 9) und t, aus den Stromkurven zu ent- 
nehmen ist (Werte s. Tabelle 1, Spalte 13), können die 
Abschaltzeiten angegeben werden. Sie nehmen mit 
dem Drahtdurchmesser merklich zu und liegen zwi- 
schen 0,1 und 0,5 usee (s. Tabelle 2, Spalte 6). 

3.2.2. Berechnung der Schallgeschwindigkeit im über- 
hitzten, flüssigen Kupfer 

Um zu prüfen, ob diese Abschaltzeiten mit der 
eben beschriebenen Modellvorstellung, nämlich der mit 
Schallgeschwindigkeit von der Drahtoberfläche nach 
innen laufenden Verdünnungswelle und der in ihrem 
Rücken einsetzenden Abnahme der Dichte, die schließ- 
lich zur Unterbrechung der metallischen Leitung 
führt, vereinbar sind, soll zunächst die Schallgeschwin- 
digkeit im überhitzten flüssigen Kupfer mit Hilfe der 
thermischen Zustandgleichung (18) berechnet werden. 
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Die Schallgeschwindigkeit v in Flüssigkeiten ist 
allgemein 
v2 (0P\ ‚% 
a 
wobei Z, der isotherme Volumelastizitätsmodul ist, 


c, und c, die spezifische Wärme bei konstantem Druck 
bzw. konstantem Volumen, g die Dichte und M das 


(26) 


I 


oP 
Molgewicht. Berechnet man c,, c, und (57 nr mit Hilfe 
der Gl. (18) und setzt die erhaltenen Ausdrücke in 
Gl. (26) ein, so erhält man für die Schallgeschwindig- 
keit die Beziehung 
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Abb. 7. ee) im Kupfer in Abhängigkeit 


von Temperatur 7 und Molvolumen V 


wobei Ac, die Abweichung der spezifischen Wärme im 
realen Zustand von der im idealen Gaszustand ist und 
den Wert 


2a 
AG = Tanz * 5 
Pe ne SITE ae n( (28) 
x| UTB Score Zi 


hat. Dabei sind B=V/-b, O0=T7T.:.V-c und D= 
V —d Abkürzungen und a, b,c, d, R sind die gleichen 
Konstanten wie in Gl. (18). Abb.7 zeigt die nach 
Gl. (27) berechneten Werte für die Schallgeschwindig- 
keit im überhitzten flüssigen Kupfer in Abhängigkeit 
von Temperatur 7 und Molvolumen V. 

3.2.3. Abschätzung der Laufzeit der Verdünnungswelle 

Eine Verdünnungswelle pflanzt sich bekanntlich 
mit der Schallgeschwindigkeit des ungestörten Me- 
diums fort. Für die Abschätzung der. Laufzeit, die die 
Verdünnungswelle von der Drahtoberfläche bis zur 
Achse benötigt, soll die vereinfachende Annahme ge- 
macht werden, daß sich die Verdünnungswelle mit 
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einer konstanten Schallgeschwindigkeit fortpflanzt, 
und zwar mit der, die dem Zustand zu Beginn der 
starken Expansion (Zeitpunkt t,) entspricht. Man er- 
hält unter dieser Annahme etwas zu lange Laufzeiten, 
da das Drahtinnere während des Fortschreitens der 
Verdünnungswelle bei konstantem Volumen weiter 
aufgeheizt wird und damit die Schallgeschwindigkeit 
wächst. Für den Zeitpunktt, können, wie gezeigt 
wurde, Temperatur und Molvolumen des Kupfers ab- 
geschätzt werden. In Tabellel, Spaltell und 12 
sind diese Werte für die verschiedenen Versuche zu- 


Tabelle 2. Schallgeschwindigkeit, Laufzeit der Verdünnungs- 
welle, Abschaltzeit und Dichteabnahme in der Drahtachse | 


Versuchs- 8 & 3 88 
parameter its) 5 © BE 
= .S= ER |: [25 
? | 8; = 2: | 38 | 8 8% 
Erle 58 Fi © S |88,5 
F S cH=| a er Ar ; 
A = >. | e el 
alıyı [nr | 1-| | u un 
mm mm) kV cem/sec usec ; Von usec" usec 9 
0,10 | 70| 28 3,0 | 0,017 135 0,023 10.12 
100) 40 2,6 0,019 1,4 h 
125| 50 2,6 | 0,019 1,4 
70) 50 1,8 0,028 1555 
= i | a 
0,195/100 25 2,6 0.037 1,4 
1125| 31,5] 2,6. | 0,097 | 1,4 
0,24 | 70| 28 1,8 0,067 | 1,55 
100 40 1105) 0,080 1,65 
125, 50 1,5 0,080 1,65 
70! 50 1,4 0,086 1,9 
Spalte 1 2 3 4 
sammengestellt. Für diese Werte erhält man aus \ 


Abb. 7 die in Tabelle 2, Spalte 2 eingetragenen Ge- 
schwindigkeiten in der Größenordnung von einigen 
10° cem/sec. Damit lassen sich sofort die Zeiten t, 
angeben, in denen die Verdünnungswellen Strecken 
durchlaufen, die den ursprünglichen Drahtradien ent- 
sprechen (Tabelle 2, Spalte 3). Um die Laufzeiten der 
Verdünnungswellen t, (Tabelle 2, Spalte 5) zu erhalten, 
müssen dr Zeiten t,, noch mit einem Faktor f (Ta 
belle 2, Spalte 4) multipliziert werden, der die vorn 
menausdehnung bis zum Beginn der starken Expan- 
sion berücksichtigt. Es zeigt sich, daß die Laufzeiten 
der Verdünnungswellen nur einen Bruchteil (!/, bis 1) h 
der Abschaltzeiten (Tabelle 2, Spalte 6) betragen. In E 
der restlichen Abschaltzeit müßte nach der Modell 
vorstellung die Dichte im Drahtzentrum soweit ab- 
nehmen, daß die metallische Leitung im Draht unter- © 
brochen wird. H 
3.2.4. Dichteverlauf in der Drahtachse 
Als letztes soll deshalb der räumliche und zeitliche 
Dichteverlauf im Kupferzylinder nach Beginn der 
starken Expansion betrachtet und abgeschätzt werden, | 
wann die Dichte im Drahtzentrum so stark gesunkenil 
ist, daß keine metallische Leitung mehr möglich ist. 
Abb. 8a zeigt schematisch das Weg-Zeit-Diagramm 
einer Verdünnungswelle, aus dem der räumliche 
Dichteverlauf für jeden beliebigen Zeitpunkt entnom- 
men werden kann. In Abb. Sb ist er für einen Zeit 
punkt gezeichnet, zu dem die Verdünnungswelle die 
| 


i 


ahtachse noch nicht erreicht hat. Qualitativ spielt 
a folgender Vorgang ab: Der Kopf der Verdünnungs- 
le bewegt sich mit der Schallgeschwindigkeit des 
yestörten Mediums, das sich im Gebiet I befindet 
„dius r, in Abb. 8b), in welchem die Dichte g,, die 
1allgeschwindigkeit vo, und die Teilchengeschwin- 
keit w,—0 sind. Hinter dem Kopf der Verdün- 
ngswelle befindet sich ein Gebiet II, in dem Schall- 
‚chwindigkeit vo und Teilchengeschwindigkeit u eine 
nktion von Ort und Zeit sind. Die Flächen kon- 
nter Zustandsgrößen laufen jedoch auch in diesem 
biet mit konstanter Geschwindigkeit — das ergibt 
'raden im Weg-Zeit-Diagramm — und zwar mit der 
schwindigkeit u+v. Die Drahtmasse wird in 
sem Gebiet kontinuierlich von uy=0 auf die Ge- 
© windigkeit der expandierenden Mediengrenze u —u, 
Üschleunigt und die Dichte nimmt stetig ab. Zur 
Sreinfachung soll angenommen werden, daß die 
S:diengrenze HI—IV plötzlich vom Zustand der 
he auf die konstante Expansionsgeschwindigkeit 4, 

schleunigt wird (sog. zentrierte Entlastungswelle), 
'& dies in Abb. 8a dargestellt ist. Das Ende der 
‘rdünnungswelle (Radiusr, in Abb. 8b) trennt 
ssen Bereich von einem Gebiet III, in dem wieder 
‘eim Gebiet I ein stationärer Zustand herrscht, und 
ar mit der Dichte g,, der Schallgeschwindigkeit v, 
d der Teilchengeschwindigkeit u,. Die Metall- 
mpffront, die in Abb.8b den Radius r, erreicht hat, 
Ihiebt eine Luftschicht vor sich her. In ihr bildet 
“h eine Stoßfront aus (Radius r, in Abb. Sb), die das 
“biet IV, in dem die Luft auf die Dichte 9, kompri- 
iert ist, vom Gebiet V trennt, in dem die Luft die 
'®sprüngliche Dichte g. hat. 

Für die Abschätzung des Dichteverlaufes im Ge- 
| let II, das sich, sobald der Kopf der Verdünnungs- 
jelle die Drahtachse erreicht hat, über das ganze 
h \ahtinnere erstreckt, soll ein ideales Gas angenom- 
Jen und nur der instationäre eindimensionale Fall, 
sr von CourANT und FriepricHs [1] behandelt 
irde, berechnet werden. Aus der Kontinuitätsglei- 
(nung, der Bewegungsgleichung, der Adiabatenglei- 
“hung! und der Zustandsgleichung für ideale Gase 
ıtgibt sich für die Fortpflanzungsgeschwindigkeit der 
lerdünnungswellen die Beziehung 


vwtv=vH4t an Ad a, (29) 


|lenn für das ideale Gas y =c,[c, —5/3 gesetzt wird. 
"as Ende der Verdünnungswelle pflanzt sich mit der 
‚eschwindigkeit 


1 E: 
%t%=% = = g% (30) 


"hrt. Es hängt also bei gegebenem », von der Expan- 
onsgeschwindigkeit u, ab, ob und wann das Ende der 


'erdünnungswelle die Drahtachse erreicht. Ist 


-dg, (31) 


P) 
I“, Bere 


b wird die linke Seite der Gl. (30) negativ, d.h. das 
"inde der Verdünnungswelle kann die Drahtachse 
\icht erreichen und die Dichte nimmt ständig ab. Von 
fürzer [15] wurde bei Versuchen, die den hier be- 
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"chriebenen sehr ähnlich waren, eine Expansions- 
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geschwindigkeit |u,| =3 km/sec gemessen. Bei dieser 
Expansionsgeschwindigkeit wäre für alle Werte, die 
sich für v, ergaben (1,4 km/sec<v,= 3,0 km/sec), die 
Ungleichung (31) erfüllt, d.h. die Dichte würde auch 
ohne Berücksichtigung einer Reflexion ständig ab- 
nehmen. 

Für den Dichteverlauf im Gebiet der Verdünnungs- 
welle Il erhält man — wenn man vom Einfluß der 
Reflexion in der Drahtachse absieht — aus den vorher 
erwähnten Gleichungen (s. Fußnote 1, S. 57) die Be- 
ziehung 


> 
2 ER 

g=m|t+ ) 2 z 

Da die vereinfachende Annahme einer „zentrierten‘“ 
Verdünnungswelle gemacht wurde, ist u+v —=gjt. 


(32) 
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Abb. 8a u.b. a Weg-Zeit-Diagramm einer „zentrierten‘“ Verdünnungswelle. 
» Räumlicher Dichteverlauf 


Damit erhält man aus Gl. (29) für die Teilchenge- 
schwindigkeit u als Funktion von Ort und Zeit 


af 


Setzt man diesen Ausdruck in Gl. (32) ein, so ergibt 
sich mit y=3 für den räumlichen und zeitlichen 
Dichteverlauf ohne Berücksichtigung einer Reflexion 


die Beziehung 
er ® 
Erz 


Für einige Versuche wurde der zeitliche Dichteverlauf 
in der Drahtachse nach Gl. (34) berechnet und in 
Abb. 9 dargestellt. Dabei wurde für x der Draht- 
radius zu Beginn der starken Expansion r=f ro 
eingesetzt. Die Werte für den linearen Ausdehnungs- 
koeffizienten f und die Schallgeschwindigkeit v, sind 
in Tabelle 2 enthalten. Aus Abb. 9 ist zu ent- 
nehmen, daß nach Zeiten, die den experimentell 
gefundenen Abschaltzeiten entsprechen, die Dichte 
in der Drahtachse etwa auf die Hälfte der Dichte zu 
Beginn der starken Expansion gesunken ist (8. Ta- 
belle 2, Spalte 7). Zusammen mit der Dichteabnahme 
vor der starken Expansion ergeben sich damit für den 


2 m € 
2 | = (33) 


(34) 
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Zeitpunkt der Stromunterbrechung ti, im Draht- 
zentrum Dichten von etwa 2 g/em? (Dichte des Kupfers 
bei Zimmertemperatur 8,92 g/em®). Die Dichte, bei 
der die metallische Leitfähigkeit verschwindet, liegt 
nach Davı [17] für Kupfer bei 0,5 bis 1,0 g/em?. 
Berücksichtigt man, daß die durchgeführte Ab- 
schätzung des Dichteverlaufes bestimmt eine zu lang- 
same Dichteabnahme ergibt, da der Einfluß der Re- 
flexion in der Drahtachse unberücksichtigt bliebund die 
Schallgeschwindigkeit v, nach Beginn der starken Ex- 
pansion infolge der weiteren Aufheizung des Draht- 
inneren noch ansteigt, so ist dieses Ergebnis in guter 
Übereinstimmung mit der Modellvorstellung, nach der 
bei den vorliegenden Versuchen die Unterbrechung der 
metallischen Leitung eine Folge der kontinuierlich von 
außen nach innen fortschreitenden Verdünnung des 
Metalls war. Es ist nicht anzunehmen, daß die Unter- 
brechung der metallischen Leitung durch Einschnürung 
der Flüssigkeitssäule und Zerfallen der Säule in einzelne 
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Zeitschri g 
angewandte P 


Geschwindigkeit 
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r, — Drahtradius zu Beginn der starken Expans) 
1,—t, = Abschaltzeit 


die man aus den experimentellen Ergebnissen bereef" 
nen kann, auf die Drahtachse zusammenschrump: 
Dann ist 


t—1,\? 
ty— be 


der Teil der gesamten Drahtmasse m,, der zur Ze 
noch weiter aufgeheizt wird. Teilt man die Absch 
zeit in Intervalle auf, so kann man für jedes die 
Intervalle die Temperatur berechnen, auf die 
Masse m(t) durch die zugeführte Energie aufgehe 
wird. Durch Extrapolation bis zum Zeitpunkt & 
Stromunterbrechung erhält man die Temperatur, di 
in der Drahtachse erreicht wird. Sie betrug na@ 
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Abb. 9 Abb. 10 
Abb.9. Zeitlicher Dichteverlauf in der Drahtachse unter verschiedenen Versuchsbedingungen 
Abb. 10. Zeitlicher Verlauf der Lichtemission (obere Kurve) und des Entladestromes (untere Kurve) während der Explosion eines Kupferdrahtes? j 


0,12 mm Durchmesser und 100 mm Länge. Kondensatoraufladespannung 25 kV 


Teile (Unduloide), die nicht mehr durch Funken über- 
brückt werden konnten, bewirkt wurde. Denn es 
läßt sich abschätzen, daß Ungleichmäßigkeiten in der 
Oberfläche und kleine Einkerbungen in den gemesse- 
nen Abschaltzeiten von einigen 10”sec durch die 
magnetischen Kräfte nicht bis zur Durchschnürung 
des Drahtes vertieft werden können. Ebenso läßt 
sich sagen, daß die stärkere Verdampfung an etwa 
vorhandenen dünneren Stellen in diesen kurzen Zeiten 
nicht zur Stromunterbrechung führen kann. 

3.3. Temperatur und Druck im Drahtinneren 

3.3.1. Abschätzung der im Drahtzentrum_ erreichten 
Temperaturen aus der aufgenommenen Energie 

Wie am Anfang dieses Kapitels erwähnt wurde, 
zieht sich der Strom nach Beginn der starken Expan- 
sion infolge der von außen nach innen fortschreitenden 
Dichteabnahme immer mehr auf die Drahtachse zu- 
sammen. Während dieser sog. Abschaltzeit wird den 
Drähten noch ein großer Energiebetrag zugeführt 
(s. Tabelle 1, Spalte 10 und 14) und das Drahtinnere 
damit weit über die kritische Temperatur hinaus auf- 
geheizt. Da die Stromverteilung in dieser Phase nicht 
genau bekannt ist, ist es nicht möglich, aus der auf- 
genommenen Energie die im Drahtzentrum erreichten 
Temperaturen zu berechnen. Man kann sie jedoch 
unter der Annahme abschätzen, daß das Gebiet, in 
dem.die Dichte noch groß genug ist, um gute metalli- 
sche Leitung zu gewährleisten, mit der konstanten 


Durchmesser und 100 mm Länge etwa 140000°] 
wenn man in erster Nährung mit der spezifischer. 
Wärme des flüssigen Kupfers bei 1500° C rechnet. 

3.3.2. Abschätzung der Strahlungsverluste 

Es soll nun noch kurz die Frage der Energieverlust® 
behandelt werden. Die Verluste durch Wärmeleitur 
und Konvektion sind bei der nur etwa 10° sec dauert 
den Aufheizung vernachlässigbar klein. Inwiewe 
Energie durch Strahlung verloren geht, wurde bereit 
in einer früheren Arbeit [20] in vereinfachten Betrad 
tungen untersucht. Unter der Annahme einer unen 
lich großen Wärmeleitfähigkeit im Kupferzylinder er 
gab sich, daß bei Zufuhr einer konstanten Leistung 
von 108 bis 10° W — diese Leistung wurde den Drähte 
bei den hier. beschriebenen Versuchen im zeitliche 
Mittel zugeführt — die Strahlungsverluste im Bereie 
von 10? und 105 Grad merklich werden, Berücksicl 
tigte man die endliche Wärmeleitfähigkeit, so erga 
sich für die bei den Versuchen erzielten Leistung 
dichten im stationären Fall ein parabelförmig von de 
Drahtachse nach außen abfallender Temperaturver 
lauf. Obwohl der stationäre Zustand in der zur Vei 
fügung stehenden Zeit keineswegs erreicht wird, wir 
sich dennoch nach Beginn der starken Expansio) 
ein Temperaturgefälle ausbilden, da sich der Stron 
immer mehr auf die Drahtachse zusammenzieht un 
sich die Randgebiete infolge der adiabatischen Ex 
pansion stark abkühlen. Einen deutlichen Hinwe: 


\ 
£ dieses Auftreten eines starken Temperaturgefälles 
\fert der experimentell gewonnene zeitliche Verlauf 
'r Lichtemission (s. Abb. 10). Es zeigt sich nämlich 
stens, daß das Maximum der Lichtemission vor dem 
rommaximum liegt, die Liehtemission also bereits 
ı einem Zeitpunkt abzunehmen beginnt, zu dem das 
rahtmaterial sicherlich noch weiter aufgeheizt wird, 
d zweitens, daß der Lichtimpuls schneller abklingt, 
\s es bei einer über den Drahtquerschnitt konstanten 
emperatur der Fall wäre. Da in den Abschaltzeiten 
on einigen 10°” sec nur relativ wenig Energie aus 
sm Drahtinneren an die Oberfläche transportiert 
ird, wird das Drahtinnere durch die kühleren äußeren 
chichten gegen Strahlungsverluste abgeschirmt und 
4 können in der Achse die berechneten hohen Tem- 
eraturen entstehen. 
13.3.3. Abschätzung der im Drahtzentrum erreichten 
Irucke 
‘ Mit Hilfe der allgemeinen thermischen Zustands- 
leichung (19) läßt sich auch der Druck abschätzen, 
er im Zentrum der Drähte herrschte. Das Mol- 
olumen des Kupferdrahtes mit 0,24 mm Durchmesser 
d 100 mm Länge dürfte vor der Stromunter- 
jrechung zwischen 20 cm?/Mol — dem Wert zu Be- 
sinn der starken Expansion — und 60 em?/Mol — das 
»ntspricht der Dichte von 1 g/em?, bei der etwa die 
inetallische Leitfähigkeit verschwindet — betragen 
1aben. Das ergibt mit der Temperatur von 140 000° K 
sinen Druck zwischen 2,3 und 11,4 - 10° Atm. 


Zusammenfassung 


Kupferdrähte von 0,035 bis 0,24 mm Durchmesser 
ınd 70 bis 125mm Länge wurden durch Entladen 
einer Kondensatorbatterie (Kapazität 1,3 uF, Auflade- 
pannung 25 bis 55 kV) stoßartig verdampft. Aus 
genauen Messungen des zeitlichen Verlaufes des Ent- 
ladestromes und der am Versuchsdraht liegenden 
Spannung wurde der zeitliche Verlauf des Draht- 
widerstandes und der von den Drähten aufgenomme- 
nen Energie berechnet. Aus diesen Ergebnissen erhält 
man durch Elimination der Zeit die Energieabhängig- 
keit des Drahtwiderstandes, die bis zur Siedetempera- 
tur des Kupfers bei Atmosphärendruck (etwa 2570 KR) 
gut mit den aus Tabellen entnommenen Werten über- 
einstimmt. Das Drahtmaterial wurde jedoch durch 
den nach innen gerichteten magnetischen Druck, der 
durch das Feld der Entladeströme (sie betrugen im 
Strommaximum bis zu 35 kA) hervorgerufen wurde 
und in der Größenordnung von 10* Atm lag, weit 
“über diese Temperatur hinaus zusammengehalten. 
“Unter Benutzung einer allgemeinen thermischen Zu- 
"standsgleichung von HIMPAN, die sich bereits in ande- 
‘ren Fällen bei hohen Drucken gut bewährt hat, er- 
gibt sich aus den experimentellen Daten, daß sich 
der flüssige Kupferzylinder je nach den Versuchs- 
bedingungen bis zu Temperaturen zwischen 8000 und 
9300° K relativ wenig ausgedehnt hat. Erst bei diesen 
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Temperaturen setzte eine starke Expansion des Draht- 
materials ein. Eine Verdünnungswelle pflanzte sich 
mit Schallgeschwindigkeit von der Drahtoberfläche 
nach innen fort, in deren Rücken Dichte und Leit- 
fähigkeit abzunehmen begannen. Während ein noch 
leitender Kupferzylinder, der sich immer mehr auf die 
Drahtachse zusammenzog, weiter aufgeheizt wurde, 
kühlten sich die Randgebiete durch Strahlung und 
infolge der adiabatischen Expansion ab. Das ent- 
stehende Temperaturgefälle schirmte die inneren Ge- 
biete stark gegen Strahlungsverluste ab und ermög- 
lichte es, daß das Drahtmaterial dort durch die bis 
zur Stromunterbrechung aufgenommene Energie auf 
hohe Temperaturen gebracht wurde. Eine Abschät- 
zung ergibt, daß in der Drahtachse Temperaturen bis 
etwa 150000°K und Drucke von einigen 10° Atm 
erreicht wurden. 


Herzlichen Dank schulde ich Herrn Professor 
Dr. MArer-Leisniz, an dessen Institut ich diese Ar- 
beit anfertigen durfte. Herrn Professor Dr. FÜNFER 
möchte ich vielmals für seine Unterstützung und für 
wertvolle Hinweise, Herrn Dipl.-Phys. LEHNER für 
anregende Diskussionen danken. 
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Temperatur und Zündfähigkeit von Schleiffunken in verschiedenen Gasgemischen 


Von Hans WAHL 
Mit 6 Textabbildungen 
(Eingegangen am 12. Oktober 1959) 


RIEZLER und HArDr haben eine Methode beschrie- 
ben, die es erlaubt, Messungen an Schleiffunken durch- 
zuführen [1]. Diese Methode wurde noch etwas ver- 
bessert und dazu benutzt, um die Abhängigkeit der 
Funkentemperatur von der Art und Zusammen- 
setzung der den Funken umgebenden Atmosphäre zu 
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Abb.1. Schematische Anordnung des optischen und elektrischen Teiles der 
Apparatur 


messen. Dabei wurden insbesondere explosive Gas- 
gemische untersucht. 

Die Funken wurden erzeugt durch Anpressen 
eines Stahlstiftes gegen eine mit 2800 Umdrehungen 
pro Minute rotierende Schleifscheibe aus Edelkorund, 
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Abb. 2. Mittelwerte der Funkentemperaturen bei Zufügung von Sauerstoff 
oder Stickstoff zu Luft in Abhängigkeit von der Konzentration des der Luft 
zugemischten Gases 


Körnung 36, Härtegrad M, keramische Bindung, Durch- 
messer 115mm. Sie befand sich in einer vakuum- 
dichten Kammer, in die man Gasgemische einstell- 
barer Konzentration einfüllen konnte und die auch bei 
Zündung explosiver Gasgemische nicht zerstört wurde. 
Die Funken flogen an 1 cm breiten Fenstern in ver- 
schiedenen Abständen von der Schleifscheibe vorbei, 
von wo ihr Licht durch einen Lichtleiter auf die Photo- 
kathoden zweier Sekundärelektronenvervielfacher ge- 
langte. Der optische Teil der Apparatur wurde gegen- 
über dem in [1] benutzten aus Gründen der höheren 
Lichtintensität und besseren Geometrie etwas ab- 
geändert. Er ist zusammen mit dem elektrischen Teil 
in Abb. I schematisch dargestellt. Die Änderung be- 
stand darin, daß das aus dem Lichtleiter Z austretende 
Funkenlicht durch einen halbdurchlässigen Spiegel S 
in zwei senkrecht zueinander stehende Bündel auf- 
geteilt wurde. Durch die Linsen L, und Z, wurde das 
hintere Ende des Lichtleiters verkleinert auf die Photo- 
kathoden der SEV S, und S, abgebildet. Die Linse Z, 
diente als Kondensor, um möglichst alles den Lichtleiter 
verlassende Licht auf die Photokathoden zu bringen. 


Zeitschrift für 
angewandte Phys 


Vor beiden Photokathoden befanden sich Metall-Inteı 
ferenzfilter M, und M,, deren Schwerpunktswellen 
längen A,=1769,9mu. und /, —=601,5 mu betrugen. Die 
beiden Vervielfacher waren an zwei Breitbandverstärke 
V, bzw. V, angeschlossen, deren Ausgänge zu je einem 
System eines Zweistrahloszillographen Z führten, des= 
sen Meßplatten entgegengesetzt gepolt waren. Die Os“ 
zillogramme wurden mit einer Laufbildkamera photo 
graphiert. Aus dem Verhältnis der sich entsprechen 
den Impulshöhen beider Oszillogramme kann mai 
nach Eichung mit einer Wolfram-Bandlampe auf die 
Temperatur der Funken schließen. 

Die Messungen zeigten, daß der Gang der Abhän 
gigkeit der Funkentemperatur von der Art und 
Konzentration der den Funken umgebenden Gase bei 
verschiedenen Stahlsorten ähnlich ist. Daher sind im 
folgenden aus zahlreichen Meßreihen nur die Meß 
ergebnisse als Beispiel aufgeführt, die mit einer be: | 
stimmten Stahlsorte gemacht wurden, nämlich einem F 
gehärteten unlegierten Kohlenstoffstahl mit 1% Koh- 
lenstoff, 0,18% Silizium und 0,19% Mangan. | 

Alle hier aufgeführten Messungen wurden in einem 
Abstand von 8 cm von der Schleifscheibe gemacht, da” 
sich dort in Luft das Maximum der Temperaturmitt 
werte der Funken befand. Bei geringerer Sauerstoff 
konzentration verschob sich das Maximum zu größeren 
Abständen hin und umgekehrt, da das Leuchten des” 
Funkens von einer Verbrennung des abgerissenen Teil-” 
chens herrührt. | 

Abb. 2 zeigt die Mittelwerte der Funkentempera= 
turen bei Zufügung von Sauerstoff oder Stickstoff zu 
Luft. Die Abszisse zeigt dabei die Konzentration des 
der Luft zugefügten Sauerstoffs oder Stickstoffs in 
Volumprozenten an. Wie nicht anders zu erwarten, 
nimmt die Funkentemperatur mit wachsendem Saue 
stoffgehalt zu, während sie mit zunehmendem Stie 
stoffgehalt, also abnehmendem Sauerstoffgehalt, stark 
abnimmt. In einer reinen Stickstoffatmosphäre sind 
überhaupt keine Funken sichtbar, ein Beweis dafür, 
daß das abgerissene Stahlteilchen erst durch die be 
seiner Verbrennung entstehende Energie eine so hohe”) 
Temperatur erreicht, daß es für das Auge sichtbar 
leuchtet. 

Weiter wurde der Einfluß folgender, mit Luft ge- 
mischter brennbarer Gase und Dämpfe untersucht: 
Methan, Wasserstoff, Leuchtgas (Bonner Stadtgas), ® 
Äthylen, Kohlenmonoxyd, Äthyläther, Methylalkohol, 
Aceton und Benzol. Abb. 3 enthält die Meßergebnisse 
bei brennbaren Gasen, Abb.4 die bei brennbaren 
Dämpfen. In den Bereichen, in denen das Gemisch 
zündete, konnte natürlich nicht gemessen werden; 
diese Bereiche sind in den Abbildungen punktiert. 

Auch hier fällt erwartungsgemäß die Funken 
temperatur mit abnehmendem Sauerstoffgehalt, wenn 
auch bei den verschiedenen Gasen unterschiedlich. 
Die Abnahme der Funkentemperatur ist bei den brenn- 
baren Gasen und Dämpfen mit abnehmender Sauer 
stoffkonzentration viel stärker als bei Stickstoff oder 
Argon. Abb. 5 zeigt die Meßergebnisse bei einigen der 
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tersuchten Gase und Dämpfe, gemischt mit Luft, Schleiffunken zur Zündung ausreichen müßte, kann 
zogen auf gleiche Sauerstoffkonzentration. Tat- der Grund für diesen Effekt nur der sein, daß die Fun- 
hlich weichen die Kurven für Stickstoff und Argon ken so schnell durch ein sonst zündbares Gemisch 
n dem schraffierten Gebiet, in dem die Kurven fliegen, daß sie an einer Stelle nicht genügend Energie 
mtlicher untersuchter brennbarer Gase und Dämpfe abgeben. Große Funken müßten, obwohl sie im all- 
‚gen, erheblich ab. Die zunächst naheliegende Ver- gemeinen geringere Temperatur haben, demnach eher 
utung, daß die Wärmeleitfähigkeit der Gase eine zur Zündung führen. Es wurde tatsächlich festgestellt, 
itscheidende Rolle spiele, trifft nicht zu, da die daß nicht jeder Funken zündet, denn auch bei bri- 
ärmeleitzahlen der untersuchten Gase mit Aus- santem Gemisch dauerte es bei starker Funkenbildung 
hme von Wasserstoff alle in derselben Größenord- bis zu 10 sec, ehe das Gemisch gezündet wurde. 
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Abb. 3 Abb. 4 Abb. 5 

turen bei Zufügung brennbarer Gase (Abb. 3) und brennbarer Dämpfe (Abb. 4) zu Luft in Abhä 
von der Konzentration des der Luft zugemischten Gases 

en bei Zufügung einiger Gase und Dämpfe zu Luft in Abhängig 


stoffs im Gemisch 
ung liegen. Die Erklärung ist offenbar folgende: An Um nun auch den Einfluß der verschiedenen Wär- 
er heißen Oberfläche der Funken wirken die brenn- meleitfähigkeiten der Gase zu untersuchen, wurden 


aren Gase reduzierend, d.h. sie entziehen durch ihre B - 
igene Oxydation der unmittelbaren Umgebung des 2200 -z 
Fa 


bb. 3 u. 4. Mittelwerte der Funkentempera ngigkeit 
bb. 5. Mittelwerte der Funkentemperatur keit von der Konzentration des Sauer- 


kens einen Teil des Sauerstoffs, so daß das be- 6 A 
reffende Eisenteilchen viel langsamer und daher bei I | 
ieferer Temperatur verbrennt. Dabei braucht die 200 WR 
)xydation des brennbaren Gases nicht unbedingt HE TE 
chnell genug zu sein, um das Gasgemisch zu zünden, 1800 22 £ 
umal der Funke ja sehr schnell durch das Gasgemisch £ [ 


liegt. Fi 


Tabelle 1 gibt den Bereich in Volumprozenten an, j 
_ dem das vorliegende Gas oder der Dampf bei ee 
ischung mit Luft durch Schleiffunken gezündet a —* 2 
wurde: Bie) ES set 
| 
Tabelle 1 Be eo el 
2 Gas bzw. Dampf | Zündbar zwischen | j | 
He. "ER 5—15% =, ex 5 
Abb. 6. Mittelwerte der Funkentemperaturen bei Zufügung von Stickstoff, 
Wasserstoff . .. 0 - 7,5—18% Argon uder sa Aunknsieneit ler he iesibe 
10— 35% des Sauerstoffs im Gemisch 
EA 6,5—20% 
Kohlenmonoxyd . ..... . > 15—30% Stickstoff, Argon und Helium direkt mit Sauerstoff 
ee ro: © hichk zündbar! statt mit Luft gemischt. Diese Gase sind nicht brenn- 
Methylalkohol . ... + - ab 7% = e B 
ton. rel 3_15% bar und verändern deshalb nicht durch eigene Oxy- 
2,5—8% dation die Funkentemperatur, außerdem haben sie 


sehr verschiedene Wärmeleitzahlen, wie aus Tabelle 2 
zu ersehen ist: 


Auffallend ist, daß der Bereich, in denen die Ge- 


mische durch Schleiffunken gezündet wurden, im Tabelle 2 

allsemeinen verhältnismäßig eng ist. So wurde z.B. 
en uk lether.-Lnfe- Gemisch überhaupt nicht ge- en ser 
zündet, obwohl dessen Zündtemperatur bei einer ande- | 

ren Art der Zündung bei 180° C liegt. Auch ein Was- Sauerstoff . - 0,587 - 107 
serstoff-Luft-Gemisch wird normalerweise bei einer Be a en i IR 
Wasserstoffkonzentration von 4 bis 75% bei einer a 3.430 - 10 


Zündtemperatur von 510° gezündet, während im 
Falle der Schleiffunkenzündung der Zündbereich nur 
zwischen 7,5 und 18% liegt. Da die Temperatur der 


Die Ergebnisse der Messungen zeigt Abb.6. Man 
sieht, daß sich bei höherer Konzentration der inaktiven 


Ss bye 1 A 
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Gase ein Einfluß der Wärmeleitfähigkeit bemerkbar 
macht, und zwar ist, wie vorauszusehen, die Funken- 
temperatur bei gleicher Sauerstoffkonzentration um so 
geringer, je größer die Wärmeleitfähigkeit des Gas- 
gemisches ist, weil die Wärme schneller abgeleitet wird. 


Zusammenfassung 


Es wurde die Temperatur von Schleiffunken in 
verschiedenen Gas- und Dampfatmosphären unter- 
sucht, indem das Verhältnis der Intensitäten des 
Funkenlichtes bei zwei verschiedenen Wellenlängen 
bestimmt wurde. Dabei ergab sich folgendes: 

Die Schleiffunken heizen sich während ihres Fluges 
durch Verbrennung auf; in Atmosphären ohne Sauer- 
stoff gibt es keine sichtbaren Funken. Die Maxima 
der Temperaturmittelwerte liegen bei niedrigen Tem- 
peraturen weiter weg von der Schleifscheibe als bei 
hohen. 

Mischungen von Luft mit brennbaren Gasen in 
geringen Konzentrationen zeigen eine viel stärkere 
Abnahme der Funkentemperatur mit abnehmender 
Sauerstoffkonzentration als Mischungen von Luft mit 
nicht brennbaren Gasen. Die Wärmeleitfähigkeit der 


Über den Einfluß von Querschnittsänderungen und Gittern in Kanälen auf Stoßwellen 


Von Heınz ReicHEnBAcH und HELMUT DREIZLER 
Mit 13 Textabbildungen 
(Eingegangen am 7. September 1959) 


1. Einleitung 


Der Einfluß von Gittern auf Stoßwellen wurde 
bisher von 1.1. GLass, W. MArTtm und G.N. PATTER- 
son [1], [2], von I.E. BowmAn und G.B.F. Nıs- 
LETT [3] und von D.S. Dosanst# [4], [5] behandelt. 
Im folgenden sollen diese Veröffentlichungen ergänzt 
und erweitert werden. 

Die theoretische Untersuchung über den Einfluß 
von Querschnittsänderungen und Gittern auf schwa- 
che Stoßwellen wird auf eine Näherungslösung füh- 
ren, die nicht auf einer Lösung der nichtlinearen Glei- 
chungen der Gasdynamik bei entsprechenden Rand- 
bedingungen beruht, sondern direkt von den Erhal- 
tungssätzen für Masse, Impuls und Energie ausgeht. 

Das Hindurchlaufen einer Stoßwelle durch das 
Gitter denken wir uns aus zwei nacheinander ab- 
laufenden Vorgängen zusammengesetzt, nämlich dem 
Eintritt der Stoßwelle in einen verengten Kanal und 
dem Austritt aus diesem verengten in einen erweiter- 
ten Kanal. Die Wirkungen beider Teilvorgänge sum- 
mieren sich zu der Gesamteinwirkung eines Gitters 
endlicher Dicke auf die Stoßwelle. Da sich aber viele 
Hindernisse als eine Folge von Querschnittsänderun- 
gen auffassen lassen, sind die oben genannten Teil- 
vorgänge über die in dieser Arbeit beabsichtigte 
spezielle Anwendung hinaus wichtig. 

Durch zweckmäßig angelegte Experimente prüfen 
wir die theoretischen Überlegungen auf ihre Trag- 
weite und versuchen, Abweichungen der experimen- 
tellen von den theoretischen Ergebnissen zu deuten. 


1 Die Ziffern beziehen sich auf das Literaturverzeichnis. 
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zu Luft zugemischten Gase hat bei diesen geringer 
Konzentrationen keinen merkbaren Einfluß auf di 
Funkentemperatur. Erst bei hohen Konzentrationen 
ist eine Abhängigkeit der Funkentemperatur von d 
Wärmeleitfähigkeit des verwandten Gasgemisches 
erkennen, und zwar liegen bei hohen Wärmeleitzahlı 
der untersuchten Atmosphären die Funkentempera 
turen niedriger und umgekehrt. 'W 

Ob ein Gemisch von Luft mit brennbarem Ga 
oder Dampf durch Schleiffunken gezündet wird, häng? 
nicht so sehr von der Funkentemperatur, sondern vom 
Konzentrationsverhältnis ab. Es wurde festgestellt 
daß die Grenzen der Konzentrationsverhältnissef 
innerhalb deren das Gasgemisch gezündet wurde, 
Schleiffunkenzündung teilweise wesentlich enger siı 
als bei einer anderen Art der Zündung. Äthyläthe 
Luft-Gemische werden durch Schleiffunken überhaup 
nicht gezündet. 


Literatur: [1] RıezLer, W., u. L. Harpr: Z. angew. Phys 

6, 497 (1954). — [2] KOHLRAUSCH, F.: Praktische Phys 
Bd.2. Leipzig 1951. 

Dipl.-Phys. Hans WAHL, 

Institut für Strahlen- und Kernphysik 

der Universität Bonn 


2. Einige Eigenschaften der Stoßwelle 


Bevor auf das eigentliche Problem der Arbei 
eingegangen wird, sollen einige Eigenschaften de 
Stoßwelle zusammengestellt werden?. Aus einer Ver 
dichtungswelle endlicher Amplitude entsteht in einen 
Gas eine Stoßwelle, da sich die Front der Verdich 
tungswelle aufsteilt. Die Stoßwelle kann man als’ 
einen Sprung des Druckes 9, der Dichte 0, der Tem 
peratur 7’ und der Strömungsgeschwindigkeit u be 
schreiben. Nur in der Zone des Sprunges, die in deı 
Größenordnung der freien Weglänge der Molekü 
liegt, spielen innere Reibung und Wärmeleitung eine 
wesentliche Rolle. Sie haben zur Folge, daß die 
spezifische Entropie eines Massenelementes über das 
eine Stoßwelle hinweggeht, ansteigt. Dennoch ist eine 
einfache mathematische Behandlung der Stoßwelle 
möglich, welche die irreversiblen Prozesse in dem Stoß 
nicht berücksichtigen muß, da sie auf den Erhaltungs: 
sätzen von Masse, Impuls und Energie beruht, welche 
von dem irreversiblen Vorgang in der Stoßfront nich 
beeinflußt werden. Wir setzen Erhaltungsgleichungen 
in einem Koordinatensystem an, das sich mit dem 
Stoß $ bewegt, der die Geschwindigkeit w, relativ 
zum ruhenden Koordinatensystem habe (Abb. 1). (Die 
Größe $ bezeichnet keinen Vektor. Sie ist das g 
bräuchliche Zeichen für eine Stoßwelle. Der Pfeil 
gibt die Richtung der Ausbreitung an.) 

Wenn man die Massen-, Impuls- und Energie- 
stromdichten jeweils vor und hinter dem Stoß einander 


® Umfassende Einführungen geben COURANT-FRIRDRICHS 
[6] und Guass, MArrtın, PATTERsoN [1]. 


\ 


02 %2 = 01 Vı> (1) 

+) =M+t (0, %ı) dı; (2) 

Up Pa + (02 95) („B+3W} (3) 
= u, Pı + (Q1 v1) (,n+2u)} 


ei 
%=UrT Ws = L52 


Zusammen mit den Zustandsgleichungen — wir 
‚en im Rahmen dieser Arbeit ideale Gase voraus — 


P=0, RT, (4) 


damit vier Gleichungen gegeben. Wenn die 
"Ben im Gebiet 1 vor dem Stoß bekannt sind, so 
bleiben im Gebiet 2 fünf unbekannte Größen Ps, 
T,, u, und w,. Wählt man eine der Größen als 
ameter — es ist üblich ?2, zu nehmen —, so sind 
‚übrigen vier Größen durch vier Gleichungen be- 
mt. Häufig wählt man noch für den Parameter p, 
dimensionslose Form P,, =P3/Pı- 


3. Theorie über den Einfluß von Hindernissen 
in Kanälen auf Stoßwellen 


3.1. Übersicht über den Weg der Untersuchung 
j _ und die Voraussetzungen 


Wir wollen uns im folgenden Abschnitt bemühen, 
. Veränderungen theoretisch zu erfassen, die eine 
Bwelle beim Auftreffen auf einfache Hindernisse 
eidet, ohne auf den Einfluß von Grenzschichten 
zugehen. Diese Hindernisse sind Querschnitts- 
derungen, und zwar: 

1. Die Verengung des Querschnitts eines Kanals 
bb. 2), 

%, die Erweiterung des Querschnitts eines Kanals 
bb. 3), 

3. das Gitter, das den freien Querschnitt eines 
ınals vermindert (Abb. 4). 

Die Hindernisse werden in der folgenden Theorie 
ein durch das Verhältnis des Querschnitts Q,; den 
der Strömung belassen (vgl. Abb. 2, 3, 4), zum 
erschnitt Q, des freien Kanals charakterisiert. Die 
ezielle Form des Übergangs bleibt unberücksichtigt. 
Im Rahmen dieser theoretischen Betrachtungen 
eine stetige Veränderung des Querschnitts einer 
stetigen, sprunghaften Veränderung des Quer- 
hnitts gleichwertig. Die Verengung und Erweite- 
g des Querschnitts eines Kanals wurden gesondert 
ıfgeführt, obwohl man geneigt ist, sie in einer einzi- 
n Betrachtung zusammenzufassen und sie durch 
n Querschnittsverhältnis Q,/@, zu beschreiben, das 
ıtweder <1 oder >1 sein kann. Eine solche ein- 
itliche Beschreibung haben wir jedoch nicht ge- 
ählt, weil bei einer Querschnittsverengung des 
anals eine Stoßwelle, bei einer Erweiterung hingegen 
ne Verdünnungswelle reflektiert wird. Über die 
oßwelle vergrößert sich die spezifische Entropie, 
hrend sie über die Verdünnungswelle konstant 
bt [6]. Dieser wesentliche Unterschied läßt sich 
lein dureh die Angabe des Querschnittsverhältnisses 
‚[@: (größer oder kleiner als eins) nicht erfassen. 


} ‚1 Wegen der Symbole s. 8. 71. 
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Deshalb werden in den nachstehenden Untersuchun- 
gen Verengung und Erweiterung gesondert behandelt. 
Natürlich müssen beide Betrachtungen in dem Grenz- 
fall des freien Kanals Q,/Q, =1, zum gleichen Ergebnis 
führen. 

Das Gitter soll als eine Aufeinanderfolge einer 
Verengung und einer Erweiterung behandelt werden. 
Die Beschreibung des Einflusses eines Gitters auf eine 
Stoßwelle wird also eine Verknüpfung der Ergebnisse 
der beiden erstgenannten Fälle sein. Obwohl es sich 
bei den Vorgängen um zweidimensionale Probleme 
handelt, führt die gegebene „quasi-eindimensionale“ 
Theorie erstaunlich weit. Sie ist immer dann möglich, 
wenn die Größen, die den Strömungsvorgang be- 
stimmen, über den Querschnitt als annähernd kon- 
stant betrachtet werden können. Dies darf man in 


2 
ee. ee 


er 
Abb.1. »-x-Diagramm einer Stoßwelle S mit nachfolgender Kontakt- 


<> 
front C in einem abgegrenzten System S,, $, zum Zeitpunkt £ 


Gebieten, die dem Hindernis nicht unmittelbar be- 
nachbart sind, und zwar um so eher, je geringer die 
Querschnittsänderung durch das Hindernis ist. In 
der Nähe der Hindernisse muß diese vereinfachte 
Theorie versagen, denn dort sind die Vorgänge sicher 
nur in mindestens zwei räumlichen Dimensionen zu 
beschreiben. Wir beschränken uns auf die Frage nach 
der Veränderung von Stoßwellen durch starre Hinder- 
nisse. Wie Stoßwellen auf die Hindernisse selbst 
wirken, wie also der Druck am Hindernis verteilt ist, 
ist ein schwieriges Problem zu dessen Lösung die Vor- 
gänge in der Nähe der Hindernisse ergründet werden 
müssen. Es wird in dieser Arbeit nicht behandelt. 
(Einen Ansatz zur Berechnung der Gesamtkraft auf 
ein Hindernis findet man in [7].) 


3.2. Die Verengung des Querschnitts eines Kanals 


Trifft eine Stoßwelle auf eine Verengung, so kann 
man sich nach experimentellen Erfahrungen und 
theoretischen Überlegungen folgendes vereinfachtes 
Strömungsschema (Abb. 2) aufzeichnen. Es treten 
eine durchgehende Stoßwelle 8 und eine reflektierte 
Stoßwelle S auf, die in gewisser Entfernung vom 
Hindernis wieder ebene Fronten aufgebaut haben. 
Die Kontaktfront € trennt das Gas, das sich zur Zeit 
des Auftreffens der Stoßwelle auf das Hindernis im 
weiteren Kanal mit dem Querschnitt Q, befand, von 
dem Gas, das den engeren Querschnitt ausfüllte. 
Druck und Nachströmgeschwindigkeit sind vor und 
hinter der Kontaktfront © gleich, dagegen sind Dichte, 
Temperatur und Entropie verschieden, da das Gas 
im Gebiet? nur der Wirkung der durchgehenden 
Stoßwelle S, während das Gas im Gebiet 3 sowohl 
der Wirkung der einfallenden Stoßwelle S, als auch 


der reflektierten Stoßwelle N unterliegt. Die auf- 
tretenden „‚transversalen‘“ Stoßwellen [1] haben einen 
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geringen Einfluß und werden vernachlässigt. Sie 
entstehen aus den anfänglichen zylinderartig ge- 
krümmten, von den scharfen Kanten des Hindernisses 
ausgehenden Stoßwellen beim allmählichen Übergang 
zu der später ebenen Front der Stoßwelle. 

In unsere Theorie gehen noch die folgenden zwei 
Voraussetzungen ein. Wir nehmen an, daß 

1. der Fluß zwischen dem Gebiet 5 und 3 adiaba- 
tisch und stationär und 

2. die Mach-Zahl im Gebiet 3, <1 sei. 


a) 
b) 
d 
ar Yen 
3 N 
RI 
AUS 
? 30 ER 
Ä ERS 
N 
LS 
PERS 
= 70 T 
bi a0 N 
2, ARE 
SE 
20 SS 
N 
RS 
5 S 
715 S ze Ä u 
een WS 
N RS 
ne N TE OR 
a ee — -_ 
70 18 20 De. 


Abb. 23—c. Verhalten einer Stoßwelle in einem Kanal, dessen Querschnitt 
sich verengt. a Schematische Darstellung des Kanals. b Vereinfachtes 


— — 
Strömungsschema. Mit S,, Sund Ssind die Stoßwellen, mit Gdie Kontakt- 
front und mit Q,, Q. die Querschnitte bezeichnet. ce Abhängigkeit der Stoß- 


stärken der an einer Verengung reflektierten (B) und durchgehenden (R ‚) 

Stoßwelle von der Stoßstärke P, der einfallenden Welle für verschiedene 

Querschnittsverhältnisse ö6. Für die gestrichelten Kurven gilt die obere 
Skala von P;o 


Die letztgenannte Voraussetzung läßt sich auf- 
heben. Da aber die Versuche dieser Arbeit die Vor- 
aussetzung 2 stets erfüllen, wird die hier mögliche 
Erweiterung der Theorie nicht verfolgt. Sie bleibt 
einer späteren Veröffentlichung vorbehalten. 

Das vorliegende Problem ist gelöst, wenn zu den 
vorgegebenen Werten 


Pı> 01» T: u) 
Po» Q6; Te; Us, Ws 


die folgenden 14 Unbekannten bestimmt sind: 


a a 6 > Sa =, 
P5> 05: T,, Pa» 03» R, Pa, 02, R, is 
> = >> > en 

Us, Wz Ug Un, Wa. 


* Das Zeichen > soll andeuten, daß es sich um einen 
Strömungsvorgang in einem Kanal mit enger werdendem 
Qustschnitt handelt. 


Zunächst gilt, wie in Abschnitt 2 festges 
Stoßwelle 8 


<) 


wurde, über die durchgehende 
selbstverständlich auch über die reflektierte 


welle 8 hinweg jeweils der Erhaltungssatz der M 
(1), des Impulses (2) und der Energie (3). 


Kontaktfront 6 ist: 


= > 
Pa — P3» 
> > 

Ug —= Ug. 


Da wir vorausgesetzt haben, daß die Luft wie@ 
ideales Gas zu behandeln sei, gilt jeweils im Gebie 
3 und 5 die Zustandsgleichung 

RT Eh 

Mit (1), (2), (3), (5), (6) und (7) sind 11 der n@ 
wendigen 14 Bestimmungsgleichungen vorhanden. 
fehlenden drei Gleichungen erhält man durch Ve 
knüpfung der Größen in den Gebieten 3 und 5 (Abb 
An dieser Stelle ist auch das Maß Q,/Q, für die Que 
schnittsverengung einzuführen. 


Dazu verwenden wir die Kontinuitätsgleichung 
Massenstroms 


Q 2, 2 Qı u; AR 


die Bernoullische Gleichung für stationäre Strömun 


> > 

3 il Ge Y P3 1 22 

h a 
„1. 2 Ya 2 


und die Adiabatengleichung, welche Druck- ı 
Dichteänderung eines Gases bei gleichbleiben« 
spezifischer Entropie verknüpft, 


und erhalten als bisher fehlende drei Verknüpfun 


relationen 


(@r- a (ren ( 
a m \e+o-nM 
Be „(ten er 
5 \2+w-DM | 
> a | 
za - (Ze) Br | 
“ En 1 
% \2+W-)2® | 


Die Gl. (11) wird in der quasi-eindimensionalen Theoms 
der Düsenströmung häufig verwendet [8]. Beder 
man, daß die Adiabatengleichung (10) eine zeit 
konstante Entropie eines Massenelementes vora 
setzt, diese Voraussetzung bei Strömungen ideal 
reibungsfreier Gase sich aber wiederum aus der Ft 
derung von der Erhaltung der Masse, des Imp 
und der Energie ableiten läßt [1], so stellen die dr 
Gln. (11), (12) und (13) letzlich nur eine Konseque 
der Erhaltung von Masse, Impuls und Energie vi 
bunden mit der Voraussetzung der Adiabasie u 
Stationarität des idealisierten Vorganges dar. 

Da eine allgemeine explizite Auflösung der 14 Gl 
chungen (1), (2), (3), (5), (6), (7), (11), (12) und (1 
nicht möglich ist, verbleibt nur eine numerische B 
rechnung. Sie wird vereinfacht, wenn man sch 
berechnete und tabellierte Werte der Beziehunge 
über die Stoßfronten verwendet [1], [7]. 
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_ Glücklicherweise eröffnet die Art der 14 nicht- 
nearen Bestimmungsgleichungen F, (& ... %4) =, 
—_1]...14 einen einfachen Lösungsweg. In der 


abelle1l ist das Ergebnis mitgeteilt. Die Größen 
- > 
31 und P;,, die nach dem in [7] angegebenen Rech- 


ungswege zuerst erhalten werden, waren einer di- 
»kten Messung mit vorhandenen Geräten nicht zu- 


änglich, so daß in Tabelle 1 die mittels EL und C, 


> >> > == 
erechneten Werte Iz,, |W;ı| und 15), WM, aufge- 
ommen wurden, die wir mit unseren Messungen 


=; 
päter vergleichen wollen. Die Stoßstärke R,, die 


'abelle 1. Ergebnisse der numerischen Rechnung für die Vor- 
inge beim Auftreffen einer Stoßwelle auf eine Querschmitts- 


verengung 
741919191 nn — — — 
sı ö DB ER er Wal Te W., 2: 
| - | 
200 \1,4665| — 1310097 | — | — |. — 
0,1 1,390 | 0,460 | 1,262 | 0,928 1,539 | 1,314 1,850 
0,3 |1,267 | 0,408 | 1,182 | 0,877 1,472 | 1,275 | 1,733 
0,5 1,173 | 0,369 | 1,119 | 0,839 | 1,420 | 1,245 | 1,645 
0,7 |1,093 | 0,336 | 1,064 | 0,805 | 1,378 | 1,222 | 1,578 
0,9 1,030 | 0,310 | 1,020 | 0,776 | 1,346 1,204 | 1,526 
1,0 1,000 | 0,298 | 1,000 | 0,763 | 1,332 | 1,195 | 1,500 
»0| 0,0 11,874 | — |1,552 | 0,94 — == er 
0,1 |1,748 | 0,754 | 1,480 | 0,898 | 1,227 | 1,551 | 2,645 
0,3 1,533 | 0,683 1,351 | 0,817 | 1,845 | 1,492 | 2,431 
| 0,5 \1,353 | 0,629 | 1,239 | 0,744 | 1,770 | 1,448 | 2,280 
| 0,7 |1,200 | 0,581 | 1,137 | 0,677 | 1,720 | 1,408 | 2,150 
0,9 1,064 | 0,540 | 1,044 | 0,616 1,645 | 1,375 | 2,043 
‚1,0. \1,000 | 0,523 | 1,000 | 0,587 1,623 | 1,362 | 2,000 
,5| 0,0 2,234 | — |1,746 | 0,946 | — =. En 
‚0,1 2,075 | 0,982 | 1,662 | 0,891 2,277 | 1,752 | 3,420 
0,3 1,783 | 0,900 | 1,501 | 0,783 | 2,157 | 1,678 | 3,122 
0,5 1,534 | 0,831 | 1,351 | 0,685 | 2,058 | 1,618 | 2,890 
0,7 |1,318 | 0,777 | 1,215 | 0,593 | 1,978 | 1,571 | 2,715 
0,9 |1,108 | 0,727 | 1,074 | 0,497 | 1,904 1,529 | 2,560 
1,0 |1,000 | 0,707 | 1,000 | 0,447 | 1,879 | 1,511 | 2,500 


ebenfalls in Tabelle 1 angegeben ist, dient weiteren 
Betrachtungen bei der Theorie des Gitters. In Abb. 2 
ist ein Teil der Ergebnisse der Tabelle 1 in der Form 
von Diagrammen dargestellt. 

Der Sonderfall Q,/Q,=d=1 entspricht, wie es 
sein muß, der ungestörten Stoßwelle, da der Kanal 
hindernisfrei ist. Ebenso liefert der Sonderfall des 
vollständig geschlossenen Kanals, ö6=Q,/0. =, die 
Werte für die vollständige Reflexion [1]. 


3.3. Die Erweiterung des Querschnitts eines Kanals 


Bevor die Verhältnisse beim Übergang einer Stoß- 
welle in eine Erweiterung untersucht werden, sollen 
zunächst ganz allgemein die Massen- und Energie- 
stromdichte in einer Stoßwelle bestimmt werden, da 
wir diese nachfolgend benötigen. Hierzu wird ein 
System SS, (Abb. 1) abgegrenzt, in das die Stoß- 
welle hineinläuft. 


Zur Zeit t befindet sich in dem System eine Masse 
M,—= Q {2305 + %0 + % 01} (14) 
(Q = Querschnitt des Kanals). 

Zur Zeit t+dt 


M;;a: = Qflz +u,dt) 05 + 
+ (2, +w,dt — u,dt) oa + 


| Z. f. angew. Physik. Bd. 12 


15 
(& — wzdt) Q1}- Sn 


Hieraus ergibt sich die Massenstromdichte durch die 
Systemgrenze S, 


1 dM 

Im Oo dt 

oder in dimensionsloser Schreibweise auf a, und o, 
bezogen 

Din = WM, ll <- (17) 


Um die Energiestromdichte an der Stelle 5, zu 
berechnen, setzt man den Energieinhalt im System 
an und erhält 


= w, (0 — 9) + U (0 — 05), (16) 


1) + Il): 


PR w. plc, Rn+ “ u) — 0% 7) + = (0 — 0) (18) 


oder in dimensionsloser Schreibweise auf a, und 9, 
bezogen 


> 


1) Y D) A, t | 
a 


ri] 


Ist die Kontaktfront C noch nicht im System $, 8, 
(Abb. 1), so fallen in den Gln. (16) bis (19) die Aus- 
drücke weg, in denen 0,— 0, bzw. I, —/', steht. 
Wie beim Auftreffen einer Stoßwelle auf eine Ver- 
engung wird auch beim Übergang einer Stoßwelle in 
einen Kanal weiteren Querschnitts ein vereinfachtes 
Strömungsschema (Abb. 3) angenommen. In der 


Erweiterung ändert die ankommende Stoßwelle S, 
ihre Stärke, die Stoßfront krümmt sich, um sich der 
veränderten Geometrie anzupassen, geht aber nach 
Durchlaufen eines Übergangsgebietes wieder in eine 
ebene Front S über. Dem Stoß $ folgt eine von der 


Verengung ausgehende Kontaktfront ©. Sie trennte 


(19) 


1 


beim Ankommen der einfallenden Stoßwelle & den 
Gasinhalt im engeren Kanal vom Gasinhalt des wei- 
teren Kanals. Die Gase im Gebiet 2 und 3 unter- 
scheiden sich in ihrer Dichte, Temperatur und Entro- 
pie, da sie beim Durchgang der Stoßwelle verschiede- 
nen Vorgängen unterworfen waren. Beim Übergang 
der Stoßwelle in den Kanal weiteren Querschnitts 
tritt im engeren Kanal eine reflektierte Verdünnungs- 
welle auf [1], [9], die hier nicht berücksichtigt wird. 
Die Berechtigung dazu entnehmen wir letzlich der 
Übereinstimmung der Theorie mit den ausgeführten 
Messungen (Abb. 8), bei denen schwache Stoßwellen 
bis zur Stoßstärke P,, =2,5 verwendet wurden. 

Eine Abzählung der Unbekannten des vorliegenden 
Problems ergibt folgende 9 Größen 

= < 
Das 05 D, Dan Der Tara, W2, 

Bekannt sind wiederum 


= 


SZ = NE 
9, 9,4 u = 0 


< 


< < 
Pe» 06; T,; Ug, Wg- 

Die Erhaltung der Masse, des Impulses und der 
Energie über die Stoßfront S hinweg, sowie die Be- 


ziehungen < 
: nd, (20) 
üg = Up (21) 


für die Größen beiderseits der Kontaktfront © und 
die Zustandsgleichungen in den Gebieten 2 und 3 


3 
ans 2) 
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ergeben sieben der notwendigen neun unabhängigen 
Bestimmungsgleichungen. Wir verwenden nun noch 
die Gleichheit des Massen- und Energiestroms im 
Gebiet 6 und 3 und erhalten damit die zwei noch 
fehlenden Bestimmungsgleichungen für die 9 Unbe- 
kannten. 


a) 
b) 
pP 
25 30 67 35 40 
BI T Pl v2 1 
< EL Sr La 
Br Y Ph 


[4 


Abb. 3a—e. Verhalten einer Stoßwelle in einem Kanal, dessen Querschnitt 
sich erweitert. a Schematische Darstellung des Kanals. b Vereinfachtes Strö- 


> oo — 
mungsschema. S,, 5 Stoßwellen; C. Kontaktfront. 
< 


Stoßstärke P,, einer in eine Erweiterung eintretenden Stoßwelle von der 
Stoßstärke 7%, der einfallenden Welle für. verschiedene Querschnitts- 
verhältnisse ö. Für die gestrichelten Kurven gilt die obere Skala von Po 


c Abhängigkeit der 


An dieser Stelle ist es wohl nötig, einige Gedanken 
über die Verwendung des Impulssatzes einzufügen. Im 
vorigen Abschnitt wurden zur Verknüpfung der 
Größen der dortigen Gebiete 5 und 3 die drei Gln. (11) 
bis (13) benutzt, die letzlich aus den Erhaltungssätzen 
von Masse, Impuls und Energie in differentieller Form 
folgen. Das System, in dem der Impuls erhalten 
bleibt, ist das Gas zusammen mit dem Hindernis im 
Kanal. Der dort gegebene Lösungsweg gestattete es, 
die drei GlIn. (11) bis (13) aufzustellen, ohne die 
Druckverteilung am Hindernis zu kennen. Bei den 
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Strömungen dieses Abschnittes bleibt der Impuls i 
dem System, bestehend aus dem Gas und dem Hi 
dernis im Kanal, selbstverständlich auch erhalten.) 
Aber die Erhaltungssätze von Masse und Energie 
welche die Größen in den Gebieten 6 und 3 verknüpfen 
haben eine integrale Form. Wollte man in einer ent 
sprechenden Form den Impulssatz hinzuziehen, 
müßte man hier die genaue Druckverteilung am Hin 
dernis kennen. Dieser Schwierigkeit ist man aber ent 
hoben, da die obengenannten 9 Beziehungen ausrei- 
chen, die 9 Unbekannten zu bestimmen. Der Impuls” 
satz bietet eine Möglichkeit, die in dieser Art gegebene” 
Theorie zu erweitern. | 
Nach den obigen Ausführungen lautet die Gleichung 
für die Erhaltung des Massenstroms | 


< < < E= = 
Q Wil I) Mila -1)+ U, la-Izı)} (2 
und die für die Erhaltung des ‚Energiestroms 


1 U 
Werl Pa—l+ Zu 


ER A 


1, i )+ı & 
0% 
AR Sr we By}. 

Löst man die GIn. (23) und (24) nach & auf und 
setzt sie einander gleich, so erhält man die folgende 
Bestimmungsgleichung (25), in der die Größen des 
Gebietes 6 bekannt sind und alle Größen des Ge 


nur von B, abhängen [1]. 


< 
zus D2 
LEE a At ae 3) 2 


Qı 
1 U2 4 
= Wat Ba—D+Y glg — (2 
< 
— J = Dame 
IB) 4yrnl. 


= 
Die Stoßstärke RB, ist darum die als erste zu bestim 
mende der 9 Größen. Die übrigen Unbekannten lassen 


sich dann mittels B, auf einfache Art berechnen [7] 
Das Ergebnis der Rechnung ist in dem Diagramm 
(Abb. 3) und in der Tabelle 2 angegeben. Die Werte 


< 
der aus BD, berechneten Stoßwellengeschwindigkeit W;, 
dienen zum Vergleich der dargelegten Theorie mit dem 
Experiment. 


3.4. Das Gitter 


Um den Vorgang des Aufpralls einer Stoßwelle auf 
ein gitterähnliches Hindernis und ihres Durchgangs 
theoretisch zu erfassen, konstruieren wir uns im fol 
genden Abschnitt, ähnlich wie in den vorangegange- 
nen, ein vereinfachtes Strömungsschema (Abb. 4), wo- 
zu wir experimentelle Erfahrungen und Angaben aus 
[1] benutzen. Es bildet sich, genau wie früher, eine 


durchgehende Stoßwelle S und eine reflektierte Stoß 


welle S aus. Auch eine Kontaktfront O tritt auf. Alle 
anderen Erscheinungen wollen wir vernachlässigen 
Sie werden später bei der Deutung der Abweichungen 
von Theorie und Experiment 'eine Rolle spielen. Wir 
verknüpfen jetzt die Ergebnisse aus den Abschnitten 
3.2 und 3.3, indem wir annehmen, die Stoßwelle ver 
halte sich beim Eintritt in das Gitter so, als ob sie 


ı DA 4 Ps 


XI. Band 
‚eft 2 — 1960 


E 
abelle 2. Ergebnis der Rechnung für den Übergang einer 
Stoßwelle in eine Erweiterung 


Pıı 1/6 B, 7A Pu 1/6 B.. Wi 
10 | 01 |100 I1000| 33 | 0,1 | 1,20 | 1,081 
03 | 1.00 | 1,000 03 | 151 | 1199 
0,5 | 1.00 | 1,000 05 | ım7 | 1288 
0.7 | 1.00 | 1,000 07 | 2.00 | 1.362 
0.9 | 1.00 | 1,000 10 | 230 | 1.454 

1,0 [1,00 |1,000| 95 | 0,1 | 1,23 | 1,093 

13 | 01 |108 | 1012 03 | 1.605 | 1,232 
03 | 1.09 | 1,037 05 | 190 | 1,330 
05 | 116 | 1,066 07 | 216 Lau 
07 |122 | 1.089 09 | 239 | 1,480 
10 \130 | 1.120 10 250 | 151 

15 01 vossIi1oael 30 | 01 | 1,34 | 1,136 
03 | 1,17 | 1,070 03 | 18a | 1311 
05 128 |L112 05 | 223 | 1,432 
07 | 137 | 1.147 07 | 257 | 1,530 
09 1146 | 1181 09 | 2865 | 1e1ı 

10 1.50 | 1.195 10 | 300 | 1,647 

18 oı [nur Inos| 35 | 01 | 15 | 1177 
03 | 1.295 | 1.118 03 | 2.08 | 1,387 

05 | 1.46 | 1181 05 | 2,56 | 1,528 

07 | 1,60 | 1,231 07, 298 | 1.642 

10 1.80 | 1,298 09 | 334 | 1,732 
20 0111 [1,087 IK MEN FL LLEE 
03 18 |vısıl 30 | 01 | 1565 | 1,218 
05 | 1,58 | 1,224 03 | 2325 | 1.461 

07 | 176 | 1,284 05 , 290 | 1.621 
09 | 1.92 | 1,336 07 | 3.405 | 1,749 

10 2.00 | 1,362 09 | 381 | 1845 

10 | 400 | 1,887 


sine Verengung des Querschnitts vorfinde. Die Er- 


gebnisse für die reflektierte Stoßwelle S bei der Quer- 
schnittsverengung (Abschnitt 3.2) werden daher auf 


die reflektierte Stoßwelle S beim Gitter übertragen. 


Tabelle 3. Ergebnis der Rechnung für den Durchgang einer 
Stoßwelle durch ein Gitter 


Rı ö Ps Te Wil T5ı W.ı 
—— lm — 
1,5 0,0 1,4665 | 1,310 0,957 | 1,000 1,000 
0,1 1,390 1,262 0,928 | 1,080 1,048 
0,3 1,267 1.218907 0,8772, 1,181 1,106 
0,5 1,173 1,119 0,839 | 1,248 1,147 
0,7 1,093 1,064 | 0,805 | 1,285 1,168 
0,9 1,030 | 1,020 | 0,776 | 1,320 1,189 
1,0 1,000 1,000 | 0,763 | 1,332 1,195 
2,0 0,0 1,874 1,552 0,944 | 1,000 1,000 
0,1 1,748 1,480 | 0,898 | 1,178 1,104 
0,3 1,533 1,351 0,817 | 1,373 1,220 
0,5 1,353 1,239 | 0,744 | 1,487 1,284 
0,7 1,200 1,137 0,677 | 1,555 1,320 
0,9 1,064 1,044 | 0,616 | 1,600 1,349 
1,0 1,000 1,000 | 0,587 | 1,623 1,360 
2,5 0,0 2,234 1,746 | 0,946 | 1,000 1,000 
0,1 2,075 1,662 | 0,891 1,287 1,170 
0,3 1,783 1,501 0,783 | 1,567 1,331 
i 0,5 1,534 1,351 0,685 | 1,707 1,411 
0,7 1,318 1,215 | 0,593. | 1,798 1,464 
0,9 1,108 1,074 | 0,497 | 1,854 1,497 
1,0 1,000 | 1,000 0,447 | 1,879 1,511 


Die Ergebnisse für die durchgehende Stoßwelle bei 


der Querschnittsverengung, speziell B,, werden als 
" Ausgangswerte für den Vorgang genommen, der sich 
abspielt, wenn die durchgehende Stoßwelle aus dem 
Gitter wieder in den freien Querschnitt übertritt. 
Letzteres entspricht genau dem behandelten Fall der 
Querschnittserweiterung. 
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Somit erhält man ohne weiteres das Diagramm 
Abb. 4 und die Tabelle 3, welche zusammen mit dem 
Diagramm Abb. 2 die Wirkung eines Gitters auf eine 
auftreffende Stoßwelle beschreiben. 

Bei der geschilderten Aneinanderheftung von 
Querschnittsverengung und -erweiterung zur Be- 
schreibung des Vorgangs am Gitter ist vorausgesetzt, 
daß der Raum innerhalb des Gitters, also dessen 
„Dicke“, ausreicht, um aus der beim Eintritt in das 
Gitter etwas deformierten Stoßfront noch vor dem 


WILDE, 


8) 
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Abb. 4a—c. Verhalten einer Stoßwelle in einem Kanal, in den ein Gitter 
eingebaut ist. a Schmatische Darstellung des Kanals mit Gitter. b Verein- 
fachtes Strömungsschema. € Abhängigkeit der Stoßstärke RB, einer Stoß- 
welle, die durch ein Gitter hindurchgegangen ist, von der Stoßstärke R, der 
einfallenden Welle für verschiedene Querschnittsverhältnisse 6 


Austritt in den freien Querschnitt wieder eine ebene 
Front entstehen zu lassen. Erfahrungsgemäß geschieht 
dies in dem verengten Querschnitt viel schneller als 
bei der reflektierten Stoßwelle im weiteren Kanal. 

Weiter wird angenommen, daß die beiden Kon- 
taktfronten, die sich eigentlich nach den kombinierten 
Strömungsschemata (Abb. 2 und 3) bilden müßten, 
zu einer einzigen Kontaktfront zusammenfallen. Diese 
Voraussetzung erscheint berechtigt, wenn man be- 
denkt, daß das Volumen des Gases zwischen den 
Gitterstäben klein ist verglichen mit dem Gasvolumen, 
das zwischen der Ebene des Gitters und dem Ort der 
gemessenen Stoßfronten liegt. 


4. Messungen 
4.1. Beschreibung der Stoßwellenapparatur 


Für die Erzeugung von Stoßwellen für experimen- 
telle Untersuchungen benutzten wir eine Stoßwellen- 
apparatur in üblicher Anordnung [1], [2], [3], [&), [7]- 


Sie besteht aus einem Stoßwellenrohr von 92 x 92mm 
5* 
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Querschnitt. Die Kompressionskammer des Rohres 
ist 0,65 m, die Expansionskammer ist 4,30 m lang. 
Die Entfernung zwischen der Membran und Meß- 
kammer beträgt 3,25 m. Für die Interferenzaufnah- 
men stand ein Mach-Zehnderscher Interferenzrefrak- 
tor mit einem Spiegeldurchmesser von 18 cm zur Ver- 
fügung. Das Stoßwellenrohr wurde dabei mit seiner 
Meßkammer senkrecht in den einen Teilstrahlengang 
des Refraktors gebracht. Als Lichtquelle für die Ein- 
zelaufnahmen diente eine Xenonfunkenstrecke. Zu 
Schattenaufnahmen verwendeten wir die übliche An- 
ordnung [4] mit einer punktförmigen Funkenlicht- 
quelle. 


4.2. Prinzipielle Meßanordnung im Stoßwellenrohr 


Die theoretischen Überlegungen des Abschnittes 3 
beschrieben, wie Stoßwellen durch Hindernisse der 
genannten Art verändert werden. Im folgenden be- 
sprechen wir unsere in großer Zahl ausgeführten 


Abb.5. Ausschnitt aus dem Stoßwellenrohr mit eingebauter Teilwand Zi 
und Gitter@. Zustand nach dem Auftreffen einer Stoßwelle 5; auf ein 
Gitter @. Es treten e@ ine rı reflektierte Stoßwelle S 


, eine durchgehende Stoß- 
welle Ss und eine Kontaktfront C auf 


Messungen, die dazu dienen sollten, die Güte und den 
Geltungsbereich unserer Theorie durch Vergleich mit 
dem Experiment zu prüfen. 

Man läßt eine Stoßwelle, die im 'Stoßwellenrohr 
mit einer gewünschten Stärke erzeugt werden kann [1] 
auf ein Hindernis auftreffen und untersucht die dabei 
entstehenden neuen Stoßwellen. Hierzu müssen die 
Messungen jeweils eine die Stoßwelle vollständig cha- 
rakterisierende Größe zu ermitteln gestatten. Diese 
sollte sowohl meßtechnisch leicht bestimmbar als auch 
mit dem Ergebnis der Theorie leicht vergleichbar sein. 
Häufig wird der Druck p, (oft in der Form ,,) als 
charakterisierende Größe der Stoßwelle verwendet. 
Hier benutzen wir stattdessen gleichwertig die Dichteo 
(in der Form T},) oder die Stoßw ellen gene 
(in der Form 17 7), weil mit der vorhandenen Appara- 
tur der Druck hinter den Fronten der neu entstehen- 
den Stoßwellen auf direktem Wege nicht gemessen 
werden kann. 

Eine Interferenzaufnahme mit dem Mach-Zehnder- 
Interferometer in einem geeigneten Zeitpunkt ergibt 
mittels der dabei gemessenen Änderung des Brechungs- 
index n der Luft eine Aussage über die Verteilung der 
Dichte o;, speziell über das Verhältnis 7}, der Dichten 
beiderseits der auftretenden Stoßfronten. 

Um eine Dichte- mit der Geschwindigkeitsmessung 
in einer Aufnahme zu kombinieren, wurde das Stoß- 
wellenrohr durch eine Teilwand in zwei gleiche Kanäle 
(Abb. 5) geteilt. In dem einen Teilkanal wurde das 
Hindernis angebracht, während der andere Teil- 
kanal frei blieb. Die Geschwindigkeiten der am Hin- 
dernis entstehenden Stoßwellen wurde durch Ver- 
gleich mit der Geschwindigkeit der im hindernisfreien 


(Abb. 8) findet man die Abhängigkeit der Geschwin 


Zeitschrift fü 
angewandte Phy: 


Teilkanal ungestört laufenden Stoßwelle bestimm 
Es verhalten sich nämlich die Geschwindigkeiten de 
Stoßwellen wie die von ihnen durchlaufenen Strecke 
Die Geschwindigkeit der im hindernisfreien Teilkan 
entlanglaufenden Stoßwelle ist, wie een un 
theoretisch bis zu Stoßstärken P,| =2,5 geprüft wurde 
innerhalb der Meßgenauigkeit gleich der aus den An 
fangsbedingungen bekannten Stoßwellengeschwindigf 
keit vor Auftreffen auf die angeschärfte Teilwand. 
Um die Störungen in den Gebieten hinter des 
Front der einfallenden Stoßwelle in den Teilkanäleı 
möglichst gering zu halten, wurde die Teilwand auf 
der Stirnseite messerscharf unter einem Winkel von 
co —=4° symmetrisch angeschrägt. Außerdem liest da 
Hindernis etwa 25 cm von der Stirnseite der Teilwanc R 
entfernt. Auf dieser Strecke können sich kleine Stö 
rungen noch ausgleichen. Durch einen Vergleich de e 
Lage der Stoßfronten in beiden Teilkanälen vor der 
Hindernis wurde geprüft, ob beide Kanäle glei hä 
wertig sind. 


$ 


4.3. Meßergebnisse 


In den nun folgenden Tabellen und Diagrammei 
sind die Meßergebnisse zusammengestellt. In den 
Diagrammen geben die ausgezogenen Kurven jeweik 
das Ergebnis der Theorie an. Die senkrechten Strecken 
an den Meßpunkten bezeichnen den Größtfehler. Bei 
den Abbildungen der verwendeten Modelle wird stets SS 
angenommen, daß die einfallende Stoßwelle von links 
nach rechts läuft. 


4.31. Verengung des Querschnitts eines Kanals 
Messungen der reflektierten Stoßwelle 


Das Modell bestand aus 3mm starken Platten) 
mit ebener Stirnfläche und scharfen Kanten. Bei 
einem Querschnittsverhältnis ö=0,315 wurden 2], 
bei ö=0,706 wurden 9 Platten auf die Breite des 
Kanals von 92mm verteilt. Die Meßwerte für die! 
beim Eintritt einer Stoßwelle in eine Verengung‘ 
reflektierte Stoßwelle sind nicht in ein gesondertes 
Diagramm eingetragen, sondern in die Diagramme! 
(Abb. 9) für die an Gittern gleicher Querschnitts- 
geometrie reflektierten Stoßwellen aufgenommen. 


Messungen der durchgehenden Stoßwelle- 


Als Modelle dienten hier einmal eine Stufe, die 
entsprechend Abb. 2in den Teilkanal eingebaut wurde. 
Zum anderen wurde die Querschnittsverengung nicht 
durch eine Stufe mit einer senkrechten ebenen Vorder- 
seite erreicht, sondern durch eine allmähliche Ver- 
engung, wie unten in Abb. 6 angedeutet ist. Die Meß- 
werte sind in dem Diagramm (Abb. 6) in Abhängigkeit 
vom Druckverhältnis A,, der einfallenden Stoßwelle 
dargestellt. Darin gibt die Kurve W,, die Geschwindig- 
keit der jeweiligen einfallenden Stoßwelle an. 

In dem Dingen (Abb. 7) ist die Stoßwellen- 


geschwindigkeiten W bei konstantem Druckver- 
hältnis R,, in Abhängigkeit von dem Querschnitts- 
verhältnis ö dargestellt. 


4.32. Erweiterung des Querschnitts des Kanals 


Das Modell hatte hier wieder die Form einer Stufe, | 
wie sie in Abb. 3 angegeben ist. In dem Diagramm 


digkeit W,, der durchgehenden Stoßwelle von dem 


ckverhältnis P,, der einfallenden Stoßwelle. Das 
uerschnittsverhältnis dient als Parameter der Kur- 


n. 
4.33. Gitter 


Die meisten Messungen wurden an Gittern vorge- 
mmen. Die Gitter bestanden aus Stäben von qua- 
atischem Querschnitt mit einer Seitenlänge von 
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= 
'b.6. Abhängigkeit der Geschwindigkeit W,, der durchgehenden Stoß- 
lle von der Stoßstärke A, der einfallenden Stoßwelle bei konstantem 
uerschnittsverhältnis 6=0,3. Die Geschwindigkeit W,, der ungestörten 
; Stoßwelle ist zum Vergleich eingetragen 


mm. Sie waren an beiden Enden gelagert. Die An- 
hl der Stäbe, die über den freien Querschnitt des 
anals gleichmäßig verteilt wurden, war jeweils so 


15 = 


2, 


} > 
bb. 7. Abhängigkeit der Geschwindigkeit W,, der durchgehenden Stoß- 
elle vom Querschnittsverhältnis ö bei konstanter Stoßstärke P; = 1,85. 
Die Geschwindigkeit W,, ist zum Vergleich eingetragen 


ewählt, daß sich das angegebene Querschnittsverhält- 
is ö ergab. Die Diagramme (Abb. 9 und 10) geben 
och die Messungen an dem eben erwähnten Platten- 


um 75 20 25 


% 


= 
bb. 8. Abhängigkeit der Geschwindigkeit W,, der durchgehenden Stoß- 
elle von der Stoßstärke 3, der einfallenden Stoßwelle bei konstantem 
Querschnittsverhältnis 1/6 = 0,7 bzw. 1/6=0,3. Die Geschwindigkeit Wsı 
der ungestörten Stoßwelle ist zum Vergleich eingetragen 


odell und zum Vergleich die Messungen von DOSANJH 
[4] an Lochplatten mit einer gleichen Querschnitts- 
erengung an. 

Ergänzend wurde noch ein Gittermodell unter- 
sucht, das in Abb. 11 wiedergegeben ist. Man kann 
es als Ausschnitt aus einem Gitter auffassen. Im Dia- 
gramm (Abb. 11) ist die Abhängigkeit der Geschwin- 
digkeit W,, der durchgehenden Stoßwelle von der Stoß- 
stärke P,, der einfallenden Stoßwelle aufgetragen. Das 
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Abb. 9. Abhängigkeit des Dichteverhältnisse I5s und der Geschwindig- 
keit W,, der an einem Gitter mit Stäben quadratischen Querschnitts bzw. 
an einer Verengung reflektierten Stoßwelle vom Querschnittsverhältnis ö 
bei verschiedenen Stoßstärken A,. I Gitter aus Stäben mit quadrati- 
schem Querschnitt; A Lochplatte nach [4]; x Plattenmodell 
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Abb. 10. Abhängigkeit des Dichteverhältnisses I}, und der Geschwindig- 

keit W,, der von einem Gitter aus Stäben quadratischen Querschnitts 

durchgelassenen Stoßwelle vom Querschnittsverhältnis ö bei verschiedenen 

Stoßstärken Pı. 5 Gitter aus Stäben quadratischen Querschnitts; 
A Lochplatte nach [4] 


Querschnittsverhältnis 6 = 0,3 ist konstant. Für die 
Stoßstärke P,, =1,93 wurden die Messungen bei ver- 
schiedenen absoluten Drucken ausgeführt. 
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5. Diskussion der Messungen 


Das Ziel dieses Abschnittes ist es, die gegebene 
Theorie und die ausgeführten Messungen kritisch zu 
vergleichen. Es wird versucht werden, die auftretenden 
Abweichungen zu deuten. Hinweise wird man an den 
Stellen der theoretischen Ableitung erwarten dürfen, 
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Abb.11. Abhängigkeit der Geschwindigkeit 7,, der von einem Gitter der 
skizzierten Form durchgelassenen Stoßwelle von der Stoßstärke P,, der 
einfallenden Stoßwelle.. Das Querschnittsverhältnis ö= 0,3 ist konstant 


an denen die eingeführten Voraussetzungen bekann- 
terweise nicht im vollen Umfange, sondern nur als 
Näherung erfüllt sind. 

Wie aus Abb. 9 zu ersehen ist, stimmen die Mes- 
sungen an der reflektierten Stoßwelle mit den von der 
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Abb. 12. Dichteverlauf hinter einem Gitter nach Durchgang 
einer Stoßwelle 


Theorie angegebenen Werten innerhalb der Fehler- 
grenzen in vielen Punkten überein. Es ist jedoch die 
Tendenz deutlich, daß für größere Stoßstärken P,, 
bei kleineren Querschnittsverhältnissen ö das ge- 
messene Dichteverhältnis /;, zu klein wird. Der 
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Abb.13. Mittelung der Dichte zwischen einem Gitter und einer hindurch- 

gegangenen Stoßwelle S. Rı=2,5, ö= 0,315, I3ı=1,628, T}ı gemittelt = 

1,424. Bei der Mittelung wurde die Kontaktfront vernachlässigt. Bei einer 

Stoßstärke der einfallenden Stoßwelle %,=2,5 ist in der Verengung 

zwischen den Gitterstäben W,,=1,56, BR, = 2,66 (Abb. 12). Hieraus nach 
Abb. 7a, 7b, A,=1,68 und 7;, = 1,44 % T;, gemittelt = 1,424 


Grund dieser Abweichung ist unklar, zumal bei den 
Meßpunkten für die unabhängig gemessene Stoß- 
wellengeschwindigkeit W,, eine systematische Ab- 
weichung nicht zu erkennen ist. 

Aus Abb. 9 folgt ferner, daß die an einer Quer- 
schnittsverengung reflektierte Stoßwelle sich nicht 
von der Stoßwelle unterscheidet, die an der Vorder- 
seite eines Gitters reflektiert wird. Die Meßpunkte 
für beide Fälle liegen so, daß man in beiden Fällen 
den gleichen interpolierenden Kurvenzug erhalten 
würde. Die Auffassung, den Eintritt einer Stoßwelle 
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in ein Gitter dem Eintritt in eine Querschnittsverem-\ 
gung gleichzusetzen, ist also berechtigt. Daß es hie 
bei im wesentlichen nur auf das Querschnittsverhäl 
nis ö ankommt, zeigen die eingetragenen Messung: 
von DosAansH [4], bei denen Lochplatten verwende 
wurden. 

Von dem Diagramm (Abb. 6), in dem die Geschwir 


digkeit W, der durchgehenden Stoßwelle in eine 
Querschnittsverengung gegen die Stoßstärke M, de 
einfallenden Stoßwelle aufgetragen ist, kann man al 
lesen, daß die Theorie einen weit größeren als den ex 
perimentell gemessenen Effekt voraussagt, aber nu 
dann, wenn man die Ergebnisse für die stufenartig 
Querschnittsverengung betrachtet. Hingegen liege 
die Meßwerte, die an einer Querschnittsverengung m 
allmählichem: Übergang gewonnen wurden, innerha 
der Fehlergrenzen auf der theoretisch bestimmten. 
Kurve. Zieht man noch das Diagramm Abb.7, m) 
dem die Abhängigkeit der Stoßwellengeschwindigkei 


W,, von der Querschnittsverengung ö bei konstanter‘ 
Stoßstärke P,, aufgetragen ist, zur Betrachtung hinzw 
so erkennt man, daß die Theorie die Messungen nu 
dann erklärt, wenn die Versuchsbedingungen de 
bereits in Abschnitt 3.1 genannten Voraussetzun 
genügen: Die die Strömung bestimmenden Größen? 
müssen über den Querschnitt des Kanals annäherne 
konstant sein. Diese Voraussetzung ist erfüllt 
Querschnittsverengungen mit einem allmähliche: 
Übergang vom weiteren zum engeren Querschnit 
(Abb. 6) für Verengungen mit stufenartigem Übergang 
aber nur bei Querschnittsverhältnissen Q,/Q, nahe der 
Wert 1 (Abb. 7). Bei einer stufenartigen Querschnitts 
verengung treten an den Kanten noch Wirbel au 
die in einem allmählichen Übergang fehlen. 

Im Gegensatz zu dem Vorangehenden erkennt mai 
am Diagramm (Abb. 8) für den Eintritt einer Stoßwell 
in eine Erweiterung, daß die Meßergebnisse mit de 
Voraussagen der Theorie übereinstimmen. 

Mißt man die Geschwindigkeit und die Dichte de 
durch das Gitter hindurchgegangenen Stoßwelle, s 
zeigen sich besonders bei kleinen Querschnittsver: 
hältnissen Q,/Q>,.d.h. bei Gittern mit geringem Durch 
laß, Abweichungen in Richtung einer höheren als vo! 
der Theorie vorausgesagten Stoßwellengeschwindig 
keit und Dichte (Abb. 10 und 11). Das ist erstaunlich 
denn beim Eintritt der Stoßwelle in die Verengung ist 
wie oben festgestellt wurde, eine geringere Stoßwellen 
geschwindigkeit und daher auch Dichte zu erwarten 
während beim Austritt aus den Gittern eine der Theo: 
rie entsprechende Geschwindigkeits- und Dichte- 
änderung zu erwarten steht; insgesamt sollten wi 
also hinter dem Gitter eine geringere als die theoretisch 
zu erwartende Stoßwellengeschwindigkeit und Dicht: 
beobachten können. Das Gegenteil ist jedoch der Fall 
Der scheinbare Widerspruch löst sich, wenn man be 
achtet, daß hinter dem Gitter starke Wirbel entstehen, 
die übrigens viel stärker sind als bei dem Übergang 
einer Stoßwelle in eine Erweiterung mit gleichem 
Querschnittsverhältnis ö. Die Wirbel haben zur Folge, 
daß die Dichteverteilung hinter dem Gitter am Ort 
der Wirbel ein starkes Minimum besitzt (Abb. 12), 
das in der Theorie nicht berücksichtigt wurde. Die 
Masse des Gases findet sich also weiter vorn hinteı 
der Stoßfront und erhöht dort die Dichte. Mittelt 
man die Dichte in dem Gebiet zwischen dem Gitter und 


Front der durchgegangenen Stoßwelle (Abb. 13), 
erhält man einen Wert, der unter dem theoretischen 
gt. Er wird aber bestätigt, wenn man von dem ge- 
'essenen, also gegenüber der Theorie zu kleinen 


m > 

‘ert W,, der Stoßwelle zwischen den Gitterstäben 
sgeht (Abb. 6) und aus diesem die Dichte ],, hinter 
"m Gitter bestimmt (Abb. 3 oder 8). Auf diesem 
fege erhält man einen Wert für /;,, der dem aus dem 
-perimentellen Dichtefeld gemittelten sehr nahe 


Zusammenfassung 


Es wurde untersucht, wie Querschnittsverengun- 
sn, Querschnittserweiterungen und Gitter in einem 
anal Stoßwellen beeinflussen. Die theoretische Be- 
andlung des Einflusses eines Gitters beruhte auf 
ner Kombination der Ergebnisse, die bei der Unter- 
chung der Querschnittsänderungen gewonnen wur- 
en. 
Die gegebene quasi-eindimensionale Theorie baut 
uf den Erhaltungssätzen von Masse, Impuls und 
nergie auf und deutet die Ergebnisse der Experi- 
ıente, die mit schwachen Stoßwellen bis zu einer 
toßstärke von 2,5 unter Variation der Querschnitts- 
erhältnisse an Modellen der genannten Typen aus- 
eführt wurden. Es wurde versucht, die Abweichun- 
en der Theorie von den Ergebnissen der Experimente 
deuten. Diese Deutung wird möglich, wenn man 
achprüft, in welchem Umfange die in die Theorie 
ingehenden Voraussetzungen experimentell wirklich 
füllt sind, und welchen Einfluß Effekte haben, die 
n der vereinfachten Theorie nicht berücksichtigt sind. 
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Das Bundesministerium für Wohnungsbau hat durch 
inanzielle Unterstützung die Durchführung dieser 

beit ermöglicht, wofür wir auch an dieser Stelle 
Janken möchten. 


Verzeichnis der verwendeten Symbole 


Die Symbole sind in alphabetischer Reihenfolge 
angeordnet. Die Bedeutung der Indizes ist aus dem 
Text der Arbeit zu ersehen, wenn sie im Verzeichnis 
nicht erläutert ist. 

Schallgeschwindigkeit 
Symbol für Kontaktfront 
Spezifische Wärme bei konstantem Volumen 


Von einer Stoßwelle zurückgelegter Weg 
Von einer Stoßwelle zurückgelegter Weg 


1 


Mach-Zahl im Gebiet k 
p Druck 


RB; Verhältnis der Drucke zweier Gebiete k. 
und j, Stoßstärke 

Q Kanalquerschnitte 

r Von einer Stoßwelle zurückgelegter Weg 

R Spezielle Gaskonstante 

Ss Spezifische Entropie 

5, N Sy Symbole für Stoßwellen 

t Zeit 

4 Absolute Temperatur 

u; Strömungsgeschwindigkeit des Gases im 


Gebiet j relativ zu ruhenden Koordinaten 
Te Strömungsgeschwindigkeit im Gebiet k 
bezogen auf die Schallgeschwindigkeit 
eines anderen Gebietes j 
y—=U—W Strömungsgeschwindigkeit des Gases im 
Gebiet j relativ zur Stoßfront 


Wr Stoßwellengeschwindigkeit relativ zu ru- 
henden Koordinaten 

W.; Stoßwellengeschwindigkeit bezogen auf 
die Schallgeschwindigkeit eines Gebietes j 

2 Räumliche Koordinaten 

Y=c,lc, _Adiabatenexponent 

IT; Verhältnis der Dichten zweier Gebiete k 
und 

6—=Q,/9; Verhältnis der Kanalquerschnitte, ÖFE- 
nungsverhältnis 

Yn Energiestromdichte 

Dz Energiestromdichte (dimensionslos) 

Yu Massenstromdichte 

Dyr Massenstromdichte (dimensionslos) 

0 Dichte 

Se Zeichen für Größen bei Querschnittsver- 


engung bzw. -erweiterung. 
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Experimentelle und theoretische Untersuchungen an ebenen Flächenantennen 
Von SIEGFRIED BLUME 


(Fortsetzung und Schluß aus Heft 1, 1960) 


4. Anschluß der Felddarstellungen 
und theoretischer Befund 


a) Der Reflexionskoeffizient der TEM-Welle 


Wir schließen endlich die Feldkomponente #, der 
Kugelwelle elektrischen Typs im freien Raum und die 
Komponente Z, der TEM-Welle im Innenraum der 
Sektorleitung (Abb. 1) auf der Aperturfläche r=r, 
aneinander an. Dieser Anschluß dient zur Berechnung 
der Entwicklungskoeffizienten 0}, der Reihe. Dabei 
beschränken wir uns auf die Berechnung der Reihen- 
koeffizienten 01, C} und 02. Die erforderlichen Inte- 
grationen sind im Anhang durchgeführt. 

Durch den Anschluß der Komponente ZH, der 
Kugelwelle elektrischen Typs im freien Raum — mit 
den so berechneten Koeffizienten C}, — an die Kom- 
ponente H, der TEM-Welle auf der Kugelhalbfläche 
r=r, gewinmen wir schließlich eine Formel für den 
Reflexionskoeffizienten p. Es sei verwiesen auf Glei- 
chung (101). 
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Die Funktionen G,(«) und G,(«) smd in Abb. 5 
wiedergegeben. 


b) Eingangswiderstand 


Den nunmehr bekannten Reflexionskoeffizienten p 
setzen wir in Gl. (21) und (22) ein und berechnen 
nach Gl. (24) und (25) Spannung und Strom. Das 
auf r=0 bezogene Verhältnis 


| 


ze 
r-0  1-p 


ist der Eingangswiderstand der Sektorantenne. 

Eine hiernach berechnete Widerstandsortskurve 
zeigt Abb. 6 für einen Öffnungswinkel der Kreis- 
sektor-Antenne von 2x = 90°. 

Es soll hier noch einmal betont werden, daß im 
Innenraum der Sektorleitung (Abb. 1) ausschließlich 
die TEM-Welle berücksichtigt wurde. Bei einer Er- 
weiterung der Theorie müssen Nebenwellen zur 
besseren Beschreibung der Feldverteilung auf der 
Aperturfläche r—=r, herangezogen werden. Bei hin- 
reichend großer Wellenlänge werden alle Nebenwellen 
stark gedämpft; sie werden also hauptsächlich Bei- 
träge zur Bestimmung des Nahfeldes der Antenne 
liefern. Daraus ist aber zu schließen, daß diese Neben- 
wellen praktisch auch nur die Blindanteile des Ein- 
gangswiderstandes beeinflussen. Sein Realteil wird ja 
durch das Fernfeld bestimmt. 

Nicht berücksichtigt wurde ferner bei der Berech- 
nung des Eingangswiderstandes der Einfluß des Über- 
gangsgebietes zwischen Speiseleitung und Sektor- 
leitung. 


a 


(57) 
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Des leichteren Verständnisses wegen wollen wir für 
die folgende Betrachtung eine symmetrische Antenn 
zugrunde legen. 

Es seien H und H die Lösungen der Maxwellsche: 
Gleichungen für einen sektorförmigen ebenen Flächen 
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Abb.5. Entwicklungskoeffizienten G,(«) und @,(«) des Reflexions- 
koeffizienten p nach (99) und (100) 


dipol. Wir führen dann ein neues System Z’ und H 
ein: 


VAR 


et 1 
IE regen: 


Dabei ist Z, = \uo/&o. Offenbar erfüllen Z’ und 
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Abb.6. Widerstandsortskurve der sektorförmigen Flächenantenne mit 
2x = 90° 


genauso wie vorher Z und, HZ, wie man durch Ein 
setzen erkennt. Überall dort, wo vorher etwa die 
Tangentialkomponente von Z verschwinden mußte 
verschwindet jetzt die Tangentialkomponente von H'. 
Es sind Z’ und H’ die Lösungen der Maxwellschen 
Gleichungen für das zum Flächendipol komplementäre 
System. Die Theorie des sektorförmigen ebenen 
Flächendipols erfaßt also gleichermaßen den. dazu 
komplementären Schlitzdipol (Abb. 7). Nach BooKER 
[17] besteht z. B. zwischen dem Eingangswiderstand R 
eines beliebig geformten ebenen Flächendipols und 
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Eingangswiderstand N, des hierzu komplemen- 
en Schlitzdipols der Zusammenhang 
72 
N D Ns — = . 
bei wird vorausgesetzt, daß die Flächen sehr dünn 
Vergleich zur Wellenlänge sind, daß sie eine un- 
ich große Leitfähigkeit besitzen und daß der 
eisepunkt unendlich klein ist. 


Abb. 7. Ebene Kreissektor-Antenne und dazu komplementäre 
Schlitzantenne 


Allgemein nennen wir zwei Flächen komplementär, 
nn sie zusammen eine geschlossene Fläche bilden. 
Tir denken uns eine symmetrische Kreissektor- 
tenne vom halben Öffnungswinkel « ins Unendliche 


Abb.8. Komplementäre ebene Sektoren 


ortgesetzt. Der dazu komplementäre Schlitzdipol vom 
alben Öffnungswinkel « läßt sich aber auch auf- 
assen als unendlich langer ebener Flächendipol vom 
1alben Öffnungswinkel 90°—« (Abb. 8). 


Abb. 9, Betrag des Reflexionskoeffizienten als Funktion von «& für 5 Werte 
To 


von #1, = 27 n 

Für die Eingangswiderstände NR, und Ro°—a 
dieser beiden unendlich ausgedehnten Flächendipole 
dürfen wir dann nach BookEr schreiben: 


zZ. 


RN, Roo-a = 7 5 (58) 


denn der Eingangswiderstand des zum ebenen Flächen- 
dipol mit dem halben Öffnungswinkel « komplemen- 


tären Schlitzdipols ist natürlich gleich dem Eingangs- 
widerstand des ebenen Flächendipols vom halben 
Öffnungswinkel 90°— «, da sie ja eine Antenne bilden. 
Bei unendlich ausgedehnten Flächendipolen existiert. 
nur eine vom Speisepunkt ins Unendliche auslaufende 
Welle. Am Eingang und übrigens überall auf dieser 
speziellen symmetrischen elliptischen Kegelleitung 
erscheint dann ihr Wellenwiderstand 


2, Kleine)’ 
Z Zu K (cos (90° —o)) 
90°-a = 9 K(sin(90°—o)) 


Es ist also tatsächlich 


22 
2,20 =; 


ac 
1 (39) 


wie bereits von Arır [6] aufgezeigt wurde. 

Wie lautet nun die der Gl. (58) entsprechende Be- 
ziehung für endlich lange symmetrische Kreissektor- 
Antennen? Bei Berücksichtigung von Gl. (97) und 
(59) finden wir 


x 72 (1+P» 1 Pgo_e ' 
& 90°—a& 


(60) 
Dabei ist p, der Reflexionskoeffizient der TEM-Welle 
bei dem charakteristischen Winkel« und pgo°—.. der 
Reflexionskoeffizient der TEM-Welle bei dem charak- 
teristischen Winkel 90°—«. Zwei symmetrische 
Kreissektor-Antennen mit den halben Öffnungs- 
winkeln « und 90°—« sowie der Länge r, bilden ein 
„quasikomplementäres System“: Sie füllen zusammen 
jeweils eine Kreisscheibe vom Radius r, aus. Die Un- 
symmetrie der Entwicklungskoeffizienten @,(«) und 
G;, (x) des Reflexionskoeffizienten p bezüglich « =45° 
kann man als Maß für den Fehler betrachten, den man 
macht, wenn man das Babinetsche Prinzip auf solche 
„quasikomplementäre Systeme‘ anwendet. Das Ma- 
ximum von G, («) liegt zwischen 46° und 47°, wie im 
Anhang aufgezeigt wird. 

In Abb. 9 ist der Betrag |p| als Funktion des 
halben Öffnungswinkels « für verschiedene Werte von 
x. r, aufgetragen. Mit wachsendem x - r, nimmt hier- 
nach |p| stark ab. Das ist auf Grund der starken 
Strahlungsdämpfung einer Breitbandantenne nach 
WoLTErR [5] auch zu erwarten. Aus diesen Kurven 
entnimmt man, daß bei festem x - r, gilt: 


Ip. & Ipo--al- 


Wenn wir von der leichten Unsymmetrie der Koeffi- 
zienten G, («) und G,(«) bezüglich « =45° jetzt einmal 
absehen, gilt auch 

Pu > P9o°—a 


in der Umgebung von « —45° bei vorgegebenem & - ?g- 
Dann finden wir aus (60) 


RN, No —a AA (61) 
Falls «—=45° ist, müssen wir ein Gleichheitszeichen 
setzen. Diese Beziehung zwischen den Eingangswider- 
ständen zweier Kreissektor-Antennen vom halben 
Öffnungswinkel « bzw. 90° — « ist bemerkenswert. Sie 
gestattet es nämlich, bei Messung von p, und damit 
von R, (R, liest man bei bekanntem p, aus dem Smith- 


Diagramm ab) den Eingangswiderstand Noo°—. der 
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Kreissektor-Antenne vom halben Öffnungswinkel 
90° — a. genähert zu berechnen. 

Aus der Beziehung (61) liest man ferner ab, daß 
der Eingangswiderstand zweier Sektorantennen mit 
den halben Öffnungswinkeln « und 90°—« bei unge- 
fähr dem gleichen Wert von x - r, reell wird. Betrach- 
ten wir etwa die erste „Spannungsresonanz‘: Bei 
einer Linearantenne («—-0°) wird die erste Spannungs- 
resonanz bei Xyres =7z, d.h., bei rores/X = 0,5 erreicht. 
Mit zunehmendem Winkel « wird die sog. Resonanz- 
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Abb. 10. Stehwellenverhältnis als Funktion von « für 4 Werte von 
Ton 
an, =2n Er 


zunächst immer größer werden infolge der zunehmen- 
den Breitbandigkeit der Kreissektor-Antennen. Wir 
haben aber gezeigt, daß die Eingangswiderstände der 
Kreissektor-Antennen mit den charakteristischen 
Winkeln « und 90°—« bei ungefähr demselben Wert 
von xr, reell werden müssen, d.h. aber, daß sie unge- 
fähr dieselbe Resonanzverkürzung haben müssen. Wir 
erwarten also ein Maximum der Verkürzung in der 
Nähe von «—=45°. In der Arbeit von Arrr [6] werden 
diese Resonanzverkürzungen gemessen. Die Messun- 
gen bestätigen genau unsere qualitative Aussage. 


Zwischen dem Verhältnis 


der stehenden Welle auf einer Leitung und dem Refle- 
xionskoeffizienten p besteht der Zusammenhang 


ee 

1— pl 
In Abb. 10 ist m als Funktion von « für vier verschie- 
dene Werte von x r, aufgetragen. Das Stehwellen- 
verhältnis m ist hiernach für eine Kreissektor-Antenne 
vom Öffnungswinkel 2x 90° bei festgehaltenem xr, 
jeweils kleiner als das aller Kreissektor-Antennen mit 
anderen Öffnungswinkeln. Das wird in Abb. 10 für 
vier verschiedene Werte von x r, angezeigt und folgt 
für den gesamten Konvergenzbereich der Entwicklung 
von p aus den Kurven G,(«) und G,(«). 

Wir ziehen daraus den Schluß, daß die Kreis- 
sektor-Antenne mit dem Öffnungswinkel 2x zz 90° die 
breitbandigste aller Kreissektor-Antennen ist, was 
übrigens auch aus unserer Diskussion der Resonanz- 
verkürzungen folgt. Eine Folgerung aus dem Babinet- 
schen Prinzip, angewandt auf Antennen dieser Art, 
liefert dasselbe Ergebnis (s. Arır [6]). 

Es ist zu beachten, daß das Minimum von m nicht 
genau bei «—=45° liegen wird, da auch das Maximum 
des Entwicklungskoeffizienten G, («) von p nicht genau 
bei « =45° liegt. 


Mm 


RT 
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c) Richtdiagramme 


Im Anhang wird aus der Vektorkomponente PR € 
für die Richtcharakteristik maßgebende Komp 
nente DH, des elektrischen Feldes im Fernfeld 
Kugelkoordinaten ausgerechnet. Die in dieser Form 
auftretenden Funktionen d 
Öffnungswinkels der A 
tenne sind in den Abb. 


/ 


u ze 
E 


Abb. 11. Funktionen D! («), Abb. 12. Funktionen von « nach (8 
Di(«), D;(«) nach (81), (91), (92) 


und 12 aufgetragen. Mit der Abkürzung k =cos« lau 
tet das Ergebnis [siehe Gl. (107)]: | 
ei #0 +4p e irn Do eier 


ler, 2 ee sind x 
il 
| Di («) T = 
halre) — a) 


3 c08?d [re + cs “= (1—X2) co + 


1 
RA, Men) = 
aA A 
} (4) er 7 ke? + mE (1 1)cosp] 
3 ü 
= vr, hy (xr,) 
300820 [ Re cos] a: 
+ D}(«) .— 
5 ha (er,) — 


ee Me) 
+ (&) + Pr k2+ = (1 12)cos® | | 


Falls die Antennenabmessungen klein werden im 
Vergleich zur Wellenlänge, etwa xr,<1, liefert i 
unserer Reihenentwicklung praktisch nur das erste 
Glied einen Beitrag. Die Abhängigkeit vom Winkel 9 
fällt damit heraus. Wir erhalten dann aber ein rota 
tionssymmetrisches Horizontaldiagramm. Die Ver 
tikaldiagramme gehorchen ‚unter derselben Voraus: 
setzung einem sin #-Gesetz. Die Strahlungsverteilung 
des Hertzschen Dipols ist also in unserer Lösung ent 
halten. 

Wenn ferner «—0° geht, wird wegen 1 — k?—0 die 
y-Abhängigkeit immer kleiner, was wir ja auch e 
warten mußten. 

Es wurden in dieser Arbeit alle berechneten Richt- 
diagramme auf den Wert des Fernfeldes in der Rich 
tung 9=0°, d—=90°, d.h. Richtung der y-Achse nach 
Abb. 1, normiert. 


Die Abb. 13a—c enthalten die berechneten Ho- 
sntaldiagramme (yz-Ebene) einer Kreissektor- An- 
e mit einem Öffnungswinkel von 2x =90°. 

Wir finden für r,/A=0,250 ein fast rotations- 
nmetrisches Diagramm. Bei r,/A =0,382 ist die 
-ahlung der Antenne in Normalenrichtung (P =0°, 
90°) stark gegenüber der Seitenstrahlung ge- 
chsen. Bei r,/4=0,500 beträgt das Verhältnis von 
rmalenstrahlung zur Seitenstrahlung ungefähr 1:0,6. 
sses Verhältnis wird bei größeren Werten von ry/A 
der kleiner. Die für r,/A=0,732 und r,/A =1,000 
rechneten Diagramme zeigen, daß dort die Seiten- 
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Beide Vertikalkennlinien haben bei 9—=0° bzw. 
#=180° eine frequenzunabhängige Nullstelle. Eine 
Kreissektor-Antenne könnte also eventuell bei einem 
Breitband-Minimum-Peilverfahren Verwendung fin- . 
den. 

Die berechneten Vertikalkennlinien stimmen mit 
den von Arur [6] bei etwas anderen Werten von ry/A 
gemessenen Diagrammen befriedigend überein. Bei 
seinen Messungen benutzte ARLT eine symmetrische 
Antenne. 


30° 


Abb. 13a. 


trahlung die Strahlung in Normalenrichtung bereits 
ibertrifft. Wegen der bei r,/A—=1,000 bestimmt un- 
‚enügenden Konvergenz der Reihenentwicklung ist 
las letzte Diagramm allerdings nur qualitativ zu 
werten. 

Die berechneten Werte des Verhältnisses 


Normalenstrahlung 
Seitenstrahlung 


1 dem betrachteten Frequenzbereich stimmen gut 
iberein mit in der Literatur angegebenen gemessenen 
Werten (s. [1], [5], [6))- Insbesondere zeigen Mes- 
sungen von ARLT [6] sowie von Brown und WoonD- 
wARD [1] das starke Anwachsen der Seitenstrahlung 
bei großen Werten von rg[A. Horizontaldiagramme 
ebener Flächenantennen sind in den genannten Ar- 
beiten nicht veröffentlicht worden. 

Die berechneten Vertikaldiagramme einer Kreis- 
sektor-Antenne mit einem Öffnungswinkel von 2x =90° 
in den Abb. l4a—c weisen die sektorförmige ebene 
Flächenantenne noch einmal als eine Breitband- 
antenne aus. Die Richtcharakteristiken in der xy- 
Ebene zeigen bei wachsendem 9/4 Einschnürungen, 
die aber nicht auf Null einziehen. 

_ Die Diagramme in der x2-Ebene gehorchen im 
‚gesamten betrachteten Frequenzbereich dem sin d- 
Gesetz. 


Abb. 13b. 


Abb. 13a—c. Horizontaldiagramme der sektorförmigen Flächenantenne mit 2« = 90°. 


Abb.13c. 
yz-Ebene (nach Abb. 1) 


d) Stromdichteverteilung 


Da die Stromdichte direkt proportional ist der 
magnetischen Feldstärke, können wir aus der Feld- 
komponente H, nach Gl. (22) eine Aussage über die 
Stromdichteverteilung auf der Antennenfläche und 
auf der Grundebene machen. 

Wenn wir die Stromdichte längs eines Kreisbogens 
auf dem Sektor »—b auf die Stromdichte in der Mitte 

—q dieses Kreisbogens beziehen, dann finden wir 
aus Gl. (22) 
_ ob 


Ve — 2 


Für »=b folgt aus den Transformationsformeln (1) 


2) 
u=a) 


u? = a? cos? d®, 
so daß man also bei Beachtung von bja cos « erhält: 


29 
i@=0}) 


sin& 


\eos% — cos? 


Wir führen den komplementären Winkel y=90° —d 
ein und finden 


ip) 


Ip EIKE sin? a. — cos? y 5 


sin 


Die Wurzel nimmt reelle Werte nur dann an, wenn 
der Winkel gerade den Sektor überstreicht, d.h., 


a ur 2 N yet 


0° —- a <y<HN+a 


ist. In Abb. 15a ist diese relative Stromdichte als 
Funktion des Winkels y für einen Sektor mit 2x =90° 
aufgetragen. Die Stromdichte wird also, wie zu er- 
warten, singulär an den Kanten des Sektors. 

Für die Stromdichte längs eines Kreisbogens in der 
Grundebene »—=0 bzw. #=90°, wieder bezogen auf 
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Unter Zuhilfenahme der Transformationsformeln 
erhält man hier 
Hoyer 1 £ 
je =) Y cos?« + sin?« cos? 
Diese Funktion ist in Abb. 15b wiedergegeben 
29-904 


Abb. 14a 
Abb. 14a—c. Vertikaldiagramme der sektorförmigen Flächenantenne mit 2x = 90°, 


GE 
Yv 


Abb. 15a. Stromdichteverteilung längs eines Kreisbogens auf der scktor- 
g förmigen Flächenantenne mit 2x = 90° 


u 
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Abb.15b. Stromdichteverteilung längs eines Halbkreisbogens auf der 
Grundebene (yz-Ebene nach Abb.1) bei einer sektorförmigen Flächen- 
antenne mit 2x = 90° 


die Stromdichte im Punkte u=a auf diesem Kreis- 
bogen, finden wir aus Gl. (22) 
5 VAR N. 


ju=a u’ 


Abb. 14b 


/ 

» 
1,0 a 
Ze duun 
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Abb. 14c 
zy-Ebene; --------- xz-Ebene (nach Abb.1) 
5. Anhang 


a) Berechnung der Entwicklungskoeffizienten 


Wir berechnen zunächst die Koeffizienten ( au 
dem Anschluß der Feldkomponenten Z, auf de 
Kugelhalbfläche r=r,. Aus Gl. (21) und (54) folg 
bei geeigneter Umformung 


ei #ro +p-e7*n 1 
ET Ye) @—») 


dkKs(v , 
Cat 100) Ser) Ku a) 


Integration über » von O bis » liefert 


E er * ro +4pe7i*To 


2, To 


ji dv’ 
F y® —y’?) (a? — v2) 
3n+1) (62 
ö n 
Or [M-1er)— -hrlero)| x 


x Ko) Kn(u), 


da ja rechts alle ungeradzahligen X‘; (v)-Funktionen 
für v=0 verschwinden. Im Integral auf der linker 
Seite haben wir dabei die Integrationsvariable mit » 
bezeichnet. 

Wirmachen nun Gebrauch von den Orthogonalitäts 
eigenschaften der Lam&schen Produkte und zwar ir 
der folgenden Form (s. HEıne [10] oder Hogsox [11]) 

Es seien F° und E%Y zwei Lamesche Funktionen 
derselben Klasse vom Grade n und »’. Man kann dann 


en, daß das Doppelintegral 


a 
» [au {BOB RE WE x 
b (63) 
we—»? 
x S = 55 - 
V& — ») (a —»2) Yu? — 8?) (a? =) 
‚schwindet außer für n=n’ und s=s'. 
Wir beschränken uns auf die Berechnung der 
effizienten 01, 01 und 0%. 
m—=1 
Für n=1 findet man nur eine Funktion der 


sse K‘, (s. Hogson [11]) 
Kıb) =», 
64 
Kıl) =u; sn 


i geeignet vorgenommener Normierung. 
‚Wir multiplizieren Gl. (62) mit 
u —»? 
1 — 9) (a =) Ya — 09) (@ 
d integrieren über den gesamten Wertebereich von 
d u. Bei Berücksichtigung der Orthogonalitäts- 
ation (63) erhalten wir 


=) 


xrg 
b @ r (u 2) 
Sur NE e) e 
/ ey du Ve) @ Ve ee) 
 ee- er . (65) 
ılyı 1,1 
162 IR @ero) el ro)| x 


ev? (u? — 9?) 
»») (a2) Vu2—b) ( 


b a 
IKK: 5}: 
0 b 


E lösen zunächst das Doppelintegral auf der linken 
ite von Gl. (65). Dazu substituieren wir 


ec: (66) 


it der Bezeichnung 


I (67) 
a 


gt aus dem im Integranden stehenden unbestimmten 
Be! 


Fer: er 2 —y’?2) 


so bis auf die Konstante 1/a das elliptische Integral 
ster Gattung. 
Wir führen jetzt 


ern 1er2) ’ 


(68) 


So A_Bı-Rrd 
ihr 2?) v2) 


Is neue Variable ein. Speziell ist 


u - dt’ 
a ya-r9ı ker?) 
0 r 


as vollständige elliptische Integral erster Gattung, 


as wir im folgenden mit X bezeichnen wollen. 
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Die Umkehrfunktion des elliptischen Integrals 
erster Gattung ist die Jacobische elliptische Funktion 
t=sn(v;k). (69) 
Dafür wollen wir in Zukunft einfach sn v schreiben. 
Da nach (66) aber t=»v/b ist, hat man jetzt die ellipti- 
sche Kegelkoordinate » durch eine Jacobische ellipti- 
sche Funktion dargestellt. . 
Wir erhalten für das auf der linken Seite von 
Gl. (65) stehende Doppelintegral zunächst 


b a 
3 era _n) & af: 
et du »2) (a? — 9?) Ya b2) (a? — u2) > 
0 


ee 


AR PER (70) 
= [esmvdv Ii IB 
a J Ve — 2) (a — 1?) 
0 b 
Bor x d 
— „| vontodo | ——. wir ; 
& J Vu? — 8?) (a? — u?) 
0 b 


Auf die Darstellung der elliptischen Kegelkoordinate u 
durch eine Jacobische elliptische Funktion wollen wir 
hier verzichten, da die auftretenden Integrale über u 
leicht ausführbar sind, wenn wir u? als neue Variable 
einführen (s. Rekursionsformel bei GRÖBNER-Hor- 
REITER I). 


ud 7 
- —. 71 
f Var - 2 — b2) (a? — =) 2 Sn 


Mad a Nu. 2 > 
ee Zelte. 2) 
Die für das Rechnen mit Jacobischen elliptischen 
Funktionen notwendigen Formeln sind, soweit wir sie 
hier gebrauchen werden, im Anhang zusammen- 
gestellt nach MırLne-Tuomson [18] und GRÖBNER- 
HOorFREITER [17]. Die bestimmten Integrale in dieser 
Formelsammlung wurden selbst berechnet. 

Es ist nun 
K 


[on ode — ji (dnv— ken ol, = 
€ 0 


0 


K 
u. [in (ano — kon») dv 


(73) 
0 
K 
- I mfi-®-— ; [mtano—konvdo, 
0 
K 
Fesntodo 
_ [, envdnv 1+ ze 
= re 9,5 a (dnv— kenv)| 
K IE 
1 [envdnedo— "ya [Intnv—keno)do (74) 
0 0 
Een ——. 'ı 
a Imwne-ronnan 
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Wir setzen (71), (7 


b a 
! e, E (m 


2), (73) und (74) in (70) ein: 


uva — 9?) % 
v2) (a? — »?) \(u?—b°) (a°— 1?) 


% en. 
oe) 


x 


K 
n I N 


4% nyı-® 


2k® 


Ra fa (Anekenn)dı) 7. 


BT 


x Zar um-esi hi=k®—.,— 


Das allein ER vorkommende bestimmte Integral 
K 
[In (dnv — kenv)dv 
0) 


fällt dann heraus, und wir erhalten 
2 2 
d fi ab ur (u? — »?) 
ß 12 u e X b2 — »2) (a? — »2) Ve — 82) (a — 


er = a 


Zur Bestimmung des auf der rechten Seite von Gl. (65) 
stehenden Doppelintegrals machen wir Gebrauch von 
einem Doppelintegraltheorem für Lam&sche Produkte 
(s. Hrıse [10] oder Hogson [11]). Wir wollen es hier 
folgendermaßen anschreiben: 


Mit der Substitution (68) läßt sich das Integral 


ES (y 22 
[ii Lin Te — 92) (a? — 9?) 


durch RN zurückführen auf 


K 
f(&* —n*sn?v)dv. 
0 


Ebenso läßt sich das Integral 


(75) 


(76) 


dv 
J v (»)] V® = »») (a? = „) 


durch Reduktion zurückführen auf 


K 
S(0* — 
0 


(77) 


D* sn?v) dv. 


Wir setzen 
& md: en 


ad 
a 
9=ua9*, 
= 9*, 


(78) 


Die Aussage des Doppelintegraltheorems lautet dann: 


b a 
[Er ) Brta)P - (u? — »?) 
dv | du 
J 4 el »2) (0? — »2) Yu — b?) (a? So (79) 
— z n9- E98). 


ee 
32 Re ir: Wen y 


, 
Zeitschrift fi 
angewandte Ph 


Das Doppelintegral auf der rechten Seite von GI. ( 
ist ein Doppelintegral eben dieses Typs wegen K} (v) 
und K} (u) = 


Wir reduzieren gemäß (76) 


> A 


2 
Pre 


Wir reduzieren gemäß (77) 


b K 
dv bt 
f» ren u; IESZZ 


0 0 
K 
[a eh +1) [sn’vdo, 
0 
bBa? 
0=— > 
®=-—- a1 +2). 


Nach (79) erhalten wir also 


b a 
ur v* (u? — 9°) | 
Ni eye — #)(a® — 9°) Via? — 8°) (a? — an ( 


a WARE 


BIT? 


Die Integrale (75) und (80) setzen wir jetzt in Gl. (6@ h 
ein und finden I 


ai eRLTDOL D!/a® (8 j 
er, 1 2 
\ hö(x ro) en | i 
mit | 
1 3 Q1 | 
3: S 
n=3 Mi 


Im Falle n=3 findet man zwei Funktionen € 
Klasse X}, (s. Hogson [11]) 
Kb)=n+A», 


Kb) ="+ Br 
Kz(u) 


=uw+4An, Kalu)=u’+ Bu. 


Hier ist wieder der Koeffizient der höchsten Pote 
von » bzw. von u gleich 1 gesetzt worden. Es ist da 


I+k V 60k® 
A=a: 10 + 16 - ar 4]. 
‚I1+® 608 
B=alt | 16 Bi. 


_ = und — = sind in Abb. 12 als Funktionen vo! 


wiedergegeben. 
Multiplikation der Gl. (62) mit 


Ki(v) Ki(u) nn 
V@® — ») (a? — 9) Var — )@ —) 


und Integration über den gesamten Wertebereich v 
v und u liefert bei Berücksichtigung der Orthogona 


\ 
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Br r 
relation (63) Nun ist 


ern L perirn K K 
f envdn vsn?vdv = [sn?v]X — 2 [ envdnvsn2vdv, 
0 0 


+Av)(u? + Au) (u? — v2) also 
V@ =) (0°) 3°) (@®— 8) K i 
f envdnvsn?vdv =>. 
x — a. ° 
; Ve? — v2) (a — v?) Damit erhalten wir 
Hero) — haar] x SR: 
il 2(% 70) a 0) [eswvav = X 5 R ne | In Yı—R 
>+Av)? (u? + Au)?(u2 —R) u 

' 3 H. du! — REN EN, 
; h. Le (a? —»?) Y(u—0?) (a? n EISEN 
en zunächst wieder das auf der linken E Au Ri)2 x En 
e von Gl. (84) stehende Doppelintegral. u 2k° 
Zur Berechnung der bestimmten Integrale über u K 
‚ren wir u? als neue Variable ein. Man findet (siehe x ii In (dnv — kenv)dv. 
kursionsformel bei GRÖBNER-HOFREITER [17]) 
‘ Hi w-+Au® In (73), (74) und (88) sind die Integralumformungen 

Vie — 2%) (a — 12) du 85 für die in (87) auftretenden bestimmten Integrale 

Rn, (85) über vangegeben. Wir setzen diese Integralumformun- 
d (+ KR)? Ba— "Ba +AT(1+K)ar, gen in (87) ein und berücksichtigen die Werte für J, 
SB RP 4 4 und J, aus (85) und (86). Man stellt dann fest, daß 
das allein noch vorkommende bestimmte Integral 
| Br Naar 4”..‘ '(8e) 
e Veen 4 = F 
f[ In(dnv — kenv)dv 
er Substitution (68) erhält man für das in Frage 0 

hende Doppelintegral zunächst herausfällt. Wir ordnen den Rest nach Potenzen von 


Aja® und finden schließlich folgenden Wert für das 
auf der linken Seite von (84) stehende Integral: 


fa W+Av) (u? + Au) (u? —»?) 
di a I ee 


b @ 
Be far fau \y EN er 
} Ve y@ 7) Fk V(ß® — »*)(a® AUG b2) (a? — u?) 
Urt ee (Av) a (89) 
i (87) ’ V@- b2 — y’2)(a2 — 73) 
"Ye a? — »?) nn ( v . a—y 
ba? 2 
ER: -— II (Te + + (94 Be ee (4) | 
7 (A sn?v — b?sndv)vdv + 
Ü) K Bei der Berechnung des auf der rechten Seite von 
- 4 (J, sn v — J,b?sn?v)vdev. Gl. (84) stehenden Doppelintegrals machen wir wieder 


Gebrauch von dem zitierten Doppelintegraltheorem. 


g Gemäß der Vorschrift (76) reduzieren wir 


ir integrieren partiell 


n d 
vsnd _[, modnvsnv 3 1+K2 iR f Su 12 ea 
Pr vdv lv TE Se m venvdnv) + (v® + A») Ver — wie =») 
3 (1+22)2 K 
nn (3 ein (inv— ken) — = [Won + 2ABtanto + A? b° sn? »)do 
K 
a kenn), — Er 4 1+M af bb 8 
R =| en rer 
— ga | vn von? var — ö 
4Kk? K 
0 22 6 4 2 2 
-54% [on vanvav — ö 
) Wir lesen ab: 
3 (1+ 22)? 1 
Dr x 4 1+22 2 
2K3 ge 4 
Su y a a 
x [in (dn®» — kenv)dv. 8 (1-F AR)? Bis 41+B 
R NE ee ar Tee ea 4 
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Wir reduzieren gemäß (77) 


b 
j v2 (y + Av)? 


0 


dv 


(ar) 


Y (b? — »2 


R 
[ww sn®v + 2 Ab6 sn» 4 A?btsntv) dv 
ad 


ME=rs); 


a lat a5 —; 
K 

a ee e) 
un A ze) 
be 16 (IHM 9 104 IHR 

Klar; O9 15 M 

6 3 212 6 

b° 16 (I+R2)? a® 6 


Ir Ze 
6 kt a 5h2 


E& 5 BB 8 


HA 


K 
12 
Mr 2 nn [n®vav 


0 


und lesen wieder ab: 


5 
alkt 35 ks 
SE 
En Alan 
ad BE 108 
a 105 es 


6 el RE } 
ae 3 mM 


0,63 


— Abs: 


16 (1 
bi 


® 


Abt 7 


Nach (79) erhalten wir schließlich, geordnet nach 
Potenzen von A/a?, 


b a 
; P HA + An) ) 
iv | d = —E 
/ Br 1 “ \ Y(b? — v2) (a? — v2) (u? — 2) (a? — u2) } 
1) b 


IT 


A 
370 70 un {1504 4 36%2 (1 + %2) z P* 
+[24 (1 + k2)? + 382] Es 2% 


(90) 


+56(1 +2) (2 =) +3(3 =) }- 


Die nunmehr bekannten Integrale (89) und (90) 
setzen wir in Gl. (84) ein. Wir finden 


i#r, —ian 1/47 
Gel > en (91a) 
ns ha (x ro) — SEE ha(er,) 
Krg 
mit der Abkürzung 
De nn 
15kt + 36R2(1—+ a — 
A? (91b) 


+0+119) 5 et) (=) 
+24 MEReERTE: ) a) +3:(5)" 


Die Lameschen Funktionen K}(v) und K?(v) unter- 
scheiden sich nur durch die Konstanten A und B 
voneinander bei sonst gleichem Aufbau. Der Reihen- 
koeffizient (3 hat also auch genau dieselbe Form wie 


De a ar 


angewandte P 


C!, nur mit einem B überall dort, wo bei O1 ei 
steht. Es ist also 


Dia? 


3 
Alan) 


er en j 
Kr, 


02 — 


ha(xro) 
mit der Abkürzung 

35 7ER +kt+ 
Eh 53ER (14 Re) a 


Di = 


+02) ee 


+[24( ee ) ner, .) +) 


b) Berechnung des Reflexionskoeffizienten 


Die Reihenkoeffizienten OS, enthalten den bislam 
noch unbekannten Reflexionskoeffizienten p. Zu‘ 


ner Bestimmung schließen wir jetzt die Feldkomponer]) 


ten H, der TEM-Welle im Innenraum unserer Sektor 
leitung und der Kugelwelle elektrischen Typs im Auß 
raum der Leitung auf der Aperturfläche r =r, ane 
ander an. 

Aus Gl. (22) und (55) erhalten wir bei geeigne 
Umformung und bei Berücksichtigung von 4, =By/2 


ei*To —peri*n 1 


Kr IK — 22) (a? — »?) 
co $(n+1) 


> > 08 hi(erı) 


n=1;3;5 s=1 


Wir integrieren über » von 0 bis b und über « von I 


bis a: 


ei ro — nette 


0 


k)K 
Krg 
eek) (93 


=D rar rim 


n=1,3;5 s=1 


Auf der linken Seite En Gleichung ist mit K wied 
das vollständige elliptische Integral erster Gattung 
mit dem Modul k=cos« bezeichnet. | 

Wir beschränken uns wieder auf die Reihenglied 
n=l’und n=3. 


Es ist nun 
Kı(lb)=b, 
Kib)—=b3 + Ab = ak(R + 5). 


Ki =W4+Bb-arle 4), 


[sw )\au=aıE, 


[rwan-%(i+e+24)a-m, 
D 


[Rau = 
b 


i+e+2z)a-m. 
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ir setzen dies in Gl. (93) ein und erhalten 
u ei*To ale 
0 rg 


onen) 


| ala Rt )a+R (94) 
;C} hi(xr,) a’ n # 
Ei oehler,)a” ele+ alle at 


‚ Cl und 0% sind nach Gl. (81a), (91a) und (92a) 
ann. Um unsere Endformel übersichtlich schrei- 
en zu können, führen wir eine weitere Abkürzung ein: 


k(1-+k) 


ı@)=Di Tepe (95) 
e[e+ Z)(1+#+25)a+M 
ha) =D} a 
B B 
ler ir +2) rn R 
(«) =D; iK i (97) 


= eirn x 
s(n+1) Aa 
1-5) 2, So) nero) 
n=1;3 s=1 hn-ı(&r) = Mn(&ro) (98) 
5 Far) 5 
14H), 2, 8) m. 
1:31 kn-ılaro) — een) 


Die Koeffizienten G1, G4 und 63 hängen wegen k = cosa 
ur noch vom Winkel « ab. Wir fassen zusammen: 


+ k 
(=, (99) 
@;(&) = G3(«) + 65(«) 
A Aa 
hr TarH +12 =) x 
BR 151443612 (1419), + [24014 R2)243810] x 
x[#+ (ur 2,)a +k) 
A\2 Ans? A 
x[5) +56) +35 (4) (100) 
B B\ 
35 rare 2) |x 
32K 5p43618(1 1432) +24 (1429)2+3812] x 
B B 
x[#+ ar )arR 
B\2 B'3 Bi’ 
x (2) +56(14+1°) (2) +35 (=) 
‚Wir können jetzt endgültig schreiben: 
1 
1-5), na) a) 
n=1;3 hn-ı(&r) — SEE hn(ro) 
Bin er (101) 
14:7 G„(®) m 
n=1;3 hn-ı(&r0) — Pr hu(ero) 


p ist natürlich eine dimensionslose Größe. 
2.f. angew. Physik. Bd. 12 


In Abb.5 sind die wichtigen Funktionen 6, («) 
und G,(«) aufgetragen. Die relativ stärkste Abnahme 
der Entwicklungskoeffizienten G,(«) und G,(«) liegt 
offenbar bei BER) vor. Für einen Kreissektor mit 
dem Öffnungswinkel 2x 490° wird unsere Reihen- 
entwicklung also am besten konvergieren. Um die 
Bar der Reihenentwicklung für große Werte 
von xr, zu gewährleisten, muß man weitere Glieder 
n=5;7;9;... mitnehmen. Die Lameschen Poly- 
nome K%,K$,Kö,... sind aber noch nicht bekannt. 
Bei der Berechnung dieser Polynome wird man auf 
Gleichungen vom Grade 3 (n+1) geführt. 

Wir wollen die Lage des Extremums (mit horizon- 
taler Wendetangente) der Funktion @, (x) bestimmen. 


dc)  d R =. dk 


de dk\4 K(k)) de 
dK(k) 
a EDER 
u, Bu de 
Es ist 
i 
Km=- | ——r 
’ Ya-Bua-%?) 
also 


k2? dt 


Be an 


SH + 


(1 2)3 
Wir substituieren 
v2 = a-®)r 
ln 
und finden 
dK (k) (1-K)+Rert2 , 
er Soma )4 nl ra a 
Es ist also 
dK (k) 7 E(k) 
er ATS) 


E(k) ist das vollständige elliptische Integral zweiter 
Gattung mit dem Modul k=cos«. Wir setzen 
46, («) 


da 


und finden dann folgende Bestimmungsgleichung für 
die Lage «, des Extremums von G,(«) 


E (cos a.) — (1 + cosay — 2 cos?«,) K (cos &,)= 0. 


Eine Überschlagsrechnung liefert 


46 <un< 47°. (102) 
Bei rotationssymmetrischen Kegelantennen lauten die 
unseren G',(«) entsprechenden Funktionen (s. BORGNIS- 
Papas [3]): 


Dabei ist « der halbe Öffnungswinkel der Kegelantenne 

und P, (cos «) eine zonale Kugelfunktion vom Grade n 

(n=1;3;5..... Die beiden ersten Entwicklungs- 

koeffizienten (n—=1 und n=3) weisen ungefähr das- 

selbe Verhalten auf wie unsere Koeffizienten G,(«) 
6 


und G,(«). Allerdings liegt das Maximum des ersten 
Entwicklungskoeffizienten hier bei «, 36°. Die hö- 
heren Entwicklungskoeffizienten (n=5;7;9;....) blei- 
ben bei rotationssymmetrischen Kegelantennen klein. 


Qualitativ ist dasselbe Verhalten bei den in dieser 


Arbeit behandelten Kreissektor-Antennen zu erwarten. 

Um den Konvergenzbereich der Entwicklung des 
Reflexionskoeffizienten p in (101) unter alleiniger 
Berücksichtigung der berechneten Werte von @,(«) 
und G; («) für =45° festzustellen, wurde der Ausdruck 


hn(z ro) 


n 
hn-ı(%&r,) — > hy(zro) 
“Tg 


als Funktion von x, fürn=1,n=3 und n=5 be- 
rechnet. Wenn man annimmt, daß G,(45°) kleiner 
wird als G,(45°), kann man mit befriedigender Kon- 
vergenz bis x r,7%3,8 rechnen. 


c) Berechnung der elektrischen Feldkomponente Ey. 


Wegen sp 
uU= —_— 
ör 
ist a 
B,=gradgu=t | er (103) 


Wir fanden folgende Entwicklung für die Komponente 
P, des Hertzschen Vektors [s. Gl. (50)]: 


3(R +1) 


2125 Sayen 


n=1;,3;5 s=1 


Hn+,(«r) K,(v) Kr (u je 
Bei Berücksichtigung von (53) erhält man 
=, Maler)| x (104) 
x Kim) Kin). 
Wir beschränken uns wieder auf die Reihenglieder 
a—l und n=3. Esist 
Kı(v) Ki(u) = vu, 
Kz3b) Kz(u) =’? +A? va +Aval® +), 
K3(v) Kä(u) =v?u? + Byu + Bvu® +0), 


mit den Konstanten A und B aus (83). Aus den 
Transformationsformeln (1) gewinnt man 


vu=bacosd, 
„+ u? =a?[(l + %k2) — (sin?p + k? cos? p)sin?9]. 
Dabei wurde wieder k =b/a —= cos « gesetzt. 
Die oben angeschriebenen Produkte K5(v) - Ki,(u) 
können wir hiermit in Kugelkoordinaten 9 und 
umformen. Das Ergebnis dieser Umformung setzen 


wir in Gl. (104) ein. Aus GI. (103) erhält man dann 
schließlich 


IE =A0r, [18 (#r) — ee ba——  (c0sd) + 


+0 Ren) — in en hklun)] u 

x {b3a? RO; c0sd+ A ba? cosd X 
x [1+%®— (sin?o + k? cos? ) sin?9]} + 

+08 11 (er) — om (er) 5 

x {b3 a? cos?d + B?ba cosd + Bba3cosd x 

x [1 +%?— (sin?p + k2 cos?) sin29]}. 


(105) 


x 
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Uns interessiert das Fernfeld (<- r>1). Wir mache: 
Gebrauch von der asymptotischen Entwicklung de 
Hankelschen Funktionen (s. SOMMERFELD [14]) 


2 kn f 
Hu+,(# r)—= = eig mtl) 


Es ist also asymptotisch 


a eier 
hu(ar) = —imw+D, 
Kun=. 


lich folgende Formel: 


ÄAT, —ixn ir 
By“ =B,- en Er sind x 
a) Er 
x 1-.Die)- \ - 
ha(&r,) — un Men) 
: 3c0s?d v+5 + an = (1— k2) cos? 2 + 
+D;(a) hi(zr,) — 
A A 
+) ++, 1-2) cos g| 


3 
Far hy (ro) 


3cos?# [# + E _ 


alt 


3 
= er, hy(zr,) 


5 Zu-) cp] - 
+D;(«) 


hi(zr,) — 


ar GG 2 — k?) cos? v) 


Dabei wurden Glieder mit 1/r? gegen Glieder mit1/s 
vernachlässigt. 
d,) Formelsammlung 


Für die Jacobischen elliptischen Funktionen 
sn(v;k), en(v;k) und dn(v;k) wollen wir in aller 
Formeln einfach snv, cnv und dn» schreiben. 


sv ten®v—=1, 
dn?®v + ksn®v=1. 


Integral erster Gattung. 


dx 
KR=Kik)= —— 
(m ı ya== a2) (1 — ea) ” 


nn (sn v) = envdnv, 

d 

7 (env) =— snvdnv, 
5:7 (dn®e) =— k’snvenv, 
[nv envdav=— ginn, 


. Band 
ft 2 — 1960 


[snvdnvdv=—eonv, 
[en vdnvdv = snv, 


env dnv sn? 3» 
sn?odv — 


m) 
—2)(1 +12 
a [era 


w—'3. N n—4 ’ 
- ne [m vdv, 
4 In(dnv— keno), 
envdnv | I1+E] 

2E Tem 
envdnvsn?v 3 1+% 


(n+1) 


snvydv = 


ensydv = (dn®»— kenv), 


{-289 — ’ 
sn’vdv= 42 ER envdnv + 
2\2 
+); -. sn (dno— konn), 
K K 
3 1+R 1 
a Bere de en 
sn vdv=, 18 IE vdv a Jar. 
0 () 
2 3 F 
8 (1+%) 
6 # : = 
sn vdv=|-, Mm ZauE. vdv 
) 
: K 
4 
sl 
() 
[% k2 E 
16 (1-+ R2)3 104 af 
8 E g SER 
sn vdo=|, 1 05 sn®vdv + 
0 
K 
1:65 8 (I+M) 
OT BB fer 
0 
Teil C 


Experimentelle Untersuchungen 
an ebenen Flächenantennen 


1. Widerstandsmessung 


a) Methode zur Messung 
des Antennen-Bingangswiderstandes 


Den Versuchsaufbau zur Ermittlung des Eingangs- 
iderstandes zeigt im Prinzip Abb. 1. 

Von einem Rohde & Schwarz-Sender führt ein 
oaxialkabel zu einer gewöhnlichen geschlitzten Meß- 
eitung mit gleitender Sonde. Diese ist mit einem 
nzeigegerät verbunden, welches direkt das Ver- 
Hältnis Unin/Umax der stehenden Welle auf der Meß- 
eitung anzeigt. Der Innenleiter der Meßleitung führt 
er ein als Koaxialleitung ausgebildetes kurzes An- 

chlußstück zu der Flächenantenne. Der Außenleiter 
a Meßleitung steht über den Außenmantel des An- 
schlußstückes direkt mit einer quadratischen Alumi- 
niumplatte in Verbindung, die eine Seitenlänge von 
2 m besitzt. Der Fußpunkt der ebenen Flächenantenne 
liegt genau in der Ebene dieser Grundplatte. 

Aus der Lage des ersten Spannungsminimums auf 
er Speiseleitung, vom Antennenfußpunkt aus ge- 
schen. und den abgelesenen Werten von Unin/Umax 
de in bekannter Weiese der auf den Antennenfuß- 
punkt bezogene Widerstand aus dem Smith-Diagramm 
abgelesen. Die gesamte Leitung zwischen Sender und 
Antenne besaß einen Wellenwiderstand von 60.0. 
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b) Widerstandsortskurven ebener Flächenantennen 


Die so gemessene Widerstandsortskurve einer 
Kreissektor-Antenne mit einem Öffnungswinkel von 
2&—=90° zeigt Abb. 16. Sie ist zu vergleichen mit der 
berechneten Widerstandsortskurve in Abb.6. Der 
Realteil wird demnach durch unsere Rechnung einiger- 
maßen gut wiedergegeben. Die starke Abweichung in 
den Blindanteilen ist, wie bereits näher ausgeführt 
wurde, auf die alleinige Berücksichtigung der Haupt- 
welle. zurückzuführen. 

Die Feldtheorie der ebenen Sektorantenne zeichnete 
den 90°-Kreissektor durch ein Maximum an Breit- 
bandigkeit gegenüber allen Kreissektor-Antennen mit 
anderen Öffnungswinkeln aus. Ist es möglich, eine 
ebene Flächenantenne anzugeben, die den 90°-Kreis- 
sektor an Breitbandigkeit übertrifft ? 
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Abb. 16. Widerstandsortskurve der sektorförmigen Flächenantenne mit 
2% = 90° 


Die Analyse des zur Ausstrahlung führenden Wel- 
lenvorganges ergab, daß der Eingangswiderstand von 
Antennen mit ortsunabhängigem Wellenwiderstand be- 
stimmt wird durch die Reflexion der auslaufenden 
Welle an der Inhomogenitätsstelle, die unsere spezielle 
elliptische Kegelleitung vom freien Raum trennt. Bei 
ebenen Flächenantennen mit anderen Formen wird 
aber noch eine Reflexion infolge des jetzt ortsabhängi- 
gen Wellenwiderstandes hinzutreten. Wir über- 
tragen nun in geeigneter Weise eine Aussage der 
Leitungstheorie auf unser Problem: Es ist zu vermu- 
ten, daß insbesondere Antennen mit exponentiellem 
Wellenwiderstandsverlauf günstige Breitbandeigen- 
schaften aufweisen werden. 

Wir wollen mit r(6) den Ortsvektor der Rand- 
punkte einer ebenen Fläche bezeichnen, die zu einer 
Geraden symmetrisch ist. Symmetriegerade und Orts- 
vektor schließen den Winkel ö ein. Dann soll der orts- 
abhängige Wellenwiderstand von der folgenden Form 


sein: 
_*(6) . 


Z=(Ge'') +0,. 
Wir setzen [s. Gl. (27)] 


ze K (cos ö) 
ZEIT (ind), 

und erhalten als Bestimmungsgleichung für r(ö) 
r (8) 


2=QerW) +0,. 
6* 


K (cos ö) 


307 K (sin ö) 
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Es ergeben sich flammenförmige Kurven r(ö) mit ver- 
schiedenen Öffnungswinkeln bei r=0 je nach Wahl 
von O, und C, (C,+C,>0). Rein theoretisch ge- 
sehen, haben alle Formen eine unendlich lange Spitze. 
Auf eine Untersuchung dieser Exponential-Flächen- 
antennen wurde in dieser Arbeit verzichtet. Ins- 
besondere müßte festgestellt werden, welcher Öft- 
nungswinkel am Antenneneingang jetzt zu optimalem 
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Abb. 17. Widerstandsortskurve der halben lemniskatenförmigen Flächen- 
antenne 


Breitbandverhalten führt. Außerdem müßte der Ein- 
fluß endlich langer Spitzen untersucht werden. 
Stattdessen wurde eine halbe lemniskatenförmige 
Flächenantenne durchgemessen. Die Lemniskate be- 
sitzt einen Öffnungswinkel von 90° und wird dadurch 
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Abb. 18. Widerstandsortskurve der rhombusförmigen Flächenantenne 


vergleichbar mit einem 90°-Kreissektor. Ferner ver- 
hält sie sich in Fußpunktsnähe ungefähr wie die ent- 


sprechende Exponentialantenne, besitzt aber im 
Gegensatz zu dieser eine kontinuierlich sich schließende 
Berandung. 


Ihre Widerstandsortskurve ist in Abb. 17 wieder- 
gegeben. Sie spultim gesamten betrachteten Frequenz- 
bereich enger auf als die der Kreissektor-Antenne mit 
einem Öffnungswinkel von 90°. Analoge Überlegungen 
lassen eine gegenüber der 90°-Kreissektor-Antenne 
bessere Breitbandigkeit auch bei der von Artur [6] 
angegebenen Rhombusantenne erwarten. Ihre Wider- 
standsortskurve ist hier in Abb. 18 aufgetragen. Der 
mittlere Wellenwiderstand ist kleiner als der des 
90°-Sektors. 

Abb. 19 zeigt schließlich Widerstandsortskurven 
dreieckförmiger ebener Flächenantennen. 


als gedämpfte oszillierende Funktion des Verhältn s 


ses ee meet, Es wurden dabei Wider 

Wellenlänge 
standsänderungen bis zu 20% beobachtet. Auf diese 
resonanzartige Verhalten des Grundbleches sind offen 
bar die starken Streuungen der Meßpunkte bei de) 
vorliegenden Widerstandsortskurven zurückzuführen 


2. Richtdiagramme 


ı) Anordnung zur Aufnahme von Richtdiagrammen 


Als Meßplatz konnten wir ein Segelfluggelände be- 
tzen. Sende- und Empfangsstation waren etwa 500m 


—N 


A 
->r 
x=0 
bb. 20a. Axialsymmetrisch Abb.20b. Halbe lemniskatenförmige 
erregte ebene Flächen- ebene Flächenantenne, axialsymmetrisch 
antennen erregt 
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Abb. 21. Richtdiagramme des Systems Grundplatie— 90° Sektor. 
xy-Ebene nach Abb. 1 


‘von den an einer Bahnlinie und einer Industrieanlage 
hervorgerufenen Streufeldern, wurde als Sende- 
"antenne eine Richtantenne aus 32 Lemniskatendipolen 
mit Reflektor aufgebaut. Die Symmetrisierung der 
‘von dem unsymmetrisch arbeitenden Sender geliefer- 


oneinander entfernt. Um möglichst frei zu werden. 
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Abb. 22. Richtdiagramme des Systems Grundplatte — 90° Sektor, 


xz-Ebene nach Abb. 1 


Abb. 23. Richtdiagramme des Systems Grundplatte — 90° Sektor. 
yz-Ebene nach Abb. 1 


ten Spannung geschah über eine kontinuierlich nach- 
stellbare A/2-Umwegleitung, die vom Verfasser als 
geschlitzte Koaxialleitung in der Art einer Meßleitung 
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Abb.24. Richtdiagramme des Systems Grundplatte— Halbe Lemniskate. 
xy-Ebene nach Abb. 20b 
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Abb. 25..Richtdiagramme des Systems Grundplatte— Halbe Lemniskate. 
xz-Ebene nach Abb. 20b 


angewandte Ph} 


entworfen wurde. Die durchzumessenden Antenne 
wurden als Empfangsantennen benutzt. Das Syste 
Grundplatte— Antenne befand sich auf dem Da« 
eines Meßstandes, der um seine Mittelachse drehb: 
war. Das quadratische Grundblech hatte eine Seiteı 
länge von 2m, die Flächenantennen eine Höhe ve 
30 cm. Die Empfangsspannung an der Antenne wur 
über eine Koaxialkabel auf ein Rohde & Schwar 
Empfangsgerät gegeben und dort abgelesen. 


b) Richtdiagramme ebener Flächenantennen 


Die so gemessenen Vertikaldiagramme des Sy 
stems Grundblech—90°-Sektor zeigen die Abb. 21 un 
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Abb.26. Richtdiagramme des Systems Grundplatte — Halbe Lemniskate, 
yz-Ebene nach Abb. 20b 


22 bei verschiedenen Werten von #r,. Die Vertikal- 
kennlinien des Systems 


Grundblech— Halbe Lemniskate 


sind in Abb. 24 und 25 wiedergegeben. In unseren 
theoretischen Betrachtungen haben wir eine unendlich 
große Grundebene vorausgesetzt. Das unter dieser 
Voraussetzung berechnete Fernfeld ist aber sehr ver 
schieden von dem Feld, das man bei Berücksichtigung 
einer endlichen Grundplatte erhält. Das bewiesen 
LEITNER und SPENcE [21] durch Berechnung des 
Feldes einer A/4-Stabantenne über einer runden, un. 
endlich dünnen leitenden Platte. Die unter Benutzung 
der Wellenfunktionen des abgeplatteten Rotations 
ellipsoids und unter der Annahme einer sin-förmigen 
Stromverteilung auf der Antenne abgeleitete Lösung 
zeigt folgendes: 

1. Richtdiagramme für verschiedene Verhältnisse 


Grundplattendurchmesser 
Wellenlänge 


ıterscheiden sich ziemlich stark voneinander und von 
sm Richtdiagramm einer A/4-Stabantenne über un- 
ıdlich ausgedehnter, sehr gut leitender Ebene. 

2. Während diese die Energie am stärksten in eine 
ichtung senkrecht zur Antenne abstrahlt, strahlt das 
'ystem, das aus einem endlichen Grundblech und einer 
/4-Stabantenne besteht, in eine andere Richtung. 
3. Infolge des auch auf der Rückseite des endlichen 
\rundbleches fließenden Stromes ist eine insbesondere 
ei kleinen Werten des Verhältnisses 


Grundplattendurchmesser 
Wellenlänge 


arke Rückstrahlung vorhanden. 

Die gemessenen Vertikalkennlinien in den Abb. 21, 

2, 24 und 25 bestätigen die beiden letzten Aussagen. 
Insere berechneten Vertikaldiagramme, die das von 
em endlichen Grundblech herrührende Streufeld 
\icht berücksichtigen, sind natürlich mit diesen ge- 
nessenen Diagrammen nicht vergleichbar. 
Ein Vergleich zwischen Rechnung und Messung 
lürfte am ehesten bei den Horizontaldiagrammen 
lurchzuführen sein. Allerdings wird auch hier das 
Streufeld der quadratischen Grundplatte eine Rolle 
;pielen. Da ferner die Antenne maximal in eine Rich- 
ung schräg zur Grundebene strahlt, befindet man sich 
n dieser Ebene auf den Flanken der Hauptkeule. Es 
müssen deshalb hohe Anforderungen an die Meßgenau- 
igkeit gestellt werden. Es wurde jeweils das gesamte 
Horizontaldiagramm gemessen. Da das Feld in allen 
vier Oktanten des betrachteten Halbraumes aus 
Symmetriegründen gleich groß sein muß, kann über 
die vier Meßwerte in einander entsprechenden Raum- 
richtungen gemittelt werden. Dann erhält man die 
‘wieder symmetrisch ergänzten Horizontaldiagramme 
in Abb. 23 für eine 90°-Kreissektor-Antenne. 

Bis r,/A—=0,732 gibt unsere Rechnung (s. Abb.13) 
‚qualitativ gut diese Messungen wieder. Oberhalb von 
r9A=0,732 konvergiert unsere Reihenentwicklung 
schlecht. Sie gibt uns bei r,/A=1,000 nur noch einen 
Hinweis darauf, daß die Seitenstrahlung der Antenne 
die Vorwärtsstrahlung jetzt weit überwiegt, was durch 
die Messung bestätigt wird. 

Die Horizontaldiagramme des Systems Grund- 
blech— Halbe Lemniskate in Abb. 26 fallen bei kleinen 
Werten von x:h mit den entsprechenden Richt- 
diagrammen der 90°-Kreissektor-Antenne zusammen 
und weisen bei großen Werten von x: h geringere 
Einschnürungen auf als diese. Auch hierin kann man 
eine Verbesserung der Breitbandigkeit gegenüber der 
'90°-Kreissektor-Antenne sehen. 


Zusammenfassung 


In Anlehnung an Arbeiten von SCHELKUNOFF und 
Paras u. Kına über rotationssymmetrische Kegel- 
antennen wird die Feldtheorie einer nichtrotations- 
symmetrischen. Antenne entwickelt. Als spezielle 
Lösung der Maxwellschen Gleichungen gewinnt man 
die Lecher-Welle in einer elliptischen Kegelleitung. 
Die Sektorleitung ist eine entartete elliptische Kegel- 
leitung. Ihr Wellenwiderstand wird angegeben. 

Eine allgemeine Lösung der Maxwellschen Glei- 
chungen führt zu den E-Wellen und H-Wellen der 
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Sektorleitung. Die radiale Komponente des elektri- 
schen bzw. magnetischen Strahlungsvektors wird nach 
Lameschen Produkten und sphärischen Zylinder- 
funktionen entwickelt. Das Strahlungsfeld einer. 
Kreissektor-Antenne wird durch Kugelwellen elek- 
trischen Typs dargestellt. Theoretischer und experi- 
menteller Befund stimmen bei einer Kreissektor- 
Antenne mit einem Öffnungswinkel von 90° befriedi- 
gend überein. Bei den Widerstands- und Richtdia- 
gramm-Messungen wurden unsymmetrische Antennen 
benutzt. 


Unter allen Kreissektor-Antennen besitzt die mit 
einem Öffnungswinkel von 90° optimales Breitband- 
verhalten. Leitungstheoretische Überlegungen zur 
Verbesserung der Breitbandigkeit führen zur lemnis- 
katenförmigen ebenen Flächenantenne. Die berech- 
nete Stromdichteverteilung auf Kreissektor-Antennen 
zeigt das erwartete Verhalten: Maximale Stromdichte 
an den Kanten des Sektors, minimale Stromdichte in 
der Mitte des Sektors. 


Diese Arbeit wurde in den Jahren 1956 bis 1959 im 
Institut für Angewandte Physik der Universität 
Marburg a.d. Lahn angefertigt. Ich möchte an dieser 
Stelle meinem hochverehrten Lehrer, Herrn Professor 
Dr. WorLrer, für die Überlassung des Themas, für 
wertvolle Anregungen und für die Förderung der 
Arbeit danken. Der Deutschen Forschungsgemein- 
schaft bin ich für die Bereitstellung umfangreicher 
experimenteller Hilfsmittel zu großem Dank ver- 
pflichtet. 
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Einleitung 


Das Problem der Beugung elektromagnetischer 
Wellen an Körpern beliebiger Form und Größe ist 
weitläufig und kompliziert sowohl auf analytischen 
wie auch auf meßtechnischen Wegen. Unter den ver- 
schiedenen Angriffsmethoden, die Hoffnung geben ein 
tieferes Verständnis der Beugung im allgemeinen zu 
gewinnen, ist besonders das eingehende Studium des 
einfachen Spezialfalls zu emp- 
fehlen. In diesem Sinn soll hier 
der Versuch gemacht werden, 
ein möglichst anschauliches 
und vollständiges Bild der Beu- 
gung und Streuung am leiten- 
den Kreiszylinder zu bringen. 
Entsprechendes läßt sich auch 
bei der Kugel durchführen. 


Das gebeugte Feld des leitenden 
Kreiszylinders mit E parallelzur 
Achse 


Die Beugung einer ebenen 
elektromagnetischen Welle am 
leitenden Kreiszylinder ist in 
vielen Lehrbüchern mathema- 
tisch behandelt. Das Problem läßt sich skalar lösen, 
einerseits mit Dirichletschen Grenzbedingungen, wenn 
der elektrische Vektor, andererseits mit Neumannschen 
Grenzbedingungen, wenn der magnetische Vektor zur 
Zylinderachse parallel steht. In beiden Fällen läßt sich 
die Lösung für das axiale Feld als Überlagerung des 
Streufeldes und der einfallenden ebenen Welle dar- 
stellen. Beispielsweise beim Dirichletschen Problem hat 
das Gesamtfeld folgende Form: 


E,=|E,|e%& =Ei+ BE} 


Abb.1. Ebene und Zylinder 
mit senkrecht und schräg ein- 
fallender Welle 


JIm(ka) (1) 


EIPFE Hn’(ka) 
Das einfallende elektrische Feld E! hat die Amplitude 
eins und ist parallel zur Zylinderachse gerichtet, die 
mit der zAchse des zylindrischen Koordinaten- 
systems 0, ©, z zusammenfällt. Die Welle schreitet in 
der Richtung der positiven x Achse bei O—=0 fort. 
Die entsprechende Formel für 4, beim Neumannschen 
Problem ist wie (1) aber mit den Ableitungen von 
J, (ka) und H/ (ka). 

Was besagt nun diese Formel? Physikalisch über- 
sichtlich ist sie sicherlich nicht. Man nimmt aber 
natürlich an, daß ein sinnvolles Bild des gebeugten 
elektrischen Feldes aus ihr gewonnen werden kann, 
wenn man sich die Zeit nimmt, sie numerisch auszu- 
werten. Ist das aber praktisch möglich ? Bei Zylindern 
mit kleinen Werten von ka, ist die Antwort ja; bei 
Zylindern mit großen ka hingegen, wäre ein fantastischer 
Rechenaufwand nötig, wegen der schlechten Konver- 
genz der Reihe, die das radial fortschreitende Streu- 
feld darstellt. Mit anderen Worten, die exakte Lösung 


— eikocos® _ 


H® (k 0) eimO 


des einfachsten Beugungsproblems ist praktisch v 
wenig Wert bei großen Zylindern! 
Als Einleitung zum Studium des Gesamtfeldes d( 
leitenden Kreiszylinders, vorerst bei der Dirichletsche@ 
Grenzbedingung, ist es lehrreich, die Feldverteilung 
der Mitte der beleuchteten Hälfte zu betrachten. H 
ist die Streuung eigentlich nur eine Spiegelung, beifl; 
der die Mantelfläche des Zylinders als eine Zusammen 
setzung schmaler ebener Elemente angesehen werde 
kann, die von der einfallenden Welle bei verschiedene 
Winkeln beleuchtet werden. Diese naive Darstellung 
ist in Abb. 1 gezeigt, oben für senkrechten Einfal 
unten für schrägen Einfall. In beiden Fällen bild 
die Überlagerung der einfallenden Welle, Zi, und des 
gespiegelten Welle eine stehende Welle, £,, ind 
Richtung der negativen x Achse und senkrecht z 
spiegelnden Fläche. Bei senkrechtem Einfall ist 
Ei=cos(vt—kx), E,=2sinkxsinot; die Wellene 
länge, A,, der stehenden Welle ist gleich der Welle 
länge, A,, des einfallenden Feldes; die Phasengeschwir 
digkeit, v, ‚längs der Spiegelfläche ist unendlich großS 
Bei schrägem Einfall ist Zi =cos [ot— k (= cos W- 
ysin Y), B,=2sin(kx cos Y) sin(ot— kysin Y); di 
Wellealirge der stehenden Welle ist A, —4,/|cos X 
die Phasengeschwindigkeit ist v,,—vj/|sin Y|, wo ı 
die Lichtgeschwindigkeit ist. Die stehende Welle be 
senkrechtem Einfall auf einen ebenen Spiegel in de 
Ebene kxe=—3,1 ist in Abb. 2 graphisch dargestelli 
Die einfallende Welle läuft von oben nach unten. Di 
sinusförmige Amplitudenverteilung mit Maximal 
wert 2, und die im Zickzack verlaufende Phasen 
verteilung sind bzw. gestrichelt und punktiert ge 
zeichnet. 
Auf derselben Abbildung sind auch die entspre 
chenden exakten Kurven des elektrischen Gesamt 
feldes eines leitenden Kreiszylinders mit ka=3,1, fü 
die Richtung © =180° der einfallenden ebenen Welle 
gezeigt. Die nahe Verwandtschaft der entsprechender 
Kurvenpaare bei ka=3,1 und ka= © ist sofort er 
sichtlich.. Insbesondere ist die Entfernung zwische 
benachbarten Maximalwerten dieselbe halbe Wellen- 
länge bei der Rückstreuung am Zylinder wie bei deı 
Spiegelung an der Ebene. Es ist für spätere Zwecke” 
bemerkenswert, daß die Amplituden und Phasen für 
alle Zylinder mit ka Werten größer als 3,1 zwischen 
die entsprechenden Kurven für ka=3,1 und ka = 
fallen müssen, so daß ihr allgemeiner Verlauf ohne 
weiteres als bekannt angesehen werden darf. 
Es ist nun an der Reihe, das Gesamtfeld des lei- 
tenden Zylinders im ebenen einfallenden Feld zu be 
trachten. Zu diesem Zweck ist der Wert ka—=3,1 zu 
empfehlen. Er ist klein genug, daß eine numerische 
Auswertung des Feldes mit 6 bis 9 Termen der Reihe (1) 
möglich ist, aber auch groß genug, daß der Kreis 
umfang drei Wellenlängen beträgt. Die Reihe wurde 
bei jeweils konstantem x für zunehmende Werte von % 
mit Handrechenmaschine summiert. Das errechnete 
Feld ist ausführlich in Abb. 3 dargestellt. Folgende 
Einzelheiten sind zu bemerken. 


1. Die Ebene des Bildes ist die xy Ebene 
»s rechtwinkligen Koordinatensystems 
der positiven x Achse nach unten und 
positiven y Achse nach rechts. Die 
des Querschnitts des beugenden 
inders mit Halbmesser @ ist schraffiert 
eigt; die Zylinderachse fällt mit der 
(chse zusammen. 
'2. Das einfallende Feld ist die ebene 
lle eik2, die in der Richtung der posi- 
an x Achse fortschreitet. Ihre Amplitude 
überall gleich eins; ihre Phase, O'—kx, 
"auf die y Achse bezogen, wo sie den 
rt null hat. Ihr Verlauf ist auf der Ab- 
lung durch den von oben links nach 
en rechts laufenden schrägen Strich 
ch die Mitte des Zylinderquerschnitts 
‚gestellt. 

3. Die Amplitude des Feldes längs der 
sativen x Achse, d.h., E, bei ky=0, die 
ıon im vorhergehenden Bild gezeigt 

de, ist oben in Abb. 3 eingezeichnet, 

d zwar so, daß der Wert eins mit der 
sativen x Achse zusammenfällt. Die 
ase, O5, bei ky=0 ist gestrichelt oben 
'ks gezeigt. 

4. Die Amplitude, E,, des Ge- 
mtfeldes ist auf entsprechende 
eise in ausgezogenen Strichen ent- 
ag der Geraden kx=0, 4, und 
nr gezeigt. An der Geraden, kx =0, 
ngt die Amplitude mit null am 
lindermantel an, und schwingt 
ınn mit wachsendem ky um den 
ert eins der einfallenden Welle. 

Verlauf ist ähnlich wie der schon 
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Abb.2. Amplitude und Phase des reflektierten Feldes auf der beleuchteten Seite eines 
Zylinders und einer Ebene bei senkrechtem Einfall. E:=|E;| e"@E entlang der negativen 
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schriebene entlang der negativen 
Achse, aber die Länge der stehen- 
»n Wellen ist bedeutend größer. 
it wachsendem kx nimmt diese 
ellenlänge auch beständig zu, wie 
ıs den Amplitudenverteilungen bei 
?—4r und kx—= 207 ersichtlich ist. 


5. Der Phasenunterschied, 
- 9, 

rischen dem Gesamtfeld und dem 
nfallenden Feld, ist für dieselben 
Terte von k.x gestrichelt eingezeich- 
et. Zu bemerken ist, daß außer in 
r Nähe des Zylinders, dieMaximal- 
nd Minimalwerte der Amplitude 
it den Nullwerten des Phasen- 
nterschieds zusammenfallen. Ent- 
prechenderweise sind die Maximal- 
nd Minimalwerte des Phasenunter- 
chieds an denselben Stellen wie die 
erte |E,|=|E2| =1, wo die Ge- 
mtamplitude gleich der Amplitude 
er einfallenden Welle ist. 

6. Es ist nun leicht, zwei para- 
olisch geformte Kurvenscharen ein- 
uzeichnen, die jweils Punkte ma- 
“maler oder minimaler Amplitude 
ereinigen (punktiert auf Abb. 3), 
der Punkte verbinden an welchen 
lie Gesamtamplitude beständig den 
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ungestörten Wert |Z,| =1 der einfallenden Welle hat 
(strichpunktiert auf Abb. 3). In der Nähe des Zylin- 
ders, aber nur auf seiner Schattenseite, trennt sich die 
punktierte Kurve minimaler Amplitude |Z; min. von 
der Kurvemit dem Phasenunterschied null, O9,— 0'=0, 
nahe bei ke—=4n, ky=3. Erstere erreicht den Zy- 
lindermantel bei 9 =0°, letztere (gestrichelt) bei 
9=50). 

Nach diesen Vorbemerkungen läßt sich das Ge- 
samtfeld um den beugenden Zylinder sinnvoll er- 
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Abb. 4. Gemessene Amplituden- und Phasenverteilung um einen leitenden Zylinder 
mit ka=3,1 bei ebener, einfallender Welle mit dem elektrischen Vektor parallel 


zur Zylinderachse. Messungen von KODIs, PICKARD und ALTSHULER 


klären. Von der negativen x Achse bewegen sich 
fortschreitende Wellen nach beiden Seiten parallel zu 
den parabolischen Kurvenscharen. Auf ihnen reiten 
die zu ihnen transversal verlaufenden stehenden 
Wellen, und zwar so, daß die strichpunktierten para- 
bolischen Kurven (Abb. 3) immer die unveränderte 
Amplitude der einfallenien Welle darstellen, während 
die punktierten Linien’ zwischen ihnen immer die 
maximale oder minimale Amplitude der stehenden 
Welle bezeichnen. In der Nähe der negativen x Achse 
ist die Phasengeschwindigkeit der forteilenden Wellen 
größer als die Lichtgeschwindigkeit, aber sie nimmt 
beständig ab, so daß sie sich mit wachsender Entfer- 
nung der Lichtgeschwindigskeit nähert. 

Die längs der parabolischen Kurven fortschrei- 
tenden Wellen folgen mit ihren Rändern und der 
Amplitude null dem Zylindermantel. Die Überlage- 
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rung der von beiden Richtungen in den geometrise] 
Schatten eindringenden Wellen kleiner Ampliti 
bildet stehende Wellen in transversaler Richtung, 
aus Symmetriegründen ein kleines, relatives Maxim 
längs der positiven x Achse bilden. Die Amplit 
dieses Maximums, d.h., |Z,| bei ky—=0, ist in 
linken Häfte von Abb. 3 gezeigt. Sie fängt mit 
Wert null am Zylindermantel an, und steigt 1 
wachsendem kx gegen den Wert eins empor. Die @ 
sprechende Phase fängt am Zylinder mit einem W 
von 342° an, und nähert sich asymptotisch di 
Wert, ©°, der einfallenden Welle. Sie ist 
strichelt auf Abb. 3 gezeigt, nahe bei der ( 
raden © bei ky=0. 

Das elektrische Gesamtfeld um einen leite® 
den Kreiszylinder ist von Kopıs [l] und vo 
PICKARD und ALTSHULER gemessen worden. Ih® 
Resultate sind auf Abb. 4 zusammengefaßt. Wi 
ein schneller Blick andeutet und eine gena 
Untersuchung bestätigt, stimmen sie in al 
Einzelheiten sehr gut mit der Theorie, was 
diesem Fall als Rechtfertigung des Meßverfak, 
rens und der mit ihm verbundenen Korrekture 
gedeutet werden darf. Die für die Messung 
bevorzugte Apparatur besteht aus zwei parallele 
Metallspiegeln mit absorbierender Umrand ; 
zwischen denen das strahlende Horn, der be 
gende Zylinder, und die bewegbare Sonde ang 
bracht sind. 


Das gebeugte Feld des leitenden Kreiszylinders \. 
mit H parallel zur Achse _ 1 


Wie schon früher angedeutet wurde, ist d 
Gesamtfeld, , =H! +.H}, wenn das magnetise 
Feld zur Zylinderachse parallel steht, durch eim® 
Formel ähnlich wie (1) gegeben. Errechnete Wer 
für einen Zylinder mit ka=3,1 sind auf Abb 
dargestellt. Der allgemeine Eindruck ist wie de 
der Abb. 3, aber es bestehen wichtige Untem 
schiede, besonders in der Nähe des Zylinders‘ 
In Abb.3 verschwindet das elektrische Fel@ 
überall auf der Mantelfläche des Zylinders; 
Abb. 5 hingegen, verschwindet die Ableitung def 
magnetischen Feldes in Richtung der Normale 
so daß H, auf der Mantelfläche ein Maximu 
oder ein Minimum haben muß. Wie in Abb. 5 
sehen ist, erreicht die radiale stehendeWelle b 
© =180° (also entlang der negativen x Achs 
ein Maximum auf der Oberfläche des Zylinders. Dies 
Maximum bewegt sich in beiden Richtungen um dem 
Zylinder. Es ist die Hüllfläche zweier fortschreitend! 
Wellen, die in entgegengesetzten Richtungen um dem 
Zylinder eilen. Dabei breitet sich die radiale stehend@ 
Welle aus, bis bei 9 =108° der Wellengipfel am Zylüi 
dermantel abgeflacht wird um ein kleines, relativ 
Minimum zu bilden. Auf diese Weise verliert das 
Wellenmaximum Kontakt mit dem Zylinder un 
wandert von ihm fort. Die stehende Welle bei 9 = 
(also entlang der positiven y Achse) zeigt in Abb.& 
eine kleine Vertiefung im Amplitudengipfel am Zylin 
dermantel, eben weil das Maximum, welches ve 
© =180° der Zylinderfläche folgte, sich bei O = 10 
in ein Minimum verwandelte. Dieses Minimum bleil 
nun in Berührung mit dem Zylinder wie man mi 
abnehmendem © in den geometrischen Schatten e 
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sehen. Die Mitte jeder der Abbildung ist die Mitte der 
Schattenseite. In Abb. 7 ist der elektrische Vektor 
parallel zur Zylinderachse und der axiale Strom, 
K,(9), ist in normalisierter Form dargestellt. Ind 
Abb. 8 ist der magnetische Vektor parallel zu Zylinder- 
achse, und der Querstrom, Ko(®), ist gezeigt. Die 
von Kopıs [1] gemessenen Amplituden und Phasen 
sind auch in Abb. 7 angegeben. 


ot. Dabei breitet es sich aus, und der Amplituden- 
el, in den es eingepreßt ist, senkt sich bis bei 
‚22° das relative Minimum auf einem Gipfel sich 
Sin relatives Maximum in einem Amplitudental 
vandelt. Dieses Maximum wächst in einem Tal, 
sich weiter senkt bis 9 =0°, während das Mini- 
n bei 9=22° die Berührung mit dem Zylinder- 
tel verliert und sich in radialer Richtung von ihm 
ernt. Das Maximum in |H,| bei en 


0° am Zylindermantel ist der _— A En 2 ee A Pl rl a 
ang einer Hüllfläche maximaler ma 120° -30° 0. 30° 720° ma. mo mm° 210° 2520° 
plitude, die längs der positiven ENT A KA 
hse besteht und mit wachsender rw) re 7 7 s_ 0 W049, 
fernung vom Zylinder dem Werte Qyfürky-0 N | | Leitender Zylinder 
, zustrebt. Diese Amplitude und | \/hffüryo | Aa=31 | 
‚dazu gehörende Phase sind auch WT0° | er | 
Abb. 5 gezeigt. we S Mmjfüke | |, 
Analog dem früher besprochenen | SL ; | [el \ao 
des elektrischen Feldes bei 704-0° Saul ag Rz 
Dirichletschen Grenzbedingung x | m 
, | | 
E> verlaufen die Konturen ", DES A, für ko ft | N, op 
em Amplitudenwert, |Z|= " „1 \ 9,firke | 
—], im allgemeinen parabo- 4r ni A 
h und fallen mit den Linien, u KUN %, kler ern Ir 
‘ denen Op — ©‘ ein Maximum „| u N ara 7 u u 
sr ein Minimum hat, zusammen. \ ANNE gs me 
1e Ausnahme ist die engere Um- pen: ---— 8 = IH 
hung des Zylinders. Mit derselben NSS en 
4 INN Ha] für ka=uc 
02 


schränkung sind auch die Kon- 

en, auf denen O4 = ©' para- 

lisch, und diese fallen mit den kr 
aien, auf denen |A,| ein Maximum 
er Minimum hat, zusammen. In 
Nähe des Zylinders ist das Ver- 
ten aller dieser Konturen kom- 
'ziert, wie in Abb. 6 zu sehen ist. 
»r Grund dafür ist in den stehen- 
n Wellen des Flächenstroms auf m) 


U) 


m Zylindermantel zu suchen. Diese ZUIE ga 
ben eine viel größere Amplitude, | | le 
nn H der Zylinderachse parallel 2 | } Ark ham 1-20 
sht, als wenn E zur Achse par- 9380 Be? 
el ist. 2 Ba 
5 74 60° 
Der Flächenstrom EN 1077 Mr 
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Hr=/H, 10H 
N 5 = Male & 30° | 70° | weo® muo® \n® | zuo® 2520° 
Die Quellen des Streufeldes sind a u N 
Br en d u VER 
e Flächenströme auf dem Mantel wo, 4 
3 beugenden Zylinders. Sind sie Abb. 5. Amplituden- und Phasenverteilung um einen leitenden Zylinder mit ka=3,1 bei ebener 
kannt, so läßt sich das von dem einfallender Welle mit magnetischem Vektor parallel zur Zylinderachse 


vlinder ausgestrahlte Sekundärfeld 


'enigstens prinzipiell errechnen. Umgekehrt, ist das 
'eld bekannt, so kann der Flächenstrom K (9) direkt 
on der Tangentialkomponente des magnetischen Fel- 
»s bestimmt werden. Bei der Dirichletschen Grenz- 
Sdingung gilt am Zylindermantel, E,—=0, — oE,|0o = 
»Bo—=iwu,K,(0). Es läßt sich auch nachweisen, daß 
er Phasenwinkel 9, des Flächenstroms, K,(9), sich 
m 90° von dem Phasenwinkel, O7, des elektrischen 
eldes unterscheidet, so daß 9x =9, +9’, E =. Bei 
er Neumannschen Grenzbedingung ist der Flächen- 
rom |K6(9)| =|R.|, o = «a. 

In den Abb. 7 und 8 sind die von den Reihenent- 
icklungen errechneten Kurven für die Amplituden- 
erteilung und die Phasenverteilung des Flächenstroms. 
beiden Hälften des Zylinders mit ka=3,1 zu 


Wenn die Phase, ©, des Flächenstroms, K,(0), in 
Abb. 7 in ihrer Abhängigkeit vom Winkel, ©, be- 
trachtet wird, so stellt man leicht fest, daß von der 
Schattengrenze bei O=90° bis fast zur Mitte des 
Schattengebiets sie sich nahezu linear mit © ver- 
ändert. Die Phasengeschwindigkeit im Schatten ist 
also konstant, und zwar beträgt sie etwa 0,7 der Licht- 
geschwindigkeit. Offenbar läßt sich der Flächenstrom 
als zwei fortschreitende Wellen darstellen, die von 
beiden Schattengrenzen mit konstanter Phasen- 
geschwindigkeit und schnell abnehmender Amplitude 
um den Zylinder laufen. Die Interferenz zwischen 
diesen zwei Wellenzügen, die sich in entgegengesetzter 
Richtung bewegen, dient zur Erklärung der stehenden 
Wellen mitten im Schattengebiet. Aus der Phasen- 
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verteilung im beleuchteten Gebiet folgt, daß hier die 
Phasengeschwindigkeit von einem fast unendlichen 
Wert bei @ =180° beständig abnimmt und die Licht- 
geschwindigkeit unweit von O=90° erreicht. 

Der Flächenstrom, Ke(0), bei der anderen Pola- 
risation, der in Abb. 8 dargestellt ist, kann auf ent- 
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Abb. 6. Konturen maximaler und minimaler Amplitude und Phasendiffe- 
renz, Konturen mit Amplitude oder Phase gleich der der einfallenden Welle 
in der Umgebung des Zylinders in Abb. 5 
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Abb. 7. Flächenstromverteilung auf leitendem Zylinder mit ka=3,1 bei 


Fu einfallender Welle mit dem elektrischen Vektor parallel zur 
Zylinderachse 


sprechende Weise beschrieben werden. In diesem Fall 
nimmt aber die Amplitude der zwei fortschreitenden 
Welien viel langsamer ab, so daß die stehenden Wellen 
im Schattengebiet viel ausgebreiteter sind und be- 
deutend größere Amplituden haben. Es sind diese 
stehenden Wellen, die den komplizierten Verlauf der 
Feldverteilung in der Nähe des Zylinders verursachen. 

Die Bestimmung des Flächenstroms auf dem lei- 
tenden Zylinder mit ka—=3,1 ist interessant aber un- 
nötig, da das Gesamtfeld schon bekannt ist. Bei sehr 


großen Werten von ka hingegen, ist die lang 1 
konvergierende Reihe zur Bestimmung des Fe 

nicht mehr nützlich, so daß andere Wege zu sei 
Ermittlung gesucht werden müssen. Eine Möglichl 
ist die Bestimmung des Flächenstroms, und dann 
Errechnung des von ihm erzeugten Streufeldes. 
Flächenstrom läßt sich aber natürlich auch nicht 
der Reihenentwicklung nach Eigenfunktionen 
rechnen, wiederum wegen schlechter Konverge 
Um dieser Schwierigkeit abzuhelfen, hat WETZEL 
eine asymptotische Formel für nicht zu kleine W: 
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Abb.8. Wie Abb. 7 aber mit magnetischem Vektor parallel 
zur Zylinderachse 


von ka abgeleitet. Für K,(©) hat sie folgende For 


2 


m=0 


ivrB+iv2am 


U“ 


H,% (ka) 


+ f: a) 


Als neue Variabel it ®=3n-0© eingeführt. 
merkenswert ist, daß die in (2) symmetrisch auftrete 
den unabhängigen Variabeln, die von beiden Schatte 
grenzen gemessenen Winkel ® sind. In der Form | 
läßt sich der Flächenstrom leider auch nicht leie 

auswerten. Aber folgende Näherungsformeln sind 

diesem Zweck nützlich. 

Im mittleren Teil des beleuchtenden Gebiets 
scheidet sich die asymptotische Formel (2) nur wen 
von der einfachen geometrisch-optischen Stromvert 
lung, 


An, 
kası 


K,(d) dv 


er (a D)+ivr2rm 


Hole) (ka) 


K,(0) 


m — 0050 eikacosO | { 
Im ganzen Schattengebiet, an den Schattengrei 


zen und noch ein gutes Stück in das beleuchtete Gebi 
hinein läßt sich (2) auf folgende Näherung bringe 
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+ TEN era: Aa 


Su 
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} (D+2nm), 


„(ef U(n—®) +27m], 


7 = Bed (o) 
Tı 
w(r) = [er HP dt. (7) 
Tr 


rzers Formel (4) läßt sich offenbar als eine Über- 
rung zweier fortschreitender Wellen deuten, die 
ntgegengesetzten Richtungen um den Zylinder 
1. (In Wirklichkeit sind nur die ersten Terme mit 
-0 von praktischer Bedeutung.) Die Phasen dieser 
len hängen linear von ® bzw. von n— ® ab; die 
plituden sind durch die Funktion 9 (£) gegeben. In 
ter Näherung hat diese Funktion die Form (6). 


Abb. 9. WETZELs Funktionen 9(£) und g(&) 


einer besseren, ersten Näherung, wird 9(£) durch 
- Funktion 7,(£) ersetzt, die folgendermaßen defi- 
rt ist: 3 

He = +.) (8) 


e Amplitudenfunktion 9(£) und ihre erste Korrek- 
nsfunktion ©(£) sind von WErzEL [2] ausgewertet 
rden. Erstere ist in Abb. 9 dargestellt. 

Der Formel (4) für den Flächenstrom bei der 
Polarisation entspricht folgender Ausdruck bei der 
Polarisation: 


Ke(O)n >. (&,) eika(0+22m) L 


MmM= 


(9) 
+ g(&),) eitella-+2=m]), 


) die neue Amplitudenfunktion g(£,,) wie in (6) defi- 
ert ist, aber mit w’(r) an Stelle von w(T); 


w(r) = [ ert"letdt. (10) 

T; 
ne Funktion erster Ordnung läßt sich nach dem 
uster von (8) einführen. Die Funktion g(£) ist auch 

Abb. 9 dargestellt. Sie ist von WerzeL [2] ausge- 
srtet worden zusammen mit der Korrektionsfunk- 
In c(E). : 

Man ersieht aus Abb.9, daß die Amplituden 
(E)| und |g(&)| so schnell mit wachsendem ® ab- 
»hmen, daß nur die Terme mit m =0 einen bemerk- 
ıren Beitrag bringen können. Die Funktion |s(&)] 
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nimmt aber weit weniger schnell ab als |g (&)| ; sie ist 
in Abb. 8 gestrichelt eingezeichnet mit ka=3,1. Die 
Tatsache, daß für positives ® (also im Schatten) die 


Winkel Ö, (£) und ©, (£) linear mit ® zunehmen, deutet . 


auf eine konstante Phasengeschwindigkeit der um den 
Zylinder eilenden Wellen. Sie ergibt sich als kleiner 
als die Lichtgeschwindigkeit. 

Die Funktionen nullter Ordnung, 9(E£) und g(£), 
sind mit den bekannten von Fock [3] mittels der 
Parabel eingeführten Funktionen nahe verwandt. Sie 
hängen aber von der Bogenlänge entlang des Zylinder- 
mantels anstatt wie bei Fock von der Entfernung 
längs der Tangente ab. 

Um die Genauigkeit der neuen asymptotischen 
Formeln zu prüfen, wählte WerzeL den Wert ka = 12. 


Dieser ist groß genug, um allen Annahmen der asym- 
ptotischen Darstellung gerecht zu werden, ist aber auch 
@ Grad 
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Abb. 10. Flächenstrom auf einem leitenden Zylinder mit H parallel zur Achse 


und auf’einer Kugel in der Ebene senkrecht zu H; ka=12 


klein genug, um eine allerdings etwas weitläufige Aus- 
wertung der exakten, langsam konvergierenden Rei- 
henentwicklung zu gestatten. In Abb. 10 sind die 
nach beiden Formeln errechneten Amplituden und 
Phasen des Flächenstroms, K,(©), mit dem magneti- 
schen Feld parallel zur Zylinderachse gezeigt. Die 
ausgezogenen Kurven sind mittels Werzers Formel 
erster Ordnung (also mit 9,(8)) errechnet worden, die 
kleinen Kreise sind die von RIBLET angegebenen 
exakten Werte. Die Übereinstimmung im Gültig- 
keitsbereich, O<150°, der Wetzelschen Formel ist 
ausgezeichnet. Außerhalb dieses Bereichs, also in der 
Mitte des beleuchteten Gebiets, genügen der geome- 
trisch-optische Strom und eine leicht gemachte 
Kurvenverlängerung. (Auf Abb. 10 ist auch der von 
WerzeL [2] auf ähnliche Weise bestimmte Flächen- 
strom erster Ordnung für die entsprechende Polari- 
sation und Kreisumfang der leitenden Kugel gezeigt.) 

In Abb. 11 sind die ganz ähnlichen Resultate für 
den Kreiszylinder mit dem elektrischen Feld parallel 
zur Zylinderachse dargestellt. Die ausgezogenen 
Kurven sind wieder die aus W»rzers Formel erster 
Ordnung gewonnenen Amplituden; die kleinen Kreise 
sind von RIBLET aus der exakten Theorie errechnet. 
Die geometrisch-optische Stromverteilung ist ge- 
strichelt gezeigt. Zum Vergleich ist die früher be- 
sprochene Stromverteilung auf dem viel kleineren 
Zylinder mit ka—=3,1 strichpunktiert eingezeichnet. 
Ihr Verlauf ist sehr ähnlich wie der mit ka = 12, außer 
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im tiefen Schatten, wo es bei ka =12 mehr stehende 

Wellen gibt aber mit noch kleinerer Amplitude. Die 

entsprechenden Kurven für die Phasenverteilung sind 

in Abb. 12 zu sehen. Die Übereinstimmung ist nicht 

weniger gut. 
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Abb.11. Amplitudenverteilung des Flächenstroms auf einem leitenden 
Zylinder mit E parallel zur Achse und auf einer Kugel in. der Ebene 
senkrecht zu E 


Zusammenfassend darf die Behauptung nun wohl 
gemacht werden, daß das Problem der Stromvertei- 
lung auf dem leitenden Kreiszylinder in einer ebenen 
einfallenden Welle gelöst ist. Bei ka Werten kleiner 


g kleinen Amplitude können diese Unterschied: 
m__1ı _ WM m mM _M 60 m 2 s. nu: einen lokalen Effekt auf die Feldverteilung 
S zu : Fri zo haben. : 
S = 8 =ka 00s0=6876°058 ® za Ss Ohne auf weitere Einzelheiten einzugehen 
r— — 7) ıS . . 
Sa: @u(oder 6;), erste Näherung (Metze)) 3 Be S kann man schließen, daß der allgemeine Ver 
EI en En ER MINDER RT, = 1” & auf des Gesamtfeldes bei ka=12, oder eine 
RS ü eo. +, (oderd) exakte Berechnung een = Tu, = R =: 
S 71: , oz S beliebig größeren Wert, sich genau im Schem 
Sau für Zylinder (Riblet) _£ u ! 3 aueh ı 
IN N (9-05 = S des Feldes bei ka=3,1 einreihen muß. Dies 
za ka=3] -------- (9-9) S Schema ist in Abb. 3 dargestellt. Wie sch 
N re & früher erwähnt, ist die Amplituden- und Phase 
N wi & verteilung längs der negativen x Achse für al 
3 490 Beh S .ka-Werte zwischen 3,1 und co fast diesel 
R DERRT $ (Abb. 2), so daß die Anfänge der parabolisch 
Sr ern 798 Kurvenscharen, welche die stehenden Well 
oa 0 0 0 m mM m #0 MM des gebeugten Feldes kennzeichnen, ohne we 
Grad vom Zentrum der bestrahlten Seite (180-0) 5 2 ee 
; £ teres bekannt sind. Je nach der Größe d 
Abb. 12. Phasenverteilung des Fächenstroms auf einem Zylinder mit E parallel 


zur Achse und auf einer Kugel in der Ebene senkrecht zu E 


als etwa sechs ist die Reihenentwicklung nach Eigen- 
funktionen nützlich, bei größeren Werten von ka sind 
Werzers Formeln erster Ordnung sehr gute Annähe- 
rungen. (Dasselbe kann auch von der leitenden Kugel 
behauptet werden.) 


Das gebeugte Feld bei großem ka 


Mit dem Strombelag auf dem Zylindermantel be- 
kannt, läßt sich prinzipiell das elektromagnetische 
Streufeld bestimmen. In der Praxis ist aber die Aus- 
wertung der in Betracht kommenden Integrale schwie- 
rig. Deswegen ist es von Interesse, die Möglichkeit zu 
untersuchen einiges über die Beugung an Kreis- 
zylindern beliebiger Größe aus dem schon ausführlich 
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beschriebenen Bild der Beugung am Kreiszylinder 
ka=3,l zu gewinnen. 

Das Streufeld ist vollständig durch den auf d 
Zylindermantel induzierten Flächenstrom bestimi 
Jeder Unterschied zwischen den Feldern zweier be 
gender Zylinder mit verschiedenen Werten von ka 
durch einen entsprechenden Unterschied in den Strom 
verteilungen bedingt. Wie groß ist nun der Unt 
schied zwischen den bekannten Flächenströmen 
Zylindern mit ka=3,1 und ka=12? Es ist sd 
festgestellt worden, daß die Amplitudenverteilungf 
dieselbe ist, außer tief im Schatten, wo die Amplitu@, 
auf alle Fälle klein ist. In Abb. 12 kann man & 
Phasenverteilung dieser Ströme vergleichen. 
gende Kurven sind dargestellt. Die ausgezogen 
Kurve ist mittels der Wetzelschen Formel erster Ord 
nung errechnet, und mit der optischen Stromverteilum 
gestrichelt verlängert. Die kleinen Kreise sind d 
exakten Werte. Die mit langen unterbrochenen Stz 
chen gezeichnete Kurve stellt die Phase der einfalle 
den Welle auf dem Zylindermantel dar. Die strich 
punktierte Kurvezeigt den Phasenunterschied Oz 
zwischen dem elektrischen Gesamtfeld — im Lime 
wie man sich dem Zylindermantel nähert — und de 
einfallenden Feld bei ka=12. Die punktierte Kur 
zeigt denselben Phasenunterschied bei dem viel klei 
neren Zylinder mit ka =3,1. Man sieht, daß im ganz 
beleuchteten Gebiet und in einem großen Teil de 
Schattengebiets die Phasendifferenz fast unabhängi 
vom ka Wert ist. Sogar tief im Schatten, wo d 
Amplituden klein sind, ist die Phasenverteilung äht 
lich. Bei größerem ka gibt es einige mehr stehen 
Wellen und die Phasenverschiebung ist größer, we 
der Zylinderumfang größer ist. Aber eben wegen de 


Zylinders lassen sich dann die Kurvenschare 
selbst in guter Annäherung nach dem Schem 
bei ka=3,1 bestimmen. Das anschauliche Bild 
fortschreitenden Wellen parallel zu diesen Kurve 
scharen, wie auch die zu ihnen senkrechten stehende 
Wellen, bleiben unverändert. Auchim Schattengebiet 
sind nur wenige Unterschiede. Wie bei ka=3,li 
eine stark verminderte Amplitude hinter dem Zylind 
mit ka—=12 zu erwarten, und in diesem Amplitudent 
existieren stehende Wellen geringer Amplitude un 
mit einer Verteilung und Ausdehnung, die d 
stehenden Stromwellen auf dem Zylindermantel i 
der Mitte des Schattengebiets entsprechen. 

Die obige Diskussion bezieht sich auf die Z-Pol: 
risation. Bei der H-Polarisation läßt sich dasselb 
Argument durchführen, aber der lokale Effekt de 
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nden Stromwellen auf dem Zylinder ist ausge- 
ber. 

s ist nun klar, daß der allgemeine Verlauf des am 
‚den Kreiszylinder gebeugten Feldes bei beliebig 
'n Werten von ka angenähert bestimmt werden 
‚Mittels der bekannten Feldverteilung beika —3,1 
der aus der Wetzelschen Formel errechneten 
nverteilung für den in Betracht kommenden Wert 
ka läßt sich das Feld in nicht allzu großer Nähe 
Zylinders in allen wichtigen Einzelheiten be- 
iben. 

Der Streuquerschnitt 


n Stelle der gesamten Feldverteilung um einen 
enden Körper ist es oft genügend, einen räum- 
n Mittelwert der wiederausgestrahlten Energie zu 
mmen. Ein Mittelwert dieser Art ist der Streu- 
schnitt. Für den leitenden Kreiszylinder läßt 
aber der Streuquerschnitt pro Längeneinheit auch 
für kleine ka aus der Reihenentwicklung für das 
bestimmen. Für größere ka ist es Wu [4] ge- 
en, eine Reihe nach Potenzen von ka abzuleiten. 
ersten vier Glieder (Wu hat deren sechs errechnet) 
lie zwei Polarisationen sind im folgenden ange- 
n: 
"ür die Dirichletsche Grenzbedingung mit elek- 
hem Feld parallel zur Zylinderachse, 
2 =1 40,498 (ka)-3 — 0,0112 (ka) 3 a 
— 0,0147 (ka) ?. 
die Neumannsche Grenzbedingung mit magneti- 
m Feld parallel zur Zylinderachse, 


o 


4, =1-0,432(ka)-3— 0,214(ka) + | (19) 
4-0,0557 (ka)=2. 


(Entsprechende Formeln für die Kugel sind auch von 
Wu [4] und später auf anderem Weg von BECKMANN 
und Franz [5] bestimmt worden.) Der Vergleich 
dieser Formeln mit der exakten Theorie zeigt, daß nur 
das erste Korrektionsglied in (11) nötig ist, um eine 
ausgezeichnete Übereinstimmung zu versichern, wenn 
ka so klein ist wie eins. In (12) sind zwei Korrektions- 
glieder nötig, um eine gute Annäherung zu bekommen 
für Werte von: ka so klein wie vier. (Bei der Kugel 
ist die Übereinstimmung nicht weniger gut.) 

Offenbar kann der Streuquerschnitt des leitenden 
Kreiszylinders (und der leitenden Kugel) für alle Werte 
von ka leicht errechnet werden. Für ka klein genügt 
die exakte Lösung; für größere ka sind die Wuschen- 
Reihen zu verwenden. 


Schluß 


In dem hier gebrachten Überblick ist der Versuch 
gemacht worden, ein vollständiges und anschauliches 
Bild der Beugung und Streuung elektromagnetischer 
Wellen am leitenden Kreiszylinder zu bringen. Dar- 
über hinaus sind auch die Mittel geschaffen worden, ein 
tieferes Verständnis der Beugung als physikalische 
Erscheinung zu gewinnen. Weiteres auf diesem Gebiet 
ist in der Literatur zu finden [6]. 
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(1957). — [6] Kına, R., and T.T. Wu: Scattering and Dif- 
fraction of Waves. Cambridge: Harvard University Press 1959. 
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Reimer, L.: Elektronenmikroskopische Untersuchungs- und 
arationsmethoden. Berlin-Göttingen-Heidelberg: Springer 
. 300 8. u. 135 Abb. Gzl. DM 58.—. 

Das Buch ist für den elektronenmikroskopisch arbeitenden 
wenschaftler bestimmt, dem es einen Überblick über die 
chiedenen Untersuchungs- und Präparationsverfahren der 
tronenmikroskopie geben will. Es ist kein Lehr- oder 
»ptbuch für Anfänger, das alle Kniffe und Handgriffe 
hreibt — dafür ist am Ende jedes Abschnitts umfang- 
1e Spezialliteratur angeführt, aus der Einzelheiten bei Be- 
ersehen werden können. 

In logischem Aufbau wird eingangs die Bildentstehung im 
ttronenmikroskop und: die Elektronenbeugung behandelt, 
deren Eigenarten Präparationsverfahren und Bilddeutung 
eren müssen. Hier ist besonders die klare Darstellung der 
trastentstehung — ohne den sonst üblichen mathemati- 
n Ballast — hervorzuheben. Bei den Präparations- 
hoden hat der Verf. es u.a. verstanden, die umfangreichen 
ıtions- und Einbettungsverfahren histologischer Objekte 
st ihren Vor- und Nachteilen schön herauszuarbeiten und 
it auch dem Physiker einen guten Überblick über dies 
diskutierte Spezialgebiet zu vermitteln. Die vielfachen 
weise auf Artefaktgefahren und Fehldeutungen sind im 
blick auf die immer noch anzutreffende unkritische Ein- 
lung gegenüber elektronenmikroskopischen Untersuchun- 
besonders zu begrüßen. 

Ein Bezugsquellenverzeichnis bildet den Abschluß des 
tes, in dem man leider einen Abschnitt über die photo- 


phische Verarbeitung und die Auswertung der Aufnahmen - 


mißt. Im Gegensatz zu den guten Textabbildungen und 


den schönen Bildtafeln, die alle typischen Erscheinungen 
illustrieren, steht nur die unübersichtliche, mit zahllosen De- 
tails überladene Schnittzeichnung des Siemens-Mikroskops. 
Eine kurze tabellarische Zusammenstellung der kommerziellen 
Elektronenmikroskope mit ihren charakteristischen Daten 
(inel. Objektverschiebungsbereich!) könnte eine Neuauflage 
dieses nützlichen Buches noch bereichern. E. KınDEr 


Holm, R.: Electrie Contacts Handbook (Elektrische Kon- 
takte), 3. Aufl. Berlin-Göttingen-Heidelberg: Springer 1958. 
522 8. u. 194 Abb. Geb. 52,50 DM. 


Das Buch ist entstanden aus der 1941 erschienenen Mono- 
sraphie Horms „Die technische Physik der elektrischen Kon- 
takte‘‘ der Sammlung „Einzeldarstellungen auf dem Gebiete 
der Technischen Physik‘“ (Springer-Verlag). 1946 erschien 
eine zweite Auflage dieser Monographie unter dem Titel 
„Electrie Conctacts““ bei dem Verlag H. Geber, Stockholm. 
Sie wurde in dieser Zeitschrift, Bd. II, 1950, S. 230 besprochen. 

Wie Verfasser im Vorwort betont, war die Abfassung seines 
Buches nur möglich durch ständige Verbindung mit den ein- 
schlägigen Firmen, früher besonders mit dem Siemens-Konzern, 
jetzt auch mit der Stackpole-Carbon-Company of St. Marys, 
Pa, USA, wo der Verfasser seit Jahren tätig ist. Ferner fand 
Horn auch jetzt wieder eine treue Mitarbeiterin in seiner Frau 
Erse Horm. — Die jetzige Auflage ist fast völlig neu ge- 
schrieben und in vieler Hinsicht erweitert. Dabei sind auch 
besonders neue Untersuchungen des Verfassers und seiner 
Mitarbeiter benutzt. Jedoch ist die Disposition und der Um- 
fang des gesamten behandelten Stoffes etwa derselbe wie bei 
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der 2. Auflage geblieben. — Zweifellos ist das Holmsche Buch 
immer mehr zu einem Standardwerk bezüglich der elektrischen 
Kontakte geworden. Daß die 3. Auflage in englischer Sprache 
geschrieben ist, dürfte seiner Benutzung in Deutschland 
kaum Abbruch tun. W. MEISSNER 


Stasiw, 0.: Elektronen- und Ionenprozesse in Ionen- 
kristallen. Struktur und Eigenschaften der Materie. Bd. XXII. 
Hrsg. von S. Frücee. Berlin-Göttingen-Heidelberg: Springer 
1959. 307 S. u. 107 Abb. Geb. DM 66.—. 


In den drei Jahrzehnten, die seit den grundlegenden 
Untersuchungen von HırscH und PoHL vergangen sind, ist 
eine große Zahl von Veröffentlichungen über optische Unter- 
suchungen an lonenkristallen erschienen. Eine Übersicht 
über dieses Gebiet fehlte um so mehr, als sich mit zunehmender 
Zahl der Arbeiten zeigte, daß die Erscheinungen selbst in 
den einfachsten Ionenkristallen recht verwickelt und viel- 
fältig sind. Zudem sind verschiedene Veröffentlichungen auch 
verhältnismäßig schwer zugänglich. Es ist daher sehr zu 
begrüßen, daß nun eine geschlossene Darstellung über Elek- 
tronen- und Ionenprozesse in Ionenkristallen entstanden ist, 
in der photochemische Prozesse besondere Berücksichtigung 
finden. Daß das Buch von Stasıw geschrieben wurde, dessen 
Forschungsgruppe wir einen großen Teil der neueren Ergeb- 
nisse verdanken, kommt den Ausführungen sehr zustatten. 
Auf Grund seiner umfassenden Kenntnis der Materie ist 
der Autor namentlich auch in der Lage, Vorstellungen, die 
ihm noch nicht völlig gesichert erscheinen, kritisch zu kom- 
mentieren. Gerade dies ist für Leser, die das Buch zur Ein- 
arbeitung benützen, sehr wertvoll. Auf den reichen Inhalt 
kann hier nicht ausführlich eingegangen werden. Es sei nur 
darauf hingewiesen, daß die Alkalihalogenide recht umfassend 
abgehandelt werden. Bei den Silberhalogeniden stehen natur- 
gemäß die Untersuchungen der Stasiwschen Schule im Vorder- 
grund. Am Schluß finden sich wichtige Abschnitte über 
Anwendung von Kernresonanz und adiabatischer Näherung 
auf Kristalle mit Störstellen. 

Nach Durchsicht des Buches drängt sich die Frage auf, 
wer mehr zu beglückwünschen ist, der Autor oder die Leser, 
die nun einen ausführlichen Leitfaden als Führer durch die 
Störstellenphysik und Photochemie der Ionenkristalle in die 
Hand bekommen haben. W. WAIDELICH 


Handbuch der Physik. Hrsg. von S. FLüsge. Bd. XXXIV: 
Korpuskeln und Strahlung in Materie II. Berlin-Göttingen- 
Heidelberg: Springer 1958. VIII, 316 S. u. 213 Abb. Gzl. 
DM 78.—. 


Dieser Band 34 des Handbuchs der Physik bringt Beiträge 
von R. KorzAarH (Mainz) über den „Durchgang langsamer 
Elektronen und Ionen durch Gase“, von R.D. BIRKHOFF 
(Oak Ridge) über ‚The passage of fast electrons through 
matter‘‘, von L. Sımons (Helsinki) über ‚‚Positronium‘“, von 
E. MERZBACHER und H.W. Lewis (Chapel Hill bzw. Durham) 
über „X-ray production by heavy charged particles‘‘, von 
W. WnarınG (Pasadena) über „The energy loss of charged 
particles in matter‘ und von R.D. Evans über „Compton 
Effect‘. Der erste Artikel beschränkt sich auf Bestimmungen 
der gesamten Wirkungsquerschnitte und der elastischen Streu- 
ung von langsamen Elektronen (bis etwa 100eV) und auf 
Wirkungsquerschnitte, Streuungen und Umladungen von 
langsamen Ionen (bis einige keV). Der umfassende Artikel 
von BIRKHOFF enthält unter anderem eine eingehende Dis- 
kussion der diskreten Energieverluste und ihrer teilweisen 
Deutung durch die Plasmatheorie. Der Beitrag von Sımons 
über das Positronium ist vorwiegend vom experimentellen 
Standpunkt aus geschrieben. MERZBACHER und Lewıs be- 
fassen sich in erster Linie mit der Theorie der K- und L- 
Schalen-Ionisation und mit experimentellen Ergebnissen für 
Protonen und «-Teilchen bis zu Energien von 5 MeV. W. Wna- 
LING gibt eine Zusammenstellung experimenteller Daten der 
Energieverluste und Reichweiten in verschiedensten Sub- 
stanzen, insbesondere für Protonen und «-Teilchen von we- 
niger als 10 MeV Energie. R.D. Evans vereinigt in seinem 
Artikel theoretische und experimentelle Ergebnisse über den 
Compton-Effekt in gleich umfassender Weise. Alles in allem 


vermittelt auch dieser Band des Handbuches wieder 
Fülle wertvoller Erkenntnisse in meist vorbildlich k 
trierter Weise. H. EwA 


Angerer, E. v.: Technische Kunstgriffe bei physikali 
Untersuchungen. Zwölfte Auflage. Hrsg. von H. EBERT 
Mitarbeit zahlreicher Fachwissenschaftler. Brauns: 
F. Vieweg 1959. 464 8. u. 172 Abb. Leinen DM 24 
Wieder liegt eine neue Auflage der längst zum 
Laboratoriumsbestand gehörenden Kunstgriffe vor. 
einzigartige Werk, das ursprünglich Anfängern eine 
zum Experimentieren geben sollte, wie sich E. v. An@ 
bescheiden ausdrückte, war schon von der 1924 erschi 
ersten Auflage an ein wegen seiner Originalität und Zi 
lässigkeit geschätzter Ratgeber. Oberster Grundsatz E 
GERERS war es, nichts aufzunehmen, was er nicht - 
probt hatte. Es liegt auf der Hand, daß sich heute di 
nur durch eine Gemeinschaftsarbeit einer größeren 
Experten erreichen läßt. Auf dieser Basis verstandeesH. F 
vorzüglich, die Kunstgriffe dem heutigen Stand der P} 
anzupassen und ein Werk vorzulegen, das den Bedü 
des auf den verschiedensten Gebieten tätigen Experi m 
physikers gerecht wird. 
Daß die Einheitlichkeit des Buches unter der F 
mehr als 25 Referenten teilweise etwas leidet, ist uner! 
gegenüber dem Vorzug absoluter Zuverlässigkeit. Zur 
für alle Kenner der noch zu Lebzeiten E. v. ANGER 
Monographie erschienenen Auflagen findet sich an ve 
nen Stellen noch wörtlich die ursprüngliche Fassun, 
sich der Autor nicht nur als Meister der Experimentie 
sondern auch der Darstellung zu erkennen gibt. Dies 
daß auf einigen Gebieten die Kunstgriffe noch dieselbe 
blieben sind (z.B. Fadenkreuz aus Spinnenkokon! 
Juliussche Aufhängung). Viele neue Kapitel konnten da 
erst jetzt von H. EBERT aufgenommen werden (etwa Sei 
tionszähler oder Transistor). / 
Beim Durchlesen staunt man immer wieder über die’ 
des dargebotenen Materials. Ein umfangreiches Schlag 
verzeichnis erleichtert das Arbeiten mit dem Werk, das k 
besonderen Empfehlung mehr bedarf. W. WAIDELI 


Sauter, F.: Differentialgleichungen der Physik. Samm 
Göschen Bd. 1070. Berlin: Walter de Gruyter 1958. 
u. 16 Abb. DM 2.40. 

Das Bändchen vermittelt eine Einführung in die Beh 
lung der hauptsächlich in der Physik vorkommenden 1 
rentialgleichungen: gewöhnliche Differentialgleichung 
Mechanik und partielle Differentialgleichungen der 
physik (schwingende Saite, Wärmeleitungsproblem 
Wellengleichung und die Potentialgleichung). Die Eige 
ten der Funktionssysteme, die durch die untersuchten B 
wertprobleme definiert werden, Hermitesche und Laguerre 
Polynome, Kugel- und Zylinderfunktionen, werden bespro 
Den Abschluß bildet eine Darstellung der Störungsrech 

Die neue Auflage unterscheidet sich von den vorhe 
den durch eine Vereinfachung einiger Beweise. Das Bä 
wird weiterhin für all jene, die sich von den Kenntn 
Differential- und Integralrechnung ausgehend in die B 
lung der Differentialgleichungen der Physik einarbeiten 
ein leicht verständlicher Führer sein. E.F 


Herforth, Liselott und Herbert M. Winter: Ultras 
Grundlagen und Anwendungen in Physik, Technik, In 
Biologie und Medizin. Leipzig: B.G. Teubner 1958. 
u. 65 Abb. Geb. 12.— DM. 

Gegenüber dem umfangreichen Handbuch BEr@ 
über Ultraschall will das vorliegende Bändchen einen 
Überblick über die Grundlagen und Anwendungsmöglie] 
des Ultraschalls geben. Die behandelten Kapitel sind 
zeugung von Ultraschallwellen. Ultraschallfeld und 
gen im Ultraschallfeld. Nachweis und Messung des 
schalls. Anwendung in Physik, Technik und In 
Ultraschall in Biologie sowie in Medizin. — Die vielen 
Abbildungen erleichtern das Lesen. Besonders geeignet d 
das Buch zur Information über die Anwendung des [ 
schalls sein. 


